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The Reflecting and Resolving Power of Calcite for X-rays 


By Samve Kk. ALLISON 
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(Received May 20, 1932) 


With a precision double spectrometer, the rocking curves from three pairs of 
calcites in the (1, —1) positive have been observed for seven wave-lengths from 
O.21A, (WoKa;) to 2.3A, (Cr Ka,). Three quantities were observed at each wave- 
length, (1) the percent reflection (P)), defined as the ratio of the maximum ionization 
current obtainable from crystal B to the ionization current produced by the x-rays 
incident on B from Al, (2) the coefficient of reflection (2), defined as the ratioof the area 
of the rocking curve from B to the intensity incident on B from A, (3) the half width 
at half maximum of the rocking curve from B. One of the pairs (crystals II1) was obvi- 
ously more perfect than the others, and its values of Po, XR, and w are compared with 
theory. It is shown that none of the three quantities defined above is simply related 
to the analogous observations on a single crystal with monochromatic, parallel, inci- 
dent x-rays, but that numerical predictions of the values of the above quantities 
may be obtained from somewhat laborious calculations based on Darwin's theory of 
diffraction by a perfect crystal as modified by Prins to take account of absorption. 
It is found that the predicted and observed coefficients of reflection are in good agree- 
ment in the wave-length range between 0.5 and 2.3A. On the basis of this agreement 
it is possible to state that calcite surfaces may be obtained for which there is no evi- 
dence of mosaic structure from measurements of the coetticient of reflection by double 
spectrometer methods in the wave-length region between 0.5 and 2.3A, 


EXPERIMENTAL PART 
Introduction 


HE experiments reported here are similar to those performed on calcite 

in previous researches by other investigators.' The object of the experi- 
ments in the present case was primarily to provide a set of calibrated crystals 
to use in intensity measurements in this laboratory. The method consists 
essentially in comparing the intensity reflected from crystal B of a double 
spectrometer with that incident on & from crystal A, the instrument being 
set in a parallel position. In the present experiments only the (1, —1) posi- 
tion in the notation of Allison and Williams! was used. 


' A. H. Compton, Phys. Rev. 10, 95 (1917); Davis and Stempel, Phys. Rev. 17, 608 (1921) 
and Phil. Mag. 45, 463 (1923); Wagner and Kuhlenkampff, Ann. d. Physik 68, 369 (1922): 
Allison and Williams, Phys. Rev. 35, 1476 (1930); Davis and Purks, Phys. Rev. 34, 181 (1929) 
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Apparatus 


The double spectrometer used in these experiments is shown in Fig. 1. It 
was constructed by the Gaertner Scientific Corporation, Chicago, for the 
author. The principles of construction of the instrument are those previously 
described by Williams and Allison” The advantages of this design are prin- 
cipally due to the fact that axis B can be rotated about axis <1, and the ioniza- 
tion chamber can be rotated about axis B. This enables the instrument to 
operate both in the parallel and anti-parallel positions, and allows settings 
for widely different wave-lengths to be made without moving the x-ray tube. 
In a design for a double spectrometer recently described by A. H. Compton,’ 
the x-ray tube must be rotated about axis <1 to change from one wave-length 








Fig. 1. The double spectrometer. The electrometer (above) is wrapped with wadding to 
protect it from sudden temperature changes. The instrument stands on a heavy stone slab 
mounted on concrete piers. The distance from the top of this slab to the shelf supporting the 
electrometer is 54 cm. Behind the instrument is a lead covered room in which the x-ray tubes 
are operated. 


to another, which may be a considerable disadvantage. The author's instru- 
‘ment, however, is not designed for accurate absolute wave-length determina- 
tions, since angular measurements of rotation through large angles about axis 
B cannot be made with great accuracy. The present instrument is best 
adapted to the investigation of rocking curves, fine-structures, etc. 

Although the instrument of Fig. 1 does not differ in principle from the 
previous one of Williams and Allison, much more care was taken in its con- 
struction. Because of the wide range of positions that axis B and the ioniza- 
tion chamber may assume with respect to axis A the problem of properly 


? Williams and Allison, J.O.S.A. and R.S.1. 18, 473 (1929). 
* A. H. Compton, Rev. of Sci. Instr. 2, 365 (1931) 
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balancing the weight on the various members is complicated. This problem 
was much more successfully handled in the present instrument than in its 
prototype. In particular, the essential requisite of the instrument namely 
that axis A be parallel to axis B was carefully worked out by scraping the 
face plates which connect the axes until as tested by a sensitive spirit level 
the two were parallel to less than one minute of arc. 

The crystal mountings are somewhat more elaborate than those of the 
previous instrument in that two controlled horizontal translations are possi- 
ble. In addition to the motion perpendicular to the reflecting face of the cry- 
stal, it is possible to shift the entire mounting parallel to the reflecting face. 
In this way the x-ray beam may be shifted along the surface of the crystal 
without moving that surface out of the axis of rotation. An attempt to design 
a plate for mounting the crystals which could readily be removed from the 
instrument with the crystal on it, stored in a dessicator, and then replaced, 
was unsuccessful. No scheme was devised such that the plate would return to 
its position on the axis to the required degree of accuracy (less than 10 seconds 
of arc) in a series of trials. The entire mounting had to be slid out of the 
grooved ways provided for the translatory motion in order to remove the 
crystals from the instrument. 

Some distances on the instrument are: between the centers of the two col- 
limating slits, 14.5 cm; from a point halfway between the two slits to axis 
A, 13.2 cm; from axis A to axis B, 13.8 cm; from axis B to the front window 
of the ionization chamber, 5.8 cm. 

The ionization chamber is of the type mentioned in the previous paper by 
Williams and Allison,? and in the present experiments was filled with methyl! 
bromide at one atmosphere pressure. The electrometer was of the Compton 
type, as manufactured by the Cambridge Instrument Company. It was 
operated at a sensitivity of 1.84 meters per volt as read on a scale 1.75 meters 
distant from the mirror of the instrument. It was shown by trial that a linear 
relationship between scale deflection and voltage on the unearthed quadrants 
existed at all points of the scale. 

Commercial x-ray tubes were used for studies with the tungsten and chro- 
mium K series. The tube with the chromium target had a special glass window 
to transmit the radiation. The other wave-lengths used were obtained from 
tubes of the type previously described by the author.‘ The high voltage for 
operating the tubes was obtained from the generating plant described by Alli- 
son and Andrew.® In the experiments in which the Ka; lines of Mo, Cu, Fe, 
Ag, and the Z8, line of Ir were used, the x-ray tubes were run at steady 
voltages between 30 and 42 k.v., higher values being used for Mo and Ag. 
Experiments with the Cr Ka line were conducted at 38.4 k.v. and 26.8 k.v. 
With this radiation (wave-length 2.28A), about 86 percent was absorbed in 
the 58 cm air path alone, according to the absorption coefficient found by 
Spencer.® This raises the possibility that if the voltage on the tube is suffi- 
‘ Allison, Phys. Rev. 30, 245 (1927). 

§ Allison and Andrew, Phys. Rev. 38, 441 (1931). 
® Spencer, Phys. Rev. 38, 1932 (1931). 
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ciently high, the radiation reflected at the glancing angle for Cr Ka; may con- 
tain little of this wave-length and be mostly second and higher orders of the 
general radiation. It was shown that this effect did not invalidate the results 
for chromium by using the instrument as a single crystal spectrometer in the 
usual way and examining the spectrum, and also by the fact that the rocking 
curves at the two voltages indicated above for Cr were in agreement. The 
tungsten A spectrum was obtained from a tube operated at 110 k.v. 


Adjustments 


The instrument was put in adjustment by the procedure previously de- 
scribed by Allison and Williams! with some minor refinements in technique. 
One large improvement over the methods of Allison and Williams resulted 
from the fact that it was found possible to use a telescope with a Gauss eve- 
piece in the adjustment of the crystals. This greatly increased the accuracy 
with which the cleavage faces could be set parallel to the axes of rotation. 
With the crystal holder removed, a piece of plane parallel glass, mounted on a 
small block equipped with levelling screws, was placed over the axis. The 
plane parallel plate and the telescope were adjusted until at two positions of 
the plate 180° apart the reflected cross-hairs coincided with those in the 
telescope. The axis of the telescope was then perpendicular to that of the 
spectrometer. The plane parallel plate was then removed, the crystal in its 
mounting attached to the axis, and the crystal tilted until the reflected cross- 
hairs from its face coincided with those in the telescope. In this way it was 
found possible to set the crystals accurately enough to obtain the minimum 
rocking curve width in parallel positions at once, and the tilting process 
previously described by Allison and Williams was not necessary. The tele- 
scope was placed roughly 2 meters from the spectrometer, and the crystals 
could be set with their faces within +15 seconds of arc of the axes of rotation 
by this method. 

In agreement with previous observers the rocking curves were found to be 
independent of the width of the slits used to limit the horizontal divergence 
of the beam of x-rays. The same slit widths, however, were not used at all 
wave-lengths. The inadvisability of such a procedure may be realised by con- 
sidering the area on the crystal face irradiated by the x-ray beam at different 
glancing angles. For a constant vertical divergence of the beam this area 
is proportional to its horizontal extent, which to a sufficient approximation 
is given by 

s = 2L.a/L sin @ (1) 


where L.. is the distance from a point midway between the slits to the axis of 
the crystal, a is the slit width (assumed the same for both slits), L is the dis- 
tance between the slits, and @ is the glancing angle. In the experiments on the 
tungsten K radiation, sin? =0.034, and the irradiated area becomes large 
relative to that used for softer wave-lengths. The range of wave-lengths 
present in the beam after reflection from crystal A can be calculated from 
the slit width and the distance between the slits by the formula 


dy = d\ cot 6d0 = (2ad cot 6)/L. 
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The various values of dX used are shown in column 4 of Table I. The vertical 
divergence of the beam was limited by a slit over the front window of the 
ionization chamber, at a distance of 5.8 cm from axis B. The other limitation 
was the finite size of the focal spot, certainly not greater than 0.8 cm in dia- 


TABLE I. Adjustments. 








Radiation Crystals a dy h Area irradiated 
on crystal B. 


6.2X.U. 0.8cm 1.0 cm? 


W Ka I, V 0.012cm 8 
“ il 0.012 6.2 0.5 0.65 
Ag Ka: 1 0.030 26 0.2 0.24 
“ HI 0.040 34 0.2 0.32 
“ V 0.015 13 0.8 0.48 
Mo Ka 1] 0.030 25 0.2 0.19 
" II 0.040 33 0.2 0.25 
“ V 0.015 12 0.8 0.42 
Ir Li II 0.030 25 0.2 0.12 
“ Il 0.040 33 0.2 0.16 
“ V 0.015 12 0.8 0.23 
CuKa II 0.030 24 0.2 0.09 
“ III 0.040 32 0.2 0.12 
“ V 0.020 16 0.8 0.23 
Fe Ka II 0.030 24 0.2 0.07 
“ Ill 0.040 32 0.2 0.09 
“ V 0.020 16 0.8 0.18 
0.030 23 0.8 0.24 


Cr Ka, H,1I,V 








} 
| 
} 
| 
| 
| 








eter and 58 cm distant from this slit. The width of the slit in front of the 
ionization chamber limiting the vertical divergence of the beam in various 
experiments is given as /: in column 5 of Table I. In column 6 of the same ta- 
ble is given the area irradiated on crystal B in cm’, obtained roughly by multi- 
plying the values in column 5 by the appropriate value of s obtained from 
Eq. (1). 


Crystals 


Five pairs of split calcites were examined in this work. Two of the pairs, 
calcites I and IV, were examined only for a few wave-lengths and the data are 
not reported here. The three pairs II, III, and V had freshly cleaved surfaces. 
They were clear calcites, suitable for the construction of Nicol prisms, but of 
unknown origin. They were split by a method recommended to the author by 
Professor Bergen Davis. The sample to be split was mounted on the carriage 
of a milling machine, in such a way that the lateral motion of the carriage was 
accurately parallel to the cleavage direction. A very thin circular saw was 
then run at slow speed and a groove made across the specimen. By continued 
sawing in this groove a split soon started and the slight jars due to the impact 
of the saw teeth on the crystal spread it until fissure was complete. The two 
halves of the crystal were then immediately placed in a dessicator and kept 
until used. The crystals were parallelopipeds usually from 0.5 to 1 cm thick 
and some 5 cm on a side. 

The ability of the crystal surfaces to give a sharp reflected image of the 
cross hairs of the Gauss eye-piece was different in the various specimens. It 
is a rather disconcerting fact that crystals III], which proved to give the high- 








6 SAMUEL K. ALLISON 
est resolving power for x-rays, were the poorest optical reflectors. Crystals II 
gave a very clear reflected image, as good as that from a polished glass plate. 
Crystals V were intermediate between IT and IIT in this respect. 

The most extensive investigations were carried out on crystals III], and 
before these crystals were finally removed from the instrument a test was 
performed to make sure that the resolving power had not changed during the 
experiments. On Nov. 19, 1931 these crystals gave a rocking curve of 2.6” 
half width at half maximum for Mo Ka,(1, —1). After the investigation of 
the other wave-lengths was completed, on Jan. 2, 1932, the experiment with 
Mo Ka, was repeated and again 2.6” obtained. There is therefore no evidence 
that the crystals II] had changed during the period of experimentation. 


Measurements 


From the rocking curves for a given wave-length three quantities were 
determined. The half width at half maximum, w, is the angular range in sec- 
onds of are through which crystal B must be turned to reduce the power en- 
tering the ionization chamber to one-half its maximum value. By maximum 
value is meant the largest ionization current observed, corrected for base-line. 
This base-line was taken as the ionization current observed when crystal B 
was rotated from its angular position for maximum reflection through an 
angular displacement of 6 or more times the half width at half maximum. The 
percent reflection, Po, is the ratio of the maximum ionization current obtain- 
able by reflection from crystal B to that obtainable from the beam incident 
on B from A. The coefficient of reflection, R, is in this paper defined as the 
ratio of the area under the rocking curve to the ionization current obtainable 
from the beam incident on crystal B from crystal A. Values of R observed and 
calculated in this paper are based on the radian as the unit of angle. 

A regular routine was developed for the recording of rocking curves, the 
steps of which will now be described. A run was begun with crystal B out of 
the x-ray beam. The beam reflected from crystal A (which had been set at the 
proper angle) was allowed to enter the ionization chamber, and the current 
through the x-ray tube was adjusted until the rate of deflection of the electro- 
meter was between 5 and 10 scale divisions per second. When conditions had 
become steady, a series of readings, usually three, were made on the beam 
reflected from crystal A. Between these readings the ionization chamber was 
shifted slightly to make perfectly sure that the entire beam was entering the 
window. If these three readings agreed to within 1 percent, the ionization 
chamber was displaced 5° and a reading taken of the diffusely scattered radia- 
tion. The difference between the average of the three readings and the base- 
line reading was taken as proportional to the intensity incident on B from A. 
Except in experiments on W Ka, the base-line was negligible here. Crystal 
B was then advanced into the beam, and the ionization chamber placed to 
receive the reflected beam from it. The rocking curve was then taken. Care 
was taken to obtain readings far enough from the maximum to supply a reli- 
able base. Crystal B was then removed, and the power incident upon it again 
determined; the average of the initial and final values being used as the inten- 
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sity incident on B during the taking of the rocking curve. The results are 
given in Table II. Each value is the result of at least two independent obser- 
vations. 

TABLE II. Experimental results. 


w=half width at half maximum of rocking curve (seconds). 
P,\=percent reflection. 
=C 


& 


R=coefficient of reflection, in (radians). 
, — | uw 100 P, R < 105 

we (2) in wviin @ vin Wl OV 
W Ka 208.6 | 4.2 2.3 7.2 33 35 28 2.34 1.17 2.14 
AgKa 558.3 | 3.8 2.2 7.0 48 59 34 2.15 1.53 2.38 
MoKa 707.8 | 3.9 2.6 7.2 52 63 33 2.31 1.86 2.44 
IrLe, 1155 | 48 4.1 8.1 59 64 37 3.49 2.95 3.24 
Cu Ka, 1537 5.6 4.9 7.8 61 62 46 3.92 3.46 3.91 
FeKa, 1932 | 6.8 6.2 9.9 54 58 42 4.45 4.18 4.45 

2285 | 7.7 7.1 9.3 52 55 42 14.73 4.68 4.37 


CrKa 

If we examine the data on w in Table II, we see that the w values for crys- 
tals III are consistently lower than for crystals II or V. This is taken to indi- 
cate that crystals III approach the ideal of a perfect crystal more nearly than 
do the others investigated, since any mosaic structure would enlarge the an- 
gular range over which reflected can take place. The remainder of this report 
will be concerned with a comparison of the values of w, Po, and R observed 
for crystals III with the predictions of the theory of diffraction by a perfect 
crystal. 


THEORETICAL PART 
Theory of diffraction of x-rays by an absorbing, perfect crystal 


The theory of the diffraction of x-rays from the face of a perfect crystal 
has been developed by Darwin,’ and later, Ewald.* In the discussion here we 
shall consider the theory as derived by Darwin. In Darwin's derivation, the 
absorption by incoherent processes of the x-rays as they penetrate the crystal, 
is neglected. The absorption due to ejection of photoelectrons or recoil elec- 
trons would be considered as an incoherent process in this sense. The result 
of Darwin's calculations can be expressed in the following form. The ratio of 
the intensity reflected from the crystal to that incident upon it is 


T (3-0) | = q [° ( 3) 
—_ = |— . 3 
Io v 4- (v? — q*)' 2 
where 

 ,= } 2rd(B — 0) cos 6} A, (4) 6 = 05 + 6 sec Ay cosec Io, (6) 

2rd F208) . mr . 

qa —-——— j ——» (5) 0) = sin~' — - (7) 

Asin@ F(O) 2d 


7 Darwin, Phil. Mag. 27, 325 and 675 (1914). 
8 Ewald, Phys. Zeits. 26, 29 (1925). 
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In the above expressions, (8—8@) is the deviation of the glancing angle from 
the angle 6, which is the corrected Bragg angle. \ is the wave-length of the 
x-rays in air; d, the grating space of the crystal; 7, the order of the spectrum, 
and 6 the deviation of the index of refraction from unity. The sign to be used 
for the radical in the denominator of Eq. (3) is determined by the physical 
requirement that the ratio on the left hand side of this equation must be less 



































0— 
MoKa 
CukKa 
is 
Crka 
‘o- 
i. a 
0.4-— 
o2-— 
= | L f 
s) -2 “1 z i 2 3 


Fig. 2. Theoretical diffraction patterns from a single crystal. The outer curve is Darwin's 
theory, uncorrected for absorption, the inner curves illustrate Prins’ absorption correction, 
With values of £ plotted as abscissae, the curves apparently have about the same width at half 
maximum. If instead of £, we consider the actual angular deviations from @, the width increases 
with wave-length. 


than or equal to one. F(26,) is the ratio of the amplitude of the wave scat- 
tered at angle 26, from the unit cell to the amplitude incident upon it. F(0) is 
analogous ratio for the wave scattered in the forward direction of the incident 
beam. Let us first consider the simplified case in which the entire incident 
radiation is plane polarized in such a direction that the electric vector lies 
in a plane perpendicular to the plane of incidence of the radiation on the 
crystals. For the case, 

F(26))/F(0) = FZ—! (8) 


where F is the structure factor and Z the number of electrons in the unit cell 
of the crystal. In the range —q<v<gq, Eq. (3) gives rise to 100 percent reflec- 
tion. If we call this range of angle A#, we have 


A@, = 46FZ-' cosec 26. (9) 


In this and following equations, the subscript ¢ indicates that the electric 
vector of the radiation is perpendicular to the plane of incidence. We shall use 
this A@, as a unit of measure of angular deviation from the corrected Bragg 
angle 6, and introduce the quantity £, defined by 


f = (6 - 0)/A0,. (10) 























REFLECTING AND RESOLVING POWER OF CALCITE 


If this &, is introduced into Eq. (3) ,we obtain 
F(g.) = | 26 + (486% — 1)'| -* (11) 


A plot of F(£,) is shown as the outer curve of Fig. 2. 

Prins® has recently modified Darwin's theory in an attempt to take ac- 
count of the absorption of the x-rays in the crystal by incoherent processes, 
as mentioned above. It is well known that in the classical dispersion theory, 
the index of refraction is the real part of a complex number, which number we 
may represent by (1—4—78). The real part of this number, (1—4), is the 
ordinary index of refraction, i.e., the ratio of the phase velocity in vacuum to 
that in the medium. The imaginary part is related to the linear absorption 
coefficient u, of the radiation in the medium by 


B = Nui ‘der. (12) 


The essential feature of Prins’ method is that he has treated the absorption in 
the classical manner, as if it took place by coherent processes, which in gen- 
eral leads to the substitution of (6+78) for 6 when it appears in Darwin's 
treatment. 

We shall take the result of Prins’ calculations as it appears in Eq. (11) of 
his paper. This is, with a slight change of notation, 
T (3-0,) 


— (13) 


d+ ib 2 


~ | (B—6p) sin 8 cos 0)— (5+ 18) +( | (B—Bo) sin 8) cosBo— (6+ i8) \2_ (d+ ib)?)'2 





We shall first consider the application of this formula to a crystal composed 
of atoms of one kind only, in equally spaced planes. Under these conditions 
the quantity (d+7b) is related to the quantity (6+78) as the amplitude of the 
wave scattered by an atom at angle 26) to the amplitude scattered in the for- 
ward direction. In order to be more specific about the relationship between 
the two quantities, we will adopt the simpler of Prins’ two hypotheses, hy- 
pothesis 7. This states that all electrons in the atom may be treated in the 
same way as regards the amplitude scattered as a function of the scattering 
angle, irrespective of the level to which they belong. This hypothesis leads to 
the relation 


d + ib = (6 + i8)f(20)/f(0) (14) 


where f(26,) is the amplitude scattered by the atom at the angle 26, from an 
incident wave of unit amplitude. A more detailed and correct procedure would 
be to consider the K, L, M, etc. electrons in the atom separately, both with 
respect to their contributions to 6 and 6, and the angular distribution of 
the amplitude scattered from them. In this paper, however, the simplified 
procedure of Eq. (14) will be adopted. 


® Prins, Zeits. f. Physik 63, 477 (1930). 
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In a crystal such as calcite, which consists of interlaced planes of different 
kinds of atoms; instead of (d+7b) in Eq. (13), we must introduce a quantity 
fi( 260) 


D+ iB = Yi(6; + iB) et tmithrat tute, (15) 
F (0) 


In this expression, 6; is the contribution to 6 for the medium from atom 7. 
n is the order of the reflection; 4, &,/, the Miller indices of the reflecting plane, 
and x,;, y,;, 5; the coordinates of the atom 7 in the unit cell. The summation is 
to be carried out over all the atoms present in the unit cell. In more common 
notation, for the ¢ type of polarization detined above, 


fil 200) /f(0O) = FZ" (16) 


where F; is the atomic structure factor of the atom of type 7, and Z; is the 
number of electrons in the atom. 

In calculating the values of D+i7B by Eq. (15), the author has made a 
further approximation, namely, that the contributions to 6 and 8 from a given 
atom 7 are proportional to the fraction of the total number of electrons in the 
unit cell contained in the atom in question. If Z is the total number of elec- 
trons in the unit cell, this means that 


6; + i8; = Z(6 + 18), Z. (17) 


It is clear that more accurate calculations could be made, the best method 
being perhaps to subdivide the electrons in the unit cell into groups according 
to their binding energy, but in view of the approximate nature of the graphi- 
cal calculations to come later it would seem that the somewhat crude assump- 
tion of Eq. (17) is permissible. 

Substitution of Eqs. (17) and (16) in Eq. (15) gives 


D+ iB =Z7(6 + iB) SOF ett irit hata, (18) 


The summation in Eq. (18) is commonly known as the structure factor F, so 
that finally, 
D+ iB =Z"'¥F(6 + ig). (19) 


We are now ready to adapt Eq. (13) to a crystal of more than one kind of 
atom, and to put it into a form similar to Eq. (11). It should be noted that in 
Eq. (13) the reference angle is 69, the uncorrected Bragg angle, whereas in 
Iq. (10) the reference angle is @. In terms of 00, the expression for & in Eq. 
(10) becomes, after making use of Eq. (6). 

(8 — Oo) sin 09 cos 45 1 


ae a A (20) 
26FZ WZ 


Str 


a 


Putting D+7B from Eq. (19) in place of d+7) in Eq. (13), and introducing 
—, from Eq. (26), we obtain 
_ FZ—“(1 + 8/6) ‘ 
'G) * i a eet seems. (21) 
2E,F'Z-! — 18/6 + [(2¢.FZ-' — iB/6)? — F°Z-2°(1 + i8/6)?)}'? 
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If the effect of absorption is neglected, that is, 8 =0, this expression reduces to 
Eq. (11), Darwin’s original form. 

The calculation of the shape of the diffraction pattern to be expected 
from a perfect calcite crystal irradiated at the proper angles by monochromat- 
ic, plane, x-rays of wave-lengths 0.708A, (Mo Ka), 1.537A (Cu Ka), and 
2.285A (Cr Ka) has been carried out from Eq. (21). The values of 6 were ob- 
tained from Larsson.'® The values of 8 may be obtained from the linear ab- 
sorption coefficient by Eq. (12). The value of uw, was calculated for Mo Kay, 
in calcite from data collected by A. H. Compton," and extended to other 
wave-lengths by assuming the validity of the \*° law. The values used are 
shown in Table III. 

TABLE IIT. Values of 6 and 8 used in Eq. (21). 


A ri) Bp B/6 
0.708A 1.87 x 107* 1.31 <10- 0.007 
1.537 8.80 2.92 x107? 0.033 
2.285 19.5 1.43 «1076 


0.073 











The value of FZ~' for calcite in the first order from the cleavage faces, 
based on the atomic structure factors of Bragg and West" has been shown by 
Allison and Williams! to be 0.47. All the calculations of this paper have been 
based on this value, which may be somewhat low, as discussed later. Some of 
the values of F(£,) for these three cases, both from Eq. (11) (Darwin) and 
Eq. (21) (Prins) are shown in Table IV." 


TABLE IV Diffraction pattern of a perfect calcite crystal. 





F(é,) F (&,) Prins 








Eq. (10) 
to Darwin MoKa Cu Ka CrKa 
—2.13 0.014 0.014 0.014 0.014 
—1.06 0.063 0.062 0.062 0.062 
—0.85 0.106 0.105 0.105 0.105 
—0.64 0.234 0.233 0.232 0.229 
—0.43 1.000 0.959 0.874 0.709 
—0.21 1.000 0.973 0.881 0.758 
+0.00 1.000 0.971 0.870 0.735 
+0.21 1.000 0.960 0.830 0.670 
0.43 1.000 0.924 0.702 0.503 
0.53 0.501 0.486 0.426 0.341 
0.64 0.234 0.233 0.229 0.213 
0.85 0.106 0.105 0.105 0.103 
1.06 0.063 0.062 0.062 0.062 
2.13 0.014 0.014 0.014 0.014 


The calculated values of Table IV are shown graphically in Fig. 2. An interest- 
ing feature of the curves corrected for absorption by Prins’ method is that 


10 Larsson, Diss. Upsala 1929. Reported in Siegbahn, Spektroskopie der Rontgenstrahlen, 
2nd Ed., (1931) Table XIII, p. 39. 

1! A. H. Compton, X-rays and Electrons, (1926), Table VI-2 p. 182. 

2 Bragg and West, Zeits. f. Krist. 69, 118 (1928). 

13 A large part of the calculations involved in this paper were made for me by Mrs. G, S, 
Monk, computer of this laboratory. 
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they are asymmetrical. The possibility of detecting such a lack of symmetry 
by double spectrometer measurements will be discussed later. 

The areas under the curves of Fig. 2 are important in the theory of inten- 
sity of reflection. By integration of Eq. (11) it may be shown that 


| F(é,)d& = 4/3. (22) 


The unit of angular deviation in Eq. (22) is the angle A@,, as is shown by Eq. 
(10). If we adopt more fundamental units, i.e., radians, the area under the 
diffraction pattern is 446, 3. The areas under the other curves of Fig. 2 were 
obtained by use of the integraph, and are shown in Table V." 


TABLE V. Values of | F(&,)d&y from Eqs. (11) and (21). 


Line Darwin Prins 
Mo Ka; 1.333 1.250 
Cu Ka, 1.333 1.122 


CrKa 1.333 0.935 


Calculation of the shape of the (1, —1) rocking curve for plane polarized 
radiation 


A careful analysis of the operation of the double crystal spectrometer 
shows that the shape of the rocking curve obtained in parallel positions in the 
first order could not possibly be that of the curves in Fig. 2, even if perfect 
crystals were available and the instrument were in correct adjustment. 
Neither can be the area under a rocking curve be directly compared with the 
areas under the curves in Fig. 2. The fundamental physical reason for this is 
that all the theoretical curves previously described were calculated for paral- 
lel, monochromatic radiation. Due to the fact that the dispersion is zero in 
parallel positions, the spectral distribution of the radiation is not important, 
but, although the first crystal greatly reduces the horizontal divergence of 
the beam, it does not reduce it to zero, and the beam incident on crystal B 
is still a divergent one. The author has previously pointed out" that the func- 
tion governing the shape of the rocking curves is 


P(Ex.) = f F(E,)F (& me Epo)dé,. (23) 
In this expression, 
Exe = (Bn — On) / AO, (24) 


4 It is these areas which, expressed in radian measure, are usually defined as the coeffi- 
cients of reflection, and it is true that the reflecting power of a single crystal is proportional 
to them. It is one of the principal objects of this paper to make clear that there is a discrepancy 
between this definition and the experimental attempts to measure this quantity by the double 
spectrometer. In this paper, the term coefficient of reflection refers to area under the rocking 
curve of crystal B. This is quite a different thing physically from the areas in Table V. 

!& S. K. Allison, Phys. Rev. 38, 203 (1931) Eq. (8). 
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where (8 —6,) is the angular deviation of crystal B from the position in which 
its reflecting face is parallel to that of A. ®(£,,) is the intensity of the beam 
entering the ionization chamber from a single position of crystal B. F(£,) is the 
diffraction pattern for a single crystal as in Fig. 2. The derivation of Eq. 
(23) is perhaps made clearer by Fig. 3, in which the course of a ray through 
the instrument is shown. 





Fig. 3. Course of a ray through the double spectrometer in a parallel position, 
illustrating Eq. (23). 


It is the custom to express the intensity reflected from B in terms of the 
intensity of the beam incident on B from A, which is proportional to the area 
under the single crystal diffraction pattern. We thus obtain the following 
expression 


where P, is the ordinate of the rocking curve from crystal B according to the 
convention above. The percent reflection, Po, will be given by the value of P, 
when &,, =0, or 


| [F(é,) |?dé, 


Fp eras s (26) 


f F(é,)dé, 


Laue" has discussed the possibility of solving Eq. (23), that is, of deter- 
mining the shape of the single crystal diffraction pattern from the experi- 
mentally determined ®(£,;,), and has come to some important conclusions. 
In the first place, the observed rocking curve (®(£,.) or P,) will necessarily 
be an even function, whether or not the single crystal diffraction curves 
F(&,) are asymmetrical. If we assume that the single crystal diffraction pat- 
tern is an even function, then it is possible to find its shape by solution of Eq. 
(23), otherwise, not. This means that by experiments of the type reported 
here it is impossible to discover whether the single crystal diffraction patterns 
are really asymmetrical as in the theory of Prins. If asymmetrical rocking 
curves are experimentally found, the only legitimate interpretation is that the 
diffraction patterns of the two crystals used are not identical. It should be 
remembered that the preceding statement applies only to parallel positions. 


16 Laue, Zeits. f. Physik 72, 472 (1931). 
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In view of the findings of Laue, as reported above, the author has not 
attempted to solve the Eq. (23), using his observed ®(£,,), and discover the 
nature of the single crystal curve. The procedure adopted has been the reverse 
of this, namely to calculate ®(&,,) from Eq. (23), using theoretical F(é,)'s 
from Eqs. (11) and (21), and comparing the result with experiment. Calcu- 
lated values of ®(&,,) are shown in Table VI. These values were obtained 
mainly by graphical integration. 


TABLE VI. Calculated values of ®(&p,) from Eq. (23). The unit 
of angle is 40, from Eq. (9). 








ip Darwin Prins . 
SBo : Mo Ka, Cu Ka, Cr Ka, 
0 1.00608 0.970 0.765 0.521 
+(0.21 0.954 .876 694 475 
+ .43 .786 727 .566 401 
64 .590 552 .434 .308 
85 .396 .378 nus Bk. 
-1.06 . 269 .200 .161 .135 
+1.060 .077 .070 063 052 
+2.13 .040 .036 .033 .029 
+2.66 O24 .022 .020 
+3.19 .019 014 .013 O12 








Calculation of the shape of the rocking curve with unpolarized incident ra- 
diation 

Before comparing the calculations summarized in Table VI directly with 
experiment, we must consider the fact that after reflection from crystal A, 
the beam incident on crystal B is partially polarized. Since the radiation used 
was in each case a characteristic emission line of the target, it is known that 
the radiation incident on crystal A was completely unpolarized. We may then 
arbitrarily divide up the beam incident on A into two components of equal 
intensity, one polarized so that its electric vector is perpendicular to the plane 
of incidence (subscript 7) and one with its electric vector parallel to the plane 
of incidence (subscript 7). We then treat each component separately in its 
passage through the apparatus. In Darwin's treatment, the extent of the 
region of 100 percent reflection with z-polarization is 


Ad, = 46FZ—' cosec 265 cos 209 = 46/Z-! cot 265 (27) 
which is to be compared with Eq. (9). We now introduce the quantities 
£, = (8 — 0)/AO, = &, sec 20 (28) 
and 
Expr = (Bn — On)/AOe = Ene sec 20. (29) 


By use of the preceding three equations, the derivation of the rocking 
curves for radiation polarized so that its electric vector lies in the plane of 
incidence (z-polarization) may be accomplished simply by changing the 
subscript ¢ to 7 in all the equations in which it appears up to and including 
Eq. (26). In adding the contributions from ¢- and 7-components to obtain the 
unpolarized values, it must be remembered that the angular unit in which ¢ 
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and w quantities are expressed is different. (Compare Eqs. (27) and (9).) 
We shall adopt the angular unit A@, in the following equations dealing with 
unpolarized radiation, and change quantities expressed in terms of A@, to 
their values in terms of A@,. We then obtain 

P(Ez.) + P(Enx) COS 245 P(Ez,) + (cos 269) (En. sec 260) 








; - : (30) 
'] F(&,)d&, + J F(é,)dé, (1 + cos 26) J F(&,)d&. 


where P is the ordinate of the rocking curve, that is, the ratio of the ionization 
current obtained from a setting of crystal B to the current obtained from the 


cS 
ened 


0.4 p 





Or 


2 





G,-9, (seconds) 











-40 -30 -20 -10 Oo 10 20 30 40 


Fig. 4. Comparison of observed and calculated rocking curves for Cr Ka. The dotted curve 
is the rocking curve calculated from Darwin's theory, in which absorption in the crystal is 
neglected. The solid curve is the curve derived from Prins’ theory. The circles represent ex- 
perimental points. The points are not fitted to the curve in any way, for instance by an ad- 
justable constant. 


radiation incident on 5 from A. If £g,=0in Eq. (30), that is, the two crystals 
are accurately parallel, it reduces to Eq. (26), which means that the percent 
reflection, Po, is independent of the polarization. This, of course, is not true 
of the coefficient of reflection, R. 

The use of Eq. (30) will be illustrated by a sample calculation of a point on 
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the rocking curve for Cr Ka; (1, —1) according to Prins. The denominator of 
Eq. (30) can be obtained from a knowledge of the Bragg angle for Cr Ka, 
on calcite, with the values in Table V. Its value is (1+0.715) (0.935) = 1.603. 
Let us consider a setting of crystal B such that &,,=0.64. From Table VI we 
find &( &;,) =0.308. The value of &,, sec 26) is 0.64 XK 1.40 =0.896. From inter- 
polation of Table VI, we find ©(0.896) =0.194, and ©(0.896)cos 2@ is 0.139. 
Thus P for this setting of the crystal is 0.279. 

The calculated values of A@, from Eq. (9) are shown in Table VIII. They 
are computed for FZ !=0.47, using 6-values from Table III, with the excep- 
tion of W Ka, where 6 is extrapolated by assuming a variation with A*. The 
rocking curves to be expected from calcite for Cr Ka; (1, —1) with unpolarized 
radiation from Darwin's and Prins equations can be plotted from the calcu- 
lated values shown in Table VII. 


Tape VIL. Theoretical form of the calcite rocking curve Cr Kay, (1, —1). 








. 8B—OB P 

she (seconds) Darwin Prins 

0.0 0 0.798 0.557 

wae 2.27 .698 496 

43 4.064 .540 396 

64 6.91 .370 .279 

.85 9.18 .229 .183 

+ 1.06 11.4 125 .116 

+1.00 17.3 .045 045 

+2.13 23.0 .024 .025 

+ 2.66 28.7 O15 .016 

+ 3.19 34.5 O11 010 

TaBLe VIII. Values of A0,. 
AM, Ad, 
Line r cosec 20 radians seconds 

W Ka, 0.20864 14.44 4.41076 0.91 
Mo Ka, 0.7078 4.31 15.2 3.14 
Cu Ka, 1.537 2.04 33.7 6.96 


Cr Kay 2.285 1.43 52.4 10.8 


The calculated rocking curves and theexperimental points for Cr Ka; taken 
on crystals III are shown in Fig. 4. It is seen that the calculated rocking curve 
agrees within the experimental error with that obtained experimentally. 
The half width at half maximum of the theoretical curve is 6.9”, the observed 
value is 7.1”. The percent reflection Po is 0.557 from Table VII; the observed 
value is 0.554. It is of interest to note that according to the analysis of the ac- 
tion of the double spectrometer presented in this paper, even if the diffraction 
pattern from each of the two crystals could be treated by Darwin’s equation, 
the observed percent reflection would only be 0.798 as against 1.00 which has 
been expected by some writers. 

The shape of the rocking curve to be expected from Prins’ theory has been 
calculated for Mo Ka,, Cu Ka; and Cr Ka. The rocking curve to be expected 
from W Ka, has been calculated from Darwin’s equation directly. The re- 
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sultant P's and half widths at half maximum, w, are compared with experi- 
ment in Fig. 5. The resolving power is calculated from a previous definition of 


the author,'? which is 
\/dd = Dd/2w. (31) 


This is to be interpreted as the resolving power which the instrument would 
have in the (1, 1) position with the same crystals. D is the dispersion of the 
instrument in the (1, 1) position. 
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Fig. 5. The half width at half maximum (w), resolving power (A, dA), and percent reflection 
(Po) of calcites 111 compared with theoretical values deduced from Prins. The solid curves are 
theoretical values, the points, experimental ones. 


Calculation of the coefficient of reflection 


In order to calculate the coefficient of reflection as defined in this paper 
we must find the area under the P curve (rocking curve) with the radian as 
the unit of angular measure. Thus 


i) [D(Ex.) + P(Ene sec 20) cos 20 |dEp, 


RK = Ad, aemsereresrantesieaen (32) 
(1 + cos 20) f I(&)d&, 





17S. K. Allison, Phys. Rev. 38, 203 (1931), Eq. (21) 208. 
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or, finally, 


f 8Gn dfn 
™ 1 + cos? 20 
R= — —_—— ——— 4§FZ-' cosec 26. (33) 


. 1 + cos 26 
f F (éo)dé, R 


The values of the integral in the numerator of Eq. (33) have been obtained by 
graphical integration of curves plotted from the data of Table VI, and are 
shown in columns 3 and 4 of Table IX. Columns 6, 7, and 8 of Table IX show 
the calculated and observed coefficients of reflection. 


TABLE IX. Calculation of R from Eq. (33) and comparison with observations. 


P(Epsldip, | 1+cos* 26 | RX 10° 


J : 
Darwin Prins | 1+cos 2@ | Darwin Prins Obs. 


Line nN 


W Ka 0.209A) 1.631 | 0.999 | 0.530 1.17 

Mo Ka; 0.708 “ 1.537 | 0.987 | 1.85 1.84 1.86 

Cu Ka; 1.54 “ 1.238 | 0.940 | 3.87 3.49 3.46 
0.917 5 





Cr Ka, 2.28 * 0.881 .65 4.53 4.68 


A graph of the information in Table IX, together with R values measured 
by other observers is shown in Fig. 6. Only the experimental data obtained in 
these experiments on crystals III are included in Fig. 6. 














12 14 
A(A) 

Fig. 6. Calculated and observed values of the coefficient of reflection (R). The dotted curve 
represents values calculated from Darwin; the solid curve, values from Prins. The open circles 
are experimental values of crystals III]. The solid circles are points by other observers. C— 
Compton;! WK—Wagner and Kuhlenkampff;! DS—Davis and Stempel;! DP—Davis and 
Purks.! 












Discussion of results 


The most obvious comment to be made on the agreement of observed and 
calculated coefficients of reflection from 0.5 to 2.3 A as illustrated for crystals 
III on Fig. 6 is that the agreement is very surprising since no temperature 
corrections have been applied. The coefficient of reflection should decrease 
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with rising temperature. The calculation of the temperature correction for a 
complicated crystal such as calcite is difficult and has not been attempted. 
The value of (sin @)/d for all the present measurements is 0.165 which is com- 
paratively small, and the temperature correction is proportional to this. There 
is, however, some evidence that the value FZ~'!=0.47 used here is too small. 
This value was based on the structure factors of Bragg and West." More re- 
cently James and Brindley'* and Pauling and Sherman" have published tables 
of structure factors. The exact value of FZ~! depends on the particular elec- 
tronic structure assumed for the molecule CaCQ; as it appears in the crystal. 
By assuming arrangements from Ca*?, C+4, O-? to Ca®, C°, O° we get values of 
FZ~' which in the table of James and Brindley range from 0.472 to 0.500, 
and from 0.496 to 0.518 in the table of Pauling and Sherman. If these newer 
structure factors are preferable to the older values, the value of 0.47 as used in 
this paper is somewhat too low.”° It is quite possible that this choice of a low 
structure factor is compensated by failure to make a temperature correction, 
thus giving the agreement of Fig. 6. This agreement does not extend to wave- 
lengths as short as W Ka, (0.2086A), and large deviations are apparent in the 
quantities plotted in Fig. 5 at this wave-length also. The laws of scattering of 
these relatively hard wave-lengths are undoubtedly different from the laws 
applicable in the range 0.5 to 2.3A, a large fraction of the scattering taking 
place by the Compton process. It is the opinion of the author that this lack 
of agreement shows that the theory of diffraction by a perfect crystal should 
be re-investigated in this shorter wave-length region. 

If the coefficients of reflection observed by Davis and Purks and entered 
in Fig. 6 are correct,*! they present a very serious difficulty, since they are 
far below the perfect crystal values. The Davis and Purks rocking curve 
widths for W Ka; (1, —1) were 2.5 seconds full widths at half maximum. The 
author has been unable to find any calcite specimens that even approach this 
small value, the values obtained ranging from 4.6 to 14.4 seconds in various 
specimens. The high values of the coefficient of reflection observed by other 
authors and entered in Fig. 6 are probably due to the use of calcite surfaces 
considerably less perfect than those used in the present experiments. 

The general statement that these results confirm Prins’ modification of 
Darwin's theory is incorrect in a very important aspect. The present results 
show that some sort of absorption correction to Darwin's results is needed, 
but really give no information at all as to the exact shape of the diffraction 
pattern from a single crystal. In particular, the present results give no infor- 
mation for or against the asymmetrical nature of the patterns in Prins’ theory. 
The fact that the calculated and observed values of P» are appreciably dif- 


18 James and Brindley, Phil. Mag. 12, 81 (1931). 

19 Pauling and Sherman, Zeits. f. Krist. 81, 1 (1932). 

20 T am indebted to W. H. Zachariasen of this laboratory for stimulating discussions of these 
matters. 

21 Recent communications from Professor Davis state that on repeating the experiments 
of Davis and Purks under discussion, the rocking curves from the same sample of calcite were 
found to be wider than at first reported. 
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ferent for wave-lengths other than Cr Ka; (at Mo Kq; calculated, 0.776; ob- 
served, 0.63) may be taken as evidence that the shape of the single crystal ' 
diffraction pattern is not exactly that of Prins. The value of Po is quite sensi- 
tive to the shape of the pattern, and independent of the polarization. For 
instance, let us compare three functions all having maxima equal to unity. 
If the shape of the diffraction pattern for a single crystal were given by exp 
(—&), we should find P)=0.71 from Eq. (26). In the same way if F(£) 
=(1+&)-!, P, =0.500, and if F(g) =(1+ &)~*, Po =0.625. 

A further question of interest is related to the possibility of finding other 
calcite surfaces as good as those of calcites III in this research. A. H. Comp- 
ton’ has reported rocking curves for Mo Ka; (1, —1) of 5 seconds full width at 
half maximum which agrees with the present results. Hoyt” has reported 
rocking curves under the same conditions somewhat less than 5 seconds wide. 
Allison and Williams,' found a full width at half maximum of 6 seconds. In 
view of these results it does not seem that such calcite surfaces are exceed- 
ingly rare. There is no evidence from the experiments reported here that the 
width of the rocking curve is a constant for a given piece of calcite, that is, a 
large sample from a known source. In the author’s experiments different 
splits from the same large piece showed different widths. 

Darwin, and other workers following him, have interpreted high values of 
R, for instance, those of Wagner and Kuhlenkampff in Fig. 6, as evidence for 
a mosaic structure of crystals. The present research shows that it is possible 
to obtain calcite surfaces for which there is no evidence of mosaic structure 
from measurements of the coefficient of reflection by double spectrometer 
methods in the wave-length region from 0.5 to 2.3A. Presumably the effects 
of a secondary structure as postulated by Zwicky* would not influence the 
present results. 








Note added in proof: Since completion of this article more accurate calcu- 
lations of theoretical curves have been carried out by eliminating the as- 
sumption of Eq. (17). This assumption states that the absorption correction 
can be rightly applied by assuming that the absorption of an atom in the 
unit cell is proportional to the fraction of the total number of electrons in the 
unit, all contained in the atom in question. This is rather far from correct in 
the absorption by the calcium atoms in calcite for wave-lengths shorter than 
the K-absorption limit of calcium. The more rigorous calculations direct from 
Eq. (15) show slight changes in the theoretical values, but do not invalidate 
in any way the general conclusions of this article. These new calculations will 
appear in a paper to be submitted from this laboratory in the near future. 


2 A. Hoyt, Paper delivered at the Washington Meeting of the American Physical Society, 
April, 1932. 
*3 Zwicky, Phys. Rev. 40, 63 (1932). 
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The Scattering of X-rays by Gases and Crystals 
By Y. H. Woo 


Department of Physics, National Tsing Hua University, Peiping, China 
(Received May 6, 1932) 
It is shown that, in exact agreement with the wave mechanical result worked 
out by Wentzel and Waller and Hartree, the Raman-Compton-Jauncey formula for 


the scattering of x-rays by an atom deduced on the basis of classical electrodynamics 
may be written 


Ie= LAP +2 — Let, (A) 


where /, is the intensity scattered in a direction @ with the primary beam, /, is the 
scattering from a single isolated electron as calculated by Thomson, Z is the atomic 
number, E, is the average amplitude of the waves scattered by the r-th electron in 
the atom and F is the true atomic structure factor as given by F= 2,°E,. This shows 
that the coherent scattering is proportional to F* and that the introduction of the 
average atomic structure factor F’ by Jauncey is unnecessary. On the basis of Eq. (A), 
theoretical formulas for the scattering of x-rays by gases and crystals are redeveloped 
and the results turn out to explain all the disagreement between theories proposed by 
Jauncey and the writer. 


N DISCUSSING the scattering of x-rays by polyatomic gases, Jauncey' 

has recently pointed out that the theoretical formula deduced by the 
writer’ is not exactly in agreement with that obtained by himself. Since the 
writer regards the atoms in a polyatomic system (a molecule or a crystal) 
as the scattering units, this disagreement should also occur when the scatter- 
ing of x-rays by crystals is dealt with. The purpose of the present note is to dis- 
cuss this matter in detail. 





Taking account of the modification introduced by Jauncey,’® the Raman- 
Compton! formula for the scattering of x-rays by an atom may be written 


I, =1.,F°+Z—-—F"/Z}, (1) 


where /, is the intensity scattered at an angle ¢ to a distance R, J, is the scat- 
tering from a single isolated electron as calculated by J. J. Thomson, Z is 
the atomic number and F’ is the average atomic structure factor. According 
to Jauncey,’ the relation between the average atomic structure factor F’ and 
the true atomic structure factor F is expressed by 


Z , 
FP? = F2 — (z LF? - r*) / (Z — 1), (2) 
r=1 


1G. E. M. Jauncey, Phys. Rev. 39, 561 (1932). 

2 Y. H. Woo, Phys. Rev. 39, 555 (1932). 

°>G. E. M. Jauncey, Phys. Rev. 37, 1193 (1931). 

*C. V. Raman, Indian J. Phys. 3, 357 (1928). A. H. Compton, Phys. Rev. 35, 925 (1930). 
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where /, is the average amplitude of the waves scattered by the r-th electron 
in the atom and where the true atomic structure factor F is given by 


Z 
F= DE,. (3) 
r=l 

In view of the similarity between Eq. (1) and the formula originally de- 
rived by Raman and Compton,‘ the part 7, F” was taken® by analogy to repre- 
sent the coherent scattering and the other terms to represent the incoherent 
scattering. However, there exists no physical or mathematical reason for 
justifying this separation. On substituting the value of F’ given by (2) into 

Eq. (1), we find as a result of straightforward calculation 


Z 

I, = I.) r+z- > pak (4) 
\ ’ l / 

This is exactly the same expression as that obtained by Wentzel® and inde- 

pendently by Waller and Hartree’ on the basis of wave mechanics. This shows 

that the coherent part of the total intensity is proportional to F? and that 

the introduction of the average atomic structure factor F’ by Jauncey is un- 

necessary. When corrected for the change of wave-length, Eq. (4) becomes 


Z 
Z— SE? 
ly 1 + 8,—-——_————_= (5) 
(1 + y(1 — cos ¢) | 
where y=/, mcX and hh, m, ¢ and X have their usual significance. 

Recently Herzog* has actually compared Eq. (5) with experiments re- 
ported by Wollan® on absolute measurements of the scattering of Mo K, rays 
by helium, neon and argon gases. The agreement seems to be quite satis- 
factory. 

On the basis of Eq. (5) and following the arguments presented by the 
writer in previous papers,'’ the theoretical formulas for the scattering of x- 
rays by gases and crystals can be easily redeveloped. For the scattering from 
a simple cubic crystal consisting of atoms of one kind, the diffusely scattered 
intensity is given by 


> Cf. G. E. M. Jauncey, Phys. Rev. 38, 1 (1931); 39, 561 (1932); Y. H. Woo, Phys. Rev. 
39, 555 (1932). 

° G. Wentzel, Zeits. f. Physik 43, 1 and 779 (1927). 

71. Waller, Zeits. f. Physik 51, 213 (1928); Waller and Hartree, Proc. Roy. Soc. A124, 
119 (1929). The more exact formula obtained by these authors differs from Eq. (4) by including 
some negative terms arising from consideration of Pauli's principle, but the contribution to the 
total scattering due to these negative terms is negligible in comparison with that due to the 
terms included in Eq. (4). 

* G. Herzog, Zeits. f. Physik 70, 583 and 590 (1931). Cf. also Y. H. Woo, Sci. Rep. Tsing, 
Univ. AI, No. 4 (1932). 

* E. O. Wollan, Phys, Rev. 37, 862 (1931). 

1° 'Y. H, Woo, Proc. Nat. Acad. Sci. 17, 467 (1931); Phys. Rev. 38, 6 (1931), and 39, 555 
(1932), 
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X-RAY DISPERSION 


Z 
Z- SE? 
I, = NIA — MP + ——____——___ 6 
UA ) fi+y(1 — cos ¢) | (6) 


For the intensity scattered by a gas molecule containing atoms J, - - - 7, 
- n arranged at fixed distances from each other, the formula becomes 


ES sin kS;; N 
Ton = 1, tt SE ene: «86 «ae 
~ ~ kS;; 2X (1 + y(1 — cos ¢) |? (4) 


In Eq. (6) N is the number of atoms per unit volume of the crystal and 
e~*™ is the temperature factor as calculated by Debye" and Waller” and in 
Eq. (7) F; is the true atomic structure factor for the i-th atom, k=(47/\) 
sin} ¢,s;; is the distance of any atom? from any other atom j.The other quanti- 
ties have the same meaning as those employed in Eq. (5). When applied to 
the scattering by a diatomic gas consisting of two like atoms, Eq. (7) gives 
for the scattering per electron 


1 Z 
-— DE? 
; T 52 sin kS\ F? = 
S=——=(1+- =) 5 +r ——~ (8) 
2Z1, kS JZ [1+ (1 — cos ¢)}3 


where s is the distance between the two atoms in the scattering molecule. 

It will be noticed that Eq. (6) agrees with the formula deduced by Jauncey 
and Harvey" provided the latter is properly corrected for the change of wave- 
length. Also a comparison of Eqs. (7) and (8) with those recently obtained by 
the writer“ and Jauncey"™ will clearly show that, if the coherent and inco- 
herent scattering be properly separated, the disagreement pointed out by 
Jauncey no longer exists. Finally, referring the comparison between theory 
and experiment recently made by the writer," it can be readily concluded 
that Eq. (8) is supported by Wollan’s measurements" on the scattering of 
Mo K, radiation by diatomic gases. 


1 P, Debye, Ann. d. Physik 43, 49 (1914). 

2 |, Waller, Zeits. f. Physik 17, 389 (1923); Upsala Dissertation, (1925). 
8 Jauncey and Harvey, Phys. Rev. 37, 1203 (1931). 

4 Y.H. Woo, Phys. Rev. 39, 555 (1932). 

16 G. E. M. Jauncey, Phys. Rev. 39, 561 (1932). 

1 E,O. Wollan, Phys. Rev. 37, 862 (1931); Proc. Nat. Acad. Sci. 17, 475 (1931). 
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Ionization of Air by y-Rays as a Function of Pressure 
and Collecting Field 


By I. S. BowEN 
Norman Bridge Laboratory of Physics, California Institute of Technology 
(Received May 23, 1932) 

The ionization of air by y-rays was studied for pressures from 1 to 93 atmos- 
pheres and for collecting fields from 1.55 to 1009 volts per cm. Increases in ioniza- 
tion current of over 40 percent were observed when the potential gradient was varied 
over this range, thus indicating that the lack of proportionality of ionization current 
with pressure obtained by previous observers was principally due to lack of saturation. 


I. INTRODUCTION 


HE ionization of air by y-rays or by cosmic rays has been studied through 
a wide range of air pressures by Downey,! Fruth,? Broxton,’ Swann,‘ and 
Millikan.’ All of these observers agree in finding that as high pressures are 


reached the ionization per atmosphere decreases to a small fraction of its 
value at one atmosphere. In fact, the total ionization does not increase ap- 
preciably as the pressure is pushed above 100 atmospheres. 

An explanation of this decrease has been proposed by Downey and elabo- 
rated by Broxton as follows. Practically ali of the ionization in a gas is caused, 
not by the y-rays ‘or the cosmic rays themselves, but by the secondary 8-rays 
ejected by them. At atmospheric pressure in any vessel of ordinary size the 
B-rays traversing the chamber are formed largely in the walls of the vessel. 
As the pressure is increased, however, the ranges of these 8-particles decrease 
until they finally become less than the dimensions of the chamber. When this 
condition is reached the 8-rays from the walls have produced all of the ions 
that they are capable of, and therefore any further increase of pressure cannot 
increase the total ionization caused by them. This theory appears then to 
predict the ionization pressure curves obtained by all observers. 

However, it is a well established fact that the mass absorption coeffi- 
cients of y-rays and cosmic rays are not widely different for light substances 
such as air or water from their values for the metals out of which ionization 
chambers are usually made. Since the y-rays or cosmic rays are absorbed in 
the air and since the only way in which their energy can be absorbed is by 
transfering it to secondary §-rays, it is obvious that the air must have ap- 
proximately the same efficiency of production of 8-rays as the walls. This 
being the case, any 6-rays from the wall that are absorbed by the air are re- 
placed by a corresponding number of 8-rays ejected from the air itself. If 

1K. M. Downey, Phys. Rev. 16, 420 (1920); 20, 186 (1922). 

2H. F. Fruth, Phys. Rev. 22, 109 (1923). 

* J. W. Broxton, Phys. Rev. 27, 542 (1926); 28, 1071 (1926); 37, 1320 (1931). 

4W.F.G. Swann, J. Frank. Inst. 209, 151 (1930). 

®*R. A. Millikan, Phys. Rev. 39, 397 (1932). 
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this is true, the number of ions formed should be directly proportional to the 
pressure. 

In order to reconcile this conclusion with the failure of all observers to 
collect this large a number of ions at the higher pressures, one is at once led 
to the suggestion, made independently by Millikan and Bowen® and by 
Compton, Bennett and Sterns,’ that saturation was not attained at these high 
pressures, in spite of the fact that all experimenters tested for saturation in 
the usual way. 

In all of these experiments the ionization chamber was a cylinder or 
sphere with a small collecting rod or fibre in the center. With this arrange- 
ment most of the potential drop is near the collecting rod while the great 
bulk of the chamber is subject to a potential gradient that, expressed in volts 
per cm, numerically amounts to only a very few percent of the total collecting 
potential applied. For this reason the highest potentials that these observers 
found it feasible to use, gave a gradient only a little above the point at which 
the ionization current became even approximately constant as the gradient 
was varied. Thus Broxton, in testing for saturation, varied his potential 
gradient through only 20 percent while the other experimenters carried their 
tests to only two or three times the gradient where the current became ap- 
proximately constant. Obviously these tests would not bring to light varia- 
tions of the type found in the experiments described in this paper. 


Il. DESCRIPTION OF APPARATUS 


As no indication of a lack of saturation was found by previous observers 
when the collecting potential was varied through a small range, it was neces- 
sary to build an ionization chamber in which very much higher potential 
gradients could be applied if a satisfactory test of this lack of saturation was 
to be obtained. Obviously the older type of chamber was unsuited to this be- 
cause of the very high total potentials necessary. Consequently a chamber 
was constructed in which the ions were collected between parallel plates, thus 
making possible a high uniform field. The ionization chamber was a steel 
cylinder 12.5 cm in diameter and 20 cm in length, designed to hold a pressure 
of 100 atmospheres. In this were placed 8 plates 10 cm in diameter which 
could be connected to a source of high potential. Alternating with these were 
7 collecting plates 7.5 cm in diameter. These collecting plates were surrounded 
by grounded guard rings from which they were insulated by small amber lugs 
in grounded metal shields. The inside diameter of the guard rings were 7.7 cm 
and the outside 10 cm. The collecting plates were separated from the poten- 
tial plates by 1 cm. This arrangement made possible the collection of ions 
from a definite known volume in which there was a uniform potential gfadi- 
ent which was numerically equal to the potential applied to the potential 
plates. 

Ideally it should have been possible to measure directly the ionization by 


® R. A. Millikan and I. S. Bowen, Nature 128, 582 (Oct. 3, 1931), See also reference 5. 


7A. H. Compton, R. D. Bennett and J. C. Sterns, Phys. Rev. 38, 1565 (Oct. 15, 1931) 
and Phys. Rev. 39, 873 (1932). 
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connecting an electrometer to the collecting plates and observing their rate 
of change of potential with time. Unfortunately however a change in poten- 
tial of the high potential plates due to the unavoidable variation in the bat- 
tery used to maintain it, produces by induction a change in potential of the 
collecting plates. Thus because of the large electrical capacity of the system, 
the ionization at the lowest pressure used, caused a change of only 0.3 volt 
per hour in the potential of the collecting plate. It was obviously impossible 
to maintain a 1000 volt battery so constant that its variation would not be an 
appreciable fraction of this. To compensate for this variation a second set of 
plates was built. This was an exact duplicate of the first set except that it was 
enclosed in a light brass case. The high potential plates of the ionization 
chamber A and the compensating chamber B were then connected to the same 
battery as shown in Fig. 1, while the electrometer / was arranged to measure 
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Fig. 1. Diagram of connections. 


the difference in potential of the two sets of collecting plates. The electro- 
meter used was of the single string type. The two plates of the electrometer 
were connected respectively to the two sets of collecting plates, while a po- 
tential of 225 volts was applied to the string. The electrometer was adjusted 
for a sensitivity of about 500 divisions per volt. 

Because of the fact that cosmic rays at sea level cause less than } of the 
ionization in a closed vessel, the remainder being due to y-rays from radio- 
active materials in the neighborhood, and to a-rays from radioactive impuri- 
ties in the walls of the ionization chamber, it was clear that more interpret- 
able results could be obtained if the ionization of a single source such as 
-rays were studied. 

For this purpose y-rays from a small quantity of radio-thorium, filtered 
through 13 mm of lead were used. In order to obtain a well defined beam of 
y-rays the radio-thorium and filter were placed in the center of a lead cylinder 
30 cm in diameter and 30 cm long. There was a conical opening leading to 
one end of the cylinder which allowed a beam to pass that was somewhat 
larger than the ionization chamber at its position about 2 meters away. The 
block of lead was on a swivel so that the beam of y-rays could be directed 
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either at the ionization chamber or turned to one side of it. By taking read- 
ings with the y-rays passing through the chamber and then again with it 
turned to one side it was possible to determine the ionization of the y-rays 
independent of the ionization of cosmic rays or a-rays. 

The usual procedure in taking readings was to connect the collecting 
plates of the compensating chamber B to ground while the collecting plates 
of the ionization chamber A were held at some potential, usually 0.05 volts 
below ground, by means of the potentiometer P. The electrometer was then 
read. Next the connections a and } were broken, thus leaving both sets of 
collecting plates floating. The electrometer was then immediately read again, 
and the stop watch simultaneously started. After approximately enough time 
had elapsed for the collecting plates of the ionization chamber to receive 
enough ions to bring their potential to the same amount above that of the 
collecting plates of the compensating cylinder as they had started below, the 
electrometer position was again noted and the stop watch read. Switches a 
and 6} were now closed and this same potential above ground was then applied 
to the plates by means of the reversing switch S and the electrometer again 
read. In this way it was possible to correct for any zero drift or change in 
sensitivity of the electrometer. 

This gave then with a considerable accuracy the rate of change of poten- 
tial difference between the collecting plates. The relationship between this 
rate of change of potential and the ionization per cc per sec. was obtained as 
follows: Let C, be the electrical capacity between the collecting plates and 
the high potential plates of the ionization chamber, C. the capacity to 
ground of the collecting plates, leads and electrometer plates. By symmetry 
these capacities are the same for compensating chamber. Let C; be the 
capacity between the plates of the electrometer. It is then easy to show that 
if a charge Q is accumulated on the collecting plates of the ionization cham- 
ber and Q2 on these plates of the compensating chamber the potential differ- 
ence in volts across the electrometer is W where 


Wo = (300(Q; — Q2))/(C1 + C2 + 2C3). (1) 


Likewise it may be shown that if the potential plates of the compensating 
cylinder are grounded while the potential of these plates in the ionization 
chamber are varied by an amount 2, then the charges induced on the collecting 
plates (left floating) are such as to change the potential between the collecting 
plates of the two cylinders by an amount w, where 


w= 01/(Ci + C2 + 23). (2) 


w is obviously the change in reading of the electrometer when this varia- 
tion takes place and hence can be determined directly. 
From (1) it is seen that the difference of the ionization currents in the two 


chambers is 
‘ dQ, dQe dW Cy +Ce2t+ 2C; 





—— 3 
dt dt dt 300 9) 


where dIV/dt is the rate of change of potential across the electrometer, in 
volts per sec. But substituting from (2) 
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dQ, dQ» dw vw; 


-+ (+4) 
dt dt dt 300w 


Since, however, the distance between the potential and collecting plates N 
is lcm 
C; == A tr 
where A is the total effective area of the collecting plates in each cylinder 
d02 dW oo A 


dt dt dt 300w 4x 


dQ; 


The difference in the number of ions (.V,;— V2) collected per ce per sec. 
in the two chambers is then, since the volume from which ions are obtained 
is A ce and the charge on each ion Is e, 


1 /dQ; = dQ dy 2 
NV; — No = — ) = - er 
el dt dt dt w300 XK 4re 


By taking the difference of the readings with and without the y-rays, the 
N» cancels out and the number produced by the y-rays alone are obtained. 

Since in these experiments y-rays giving an ionization current about 20 
times as large as the natural ionization current caused by cosmic rays, a-rays 
etc. were used, it was advisable to test whether the form of the ionization- 
voltage curves was dependent on the intensity of the radiation. For this pur- 
pose determinations of the ionization currents were made, at the two highest 
pressures, when the y-rays were cut down to less than { of their normal in- 
tensity by being passed through 5 cm of steel. 





The pressures up to 25 atmospheres were measured on a pressure gauge 
that was calibrated, during the course of the experiments, against a gauge 
tester. Consequently these values should be in error by less than 0.1 atmos- 
phere. It was not feasible to calibrate the gauge used on the highest pressure 
and therefore there may be a rather large error in this pressure. 

In all cases a given collecting potential was maintained on the collecting 
plates for several hours before readings were started to allow any charges ac- 
cumulated on the insulators to take up a steady value. 


III. ResUuLTs 


The ionization currents observed in these experiments for a series of pres- 
sures and potential gradients are given in Table I. The ionization current is 
expressed as the number of ions per cc per sec. per atmosphere. For the two 
highest pressures two columns of values are given. The first column marked 
I, represents the ionization at the full intensity of the y-rays, that is, the 
same intensity as that used for the lower pressures. The second column 
marked 5.187 J» is 5.187 times the number of ions obtained when the intensity 
of the y-rays was reduced by passage through 5 cm of iron. Except for the 
two lowest gradients at the 93 atmospheres pressure, the two columns agree 
within the experimental error, thus indicating that with these two exceptions 
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10.; 


79. 
86. 
89. 
94. 
102. 
108. 


TABLE I. y-rays. Number of ions collected per cc per second per atmosphere. 











the form of the curves is independent of the intensity for 
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Fig. 2 displays graphically these data, the number of ions per cc per sec. 
per atmosphere being plotted against the log of the potential gradient of the 
collecting field P. It is at once seen that, in agreement with previous obser- 
vers, the ionization per atmosphere falls off rapidly at the high pressures. 
Furthermore it agrees with their observation that when the potential gradi- 
ent is of the order of 5 or 10 volts per cm, doubling this gradient produces a 
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Fig. 2. The number of ions per cc per sec. per atmosphere (V) collected as a function of the 
potential gradient in volts per cm P for the pressures in atmospheres indicated. 


However when high potential gradients are applied a very material in- 
crease of the currents is obtained thus showing that the lower gradients do not 
collect all of the ions formed. Thus if we assume that the number of ions per 
atmosphere actually formed is constant and has the value 121.1 ions per cc 
per sec., then the number that one fails to collect at a gradient of 6.2 volts 
per cm (approximately the field used by previous observers) is (121.1 — N¢.2). 
The fraction of these that one still fails to collect at some other higher poten- 
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tial P is R=(121.1—N>), (121.1—N¢.2). This fraction then may be taken as 
a measure of the success attained in pushing towards saturation by using 
higher fields. This fraction R as calculated from Table I is plotted against 
the log of P in Fig. 3. 

From the figure it is seen that at 3.74 atmospheres it is possible by increas- 
ing the gradient from 6.2 to 1009 volts per cm to collect all but 21 percent 
of the ions that are not obtained with the lower gradient. At the higher pres- 
sures the increase of the field strength is progressively less effective. 

These results provide conclusive evidence that the falling off of ionization 
per atmosphere at higher pressures is largely due to lack of saturation rather 
than to the mechanism suggested by Broxton.°* 
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Fig. 3. R=(121.1— Np)/(121.1—N¢.2) as a function of the 
potential gradient P for the pressures indicated. 


Furthermore the lack of dependence of the form of the curves on the in- 
tensity of ionization, in nearly all cases, indicates that the lack of saturation 
is caused by recombination of the ejected electron with the parent ion rather 
than by random recombination between any ions formed. Presumably at the 
high pressures the ejected electron loses its energy before it succeeds in escap- 
ing from the region where the field of the parent ion is still very strong. This 
fortunately also means that this lack of saturation does not invalidate the 
results of many cosmic ray observers who have used high pressures, since the 
lack of saturation merely cuts down all currents by a constant factor and 
in no way effects the relative values on which absorption coefficients are 
based. 

This increase of ionization current with potential gradient also seems to 
call for a modification of the theory developed by Compton, Bennett and 
Sterns’ to explain this lack of saturation. According to their theory there is a 

8’ There may still be a slight variation of the type suggested by Broxton due to a small 
difference in efficiency of production of 8-rays by the air and by the steel walls of the chamber. 
Unpublished studies on the effect of wall materials by Workman in this laboratory indicate 


however that this effect of these materials is at most 30 percent and is probably in the opposite 
direction, i.e., it should cause the ionization per atmosphere to increase slightly with pressure. 





















31 





IONIZATION OF AIR BY y-RAYS 


certain critical sphere of definite radius surrounding each ion. If an oppositely 
charged ion finds itself inside of this sphere the force of attraction over 
balances the tendency to diffuse away and the two ions always recombine. 
On the other hand the tendency to diffuse away predominates for an ion out- 
side of this sphere and the ion always escapes. Since at the surface of this 
critical sphere the field due to the ion at the center is 40,000 volts per cm any 
ordinary collecting field should have no effect. The question whether recom- 
bination with the parent ion takes place or not is determined, according to 
this theory, solely by the distance from the parent ion that the electron is 
ejected by the 6-particle, and the factors that enter into the tendency to dif- 
fuse away such as the temperature. 

Obviously this picture is much too simplified. Thus given two ions in a 
gas there is always a certain probability that one will diffuse close enough to 
the other, regardless of their original distance apart, so that recombination 
takes place. Of course this probability falls off rapidly as the distance in- 
creases. Furthermore this probability, particularly for an ion at a consider- 
able distance from the parent ion, is materially changed as the strength of an 
external collecting field is varied. This of course changes the number that 
succeed in escaping from the parent ion and therefore one would expect va- 
riation in the strength of the ionization current similar to that shown in the 
Table I and Fig. 2. 

Qualitatively this is also in agreement with the behavior indicated in 
Fig. 3. At the lower pressures all of the ions start at a great distance from the 
parent ion where its field is small and consequently the probability of recom- 
bination should be largely effected by the size of the external field, as it is 
found to be. On the other hand for higher pressures most of the ions start so 
close to the parent ion where its field is large compared to that of any feasible 
external field. This being the case very little effect of the strength of the ex- 
ternal field should be expected. 


TABLE II. Residual ionization. Number of ions collected per cc per second per atmosphere. 


Field Pressure 
in volts per cm inatm. 24.95 , 93. 








5 
3.55 3.11 1.34 
6.2 3.22 1.45 
23.0 3.38 1.55 
91.5 3.53 1.64 
367. 4.00 1.77 

1009. 4.42 1.99 





The readings made when the y-rays were turned to one side represent 
the difference in the residual ionization caused in the two chambers by cosmic 
rays and radioactivity of the chamber and room. This was large enough in 
the case of the upper two pressures to be read with a considerable accuracy. 
Also in the case of these high pressures, the ionization in the ionization cham- 
ber was large compared to that in the compensating chamber and consequently 
the fact that the latter is subtracted should not modify the values very 
greatly. These results are given in Table II. As is at once seen the effect of 
increase in collecting potential is essentially the same as with the y-rays. 
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A Search for Isotopes of Hydrogen and Helium 


By WALKER BLEAKNEY* 
Palmer Physical Laboratory, Princeton University 


(Received June 1, 1932) 


The well-known triatomic hydrogen ion (H'H'H!)* of mass 3 is formed by elec- 
tron impact at a rate proportional to the square of the pressure if the pressure is very 
low. The isotopic molecular ion (H'H?)*, also of mass 3, would be expected to vary lin- 
early with pressure. These two facts constitute a method for distinguishing between 
the two with a mass spectrograph whose resolving power is insufficient to separate 
them by magnetic analysis. The method becomes successful only when the working 
pressures are very low. With a mass spectrograph designed for low pressure work evi- 
dence was found which amply confirms the existence of an isotope of hydrogen of 
mass 2 and also gives a fair estimate of the abundance. In a sample of ordinary com- 
mercial electrolytic hydrogen the data indicate an abundance ratio 


H?/H! = 1/30000 + 20%. 
The same ratio for a sample of hydrogen which had been concentrated by Urey, Brick- 
wedde and Murphy came out to be 

H?/H! = 1/1050 + 5%. 


An effort was made to detect a difference in mass between the ions (H'H'H!)* and 
(H'H?*)* without success. Considering the resolving power of the apparatus this failure 
was interpreted as meaning that the packing fraction of H? is greater than 4X 10~i.e., 
the atomic weight is greater than 2.008. A search was also made for isotopes 3 and 5 of 
helium but none were found. It is concluded that their abundance must be less than 
one part in 50,000. 


INTRODUCTION 


HE hydrogen and helium atoms play such fundamental roles in all 

physics that any new properties or characteristics of their nuclei are of 
wide interest. There have been indications that two or more isotopes of these 
elements may exist and several attempts have been made from time to time 
to find them experimentally. The first positive evidence for the existence of 
an isotope of hydrogen of mass 2 has been obtained recently by Urey, Brick- 
wedde and Murphy' whose paper also gives a review of previous work and a 
bibliography of the subject which need not be repeated here. Their photo- 
graphs of the Balmer spectrum under very high dispersion reveal weak satel- 
lite lines attributable to hydrogen atoms of mass 2. These lines were en- 
hanced in a sample of hydrogen in which enrichment of the heavy isotope 
had been attempted by evaporation near the triple point. Further evidence 
for the existence of an isotope in hydrogen may be found in Allison's inter- 
pretation of his experiments on optical rotation? but his method is as yet little 


* National Research Fellow. 

'H. C. Urey, F. G. Brickwedde and G. M. Murphy, Phys. Rev. 40, 1 (1932); also Phys. 
Rev. 39, 164 (1932) and 39, 864 (1932). 

* F. Allison, Jour. Ind. Eng. Chem. 4, 9 (1932). 
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understood. The importance of the problem seemed to warrant an inde- 
pendent investigation from a different angle and accordingly a careful re- 
examination of the hydrogen and helium ions was undertaken with a mass 
spectrograph. 

There are many complications which arise in the magnetic analysis of 
hydrogen ions because of the various primary and secondary combinations 
involved. Let us assume that a rare isotope of mass 2 does exist and let us 
enumerate the possibilities when the ions are formed at a very low pressure by 
electron impact. It is well known* that under these conditions monatomic, 
diatomic and triatomic ions are produced. Hence we should expect the vari- 
ous combinations listed in Table I. The column labeled f(p) gives the in- 
tensity or number of ions as a function of pressure when p is so small that 


TABLE I. Possible ions resulting from a mixture of two hydrogen isotopes 
in which H' is much mere abundant than H?. 




















lon m/e f(p) Intensity 
- (H!)* 1 apthp weak 
2 (H'H!')* 2 ap v. strong 
3. (H?)* 2 asp+hsp* v. weak 
4. (H'H'!H!)* 3 bp? weak 
5. (H'H?)* 3 asp weak 
6. (H?H?)* 4 asp v.v. weak 
7. (H'H'H?)* 4 bsp? v.v. weak 
8. (H'H?H?)* 3 bep? v.v. weak 
9, (H?H?H?)* 6 bp v.v. weak 











secondary reactions are small compared with primary ones. In general the 
primary processes will vary linearly with pressure while secondary ones will 
be proportional to the square. Of the configurations listed in Table I it is to 
be expected that 6, 7, 8 and 9 will be too weak for observation. Number 1 is 
of no interest in this particular case and 3 will be completely masked by 2. 
Since 2 is a linear function of pressure and since it is very strong it serves as 
a good index of the pressure. There are left 4 and 5, two ions of comparable 
abundance but different functions of pressure and it is in this direction that 
the following studies were made. In the case of helium no such complications 
arise. 
EXPERIMENTAL PROCEDURE 


The apparatus used in this experiment has been described‘ in a previous 
number of this journal and some preliminary results have also been re- 
ported.®:* The flow method was used in every case, the gases streaming into 
the tube through fine capillaries or a palladium tube and out through diffu- 
sion pumps. The lighter the atom the faster will it leak through the capillary 
but presumably the same selective effect will occur in the diffusion pumps so 
that the relative concentrations of isotopes in the ionization chamber will 


3H. D. Smyth, Rev. Mod. Phys. 3, 347 (1931). 

4 W. Bleakney, Phys. Rev. 40, 496 (1932). 

5 W. Bleakney, Phys. Rev. 39, 536 (1932) (Letter). 

* W. Bleakney, Phys. Rev. 40. 130 (1932) (Abstract). 
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remain unchanged after equilibrium has once been established. No difference 
in the results on hydrogen were observed when the palladium tube was sub- 
stituted for the fine capillary. 

Through the kindness of Professor Urey and Dr. Murphy a sample of the 
residue of hydrogen evaporated near the triple point by Dr. Brickwedde at 
the Bureau of Standards was provided. The details of the preparation of this 
hydrogen have been described in their paper.' This particular sample was 
designated by them “Sample III.” Ordinary commercial electrolytic hy- 
drogen was used for the rest of the work. 

The procedure in the case of hydrogen was to measure the total number of 
ions of mass 3 reaching the collector, per unit electron current, as a function 
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Fig. 1. A typical e/m analysis curve for ordinary hydrogen. p 10> mm Hg. H = 400 
gauss. 1’, =60 volts electron velocity. 


of the pressure. The pressure itself could be measured in two ways. The first 
method consisted in observing the intensity of the ion peak of mass 2 which 
is proportional to p. In the second method the tube was used as an ionization 
gauge by observing the total positive ion current per unit electron current. 
The two methods gave substantially the same results. At the low pressures 
used the McLeod gauge could not be read with any accuracy. After baking 
the tube the residual pressure of foreign gases was always in the neighbor- 
hood of 10~* mm. The working pressures ranged from 10~ to 10-* mm. 


RESULTS FOR HYDROGEN 


Some idea of the resolving power of the apparatus may be had from Fig. 1 
which shows the three typical peaks obtained in ordinary hydrogen with a 
magnetic field of 400 gauss and a pressure of about 10-> mm. The electron 
velocity was 60 volts. At higher magnetic fields the peaks become somewhat 
sharper. Fig. 2 represents the number J of ions of mass 3 as a function of the 
number of ions of mass 2. Both ordinate and abscissa are plotted in the same 
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arbitrary units. Curve I was observed for ordinary hydrogen and curve III 
for Brickwedde’s concentrated sample. The difference between these two 
curves is given by the straight line IT. 


— = 


12x10" 7 ZX | 
ne ae 


Le azal 


0 20 40 60 &0 
pw (HH!)* 
Fig. 2. Total number of ions (H'H!H!)*++(H!H?)*+ of mass 3 as a function of the number 


(H'H!*)* of mass 2. Curve I, ordinary commercial electrolytic hydrogen. Curve ITI Brickwedde’s 
concentrated sample. Curve IJ, difference. 
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These data are interpreted by the author in the following manner. As in- 
dicated in Table I we should expect the number of ions of mass 3 to be given 
by 


I = ap + bp’. (1) 


Now it is significant that curves I and III may be fitted by Eq. (1) using a 
different value for the constant a but the same value for } in each case. The 
circles represent the experimental points while the curves are plotted from 
the following equations. 


Curve I Ir; =( 6.6p— p*) X 10° (2) 
Curve III Tin => (190.0p —_ p’) x 10 $ (3) 
Curve Il Jy, = Jar — J; = 183.4 X 10 5p. (4) 


The left-hand member of Eq. (4) is interpreted as the difference between 
the numbers of isotopic ions (H'H?*)* in the two samples. Since / is actually 
measured in terms of the number of (H'H')* ions it is evident that the slope 
of curve II gives the increase due to the concentration process. 


H ' H2/H' H! = 1/546. (5) 


Here we assume of course that the probabilities of ionization of the two types 
are the same. The total concentration of the H'H? molecules is given by the 
coefficient a of Eq. (1). For curve III this is 


H! H2/H' H! = 1.90 X 10-* = 1/526 (6) 
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and for ordinary hydrogen 
H'! H?/H' H' = 6.6 X 10-5 = 1/15100. (7) 


The probable error in these two results is estimated to be 5 percent and 
20 percent respectively. No attempt was made to treat the data by the 
method of least squares in detail since neither the ordinate nor the abscissa is 
known precisely and the probable error of each observation varies from point 
to point. The curves were therefore fitted in a rather arbitrary way. If we 
write the ratio of numbers of atoms instead of molecules we have 


H' H? = 1050 + 5% (8) 
for the concentrated sample and 
H!, H? = 30,000 + 20% (9) 


for ordinary hydrogen. It will be noticed that there is quite a marked dis- 
crepancy between these results and those first reported by Urey, Brickwedde 
and Murphy,' but they have largely succeeded in accounting for the differ- 
ence’ by considering the relative absorption of strong and weak lines in their 
discharge tube. 

In order to show the results in a more striking manner it is illuminating 
to divide Eq. (1) by the pressure and so obtain 


I/p=a-+t bp (10) 


an equation of a straight line whose intercept a gives the abundance of the 
isotope sought. Curves I and III of Fig. 3 plotted in this manner represent 
the same data as are shown in Fig. 2. This method of presenting the data im- 
proves one’s confidence in the evidence for the existence of the isotope in 
ordinary hydrogen but on the same scale the points for the concentrated 
sample appear scattered over a considerable range. However a little con- 
sideration will probably convince the reader that the probable error chosen 
is not too small. The curves I and III are plotted from Eqs. (2) and (3). 

It is of importance to point out that the triatomic ion (H'H'H!)* is not 
only a function of the pressure but also of the geometry of the apparatus and 
the fields applied. The coefficient b of Eq. (1) is related to the intensity of this 
ion and is given by the slope of the curves in Fig. 3. The relative number of 
isotopic molecules, however, should be independent of these conditions. 
Curve IV was taken with an altered set of electric and magnetic fields and it 
is evident that the slope } has been greatly increased while the intercept re- 
mains unchanged within the experimental error. The fact that curves I and 
IV have exactly the same intercept is not a coincidence. Since their probable 
errors overlap they were made to coincide by the author in fitting the curves. 

Some data for another sample of hydrogen are shown by curve V in Fig. 
3. This sample was prepared by Professor Urey by a diffusion method and 
represents a concentration of the lighter isotope. The data are admittedly 
rather inaccurate but the results seem to indicate less of the heavier isotope 


7H. C. Urey, F. G. Brickwedde and G. M. Murphy, Phys. Rev. 40, 464 (1932). 
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than is to be found in ordinary hydrogen. Professor Urey estimated the con- 
centration of H'H? in this sample to be about one tenth of that in ordinary 
hydrogen, a result which is not inconsistent with the somewhat higher one 
obtained here when account is taken of the rather large probable error in- 
volved in locating the intercept of curve V. 

Kallmann and Lasareff* have attempted to verify the existence of the H? 
isotope by seeking for the (H'H!H2)* ion of mass 4, but were unable to find it. 
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Fig. 3. Jp plotted asa function of p. Curve III, Brickwedde’s concentrated sample. Curves 
I and IV, ordinary commercial electrolytic hydrogen each taken under a different set of fields. 
Curve V, sample of hydrogen enriched in H' by diffusion process. 


Perhaps the probability for its formation is very small. They then tried the 
method of this experiment but at the lowest pressures possible in their ap- 
paratus (5X10) the ratio of ions of mass 3 to those of mass 2 was 1:4000 for 
ordinary and concentrated samples alike. The conclusion they reach is that 
their concentrated sample is at most one tenth as good as the one used in 
this work. If this is true then there is probably no contradiction between the 
two experiments. 

An accurate determination of the mass of the new isotope is of the great- 
est importance. Unfortunately the mass spectrograph employed has not a 
high dispersion on the mass scale. Nevertheless a trial was made to determine 
if possible a lower limit for the packing fraction. First ordinary hydrogen was 
admitted and the shape of the peak corresponding to ions of mass 3 was ex- 


8 H. Kallmann and W. Lasareff, Naturwiss. 20, 206 (1932). 
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amined very carefully. The pressure was such that this peak was made up 
almost entirely of (H'H'H')*. Then a mixture of the two samples was ad- 
mitted at such a pressure that (H'H'H!)* and (H'H?)*+ were present in equal 
numbers and again the shape of the peak was examined. If the two ions are of 
different weight one would expect the peak in the latter case to appear some- 
what broadened. No such broadening could be detected. It is estimated 
from these measurements that the packing fraction of H* lies between 
4&10°* and 12X10", 


or 2.024 > IW > 2.008 


where IV represents the atomic weight.* 

It may be well to point out that all the results of this paper simply point 
to the conclusion that a small number of normal molecules of mass 3 does 
exist in hydrogen. Conrad® has found evidence for the formation of neutral 
H; molecules in a discharge tube but it has been shown'® on quantum me- 
chanical grounds that three normal H atoms cannot combine to form a stable 
triatomic molecule. Hence if a neutral triatomic molecule exists it must be in 
an excited state in which case the ion would be formed as the result of a 
secondary or tertiary process and would not be a linear function of the pres- 
sure. Even without considering the weight of the spectroscopic evidence'! the 
most reasonable assumption which accounts for the results of this experi- 
ment is the presence of H'H? molecules in hydrogen. 


RESULTs FOR HELIUM 


Helium purified in a misch-metal arc was admitted to the tube and a 
search made for isotopes 3 and 5. None were found. It is believed that a con- 
centration of one in 50,000 or less could have been detected. This negative 
result agrees with that of Kallmann and Lasareff* but is a test not as dis- 
criminating as theirs. 

The author is indebted to Professor Urey and Drs. Brickwedde and 
Murphy for the samples of hydrogen which made this experiment possible. 
He is also grateful to the members of the Palmer Laboratory staff for their 
kindly interest and timely suggestions. 


* The author has just recently learned through private communication that Dr. K. T. 
Bainbridge has succeeded in measuring the atomic weight of H*? with his mass spectrograph. 
His result is W= 2.0135. 

® R. Conrad, Zeits. f. Physik 75, 504 (1932). 

10 WW. Heitler and F. London, Zeits. f. Physik 44, 455 (1927). 

11 The existence of the H,? line has been confirmed by Shenstone and Turner using a grat- 
ing in this laboratory and some of the “sample II1” hydrogen. 
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Elastic and Inelastic Electron Scattering in Hydrogen 
By A. L. HuGues* and J. H. McMILLEN 


Washington University, St. Louis 


(Received May 16, 1932) 


A study of both the energy distribution and the angular distribution of electrons 
scattered by hydrogen molecules was made. Elastic scattering was investigated in the 
range between 40° and 165°, for collision energies between 35 and 200 volts. The 
number of electrons scattered at any angle was smaller the greater the energy of im- 
pact. In all the curves a minimum was found between 90° and 110°, while for the 35 
and 50 volt collisions the curves rose to a definite maximum at 155° to 160°. Inelastic 
scattering was found to occur with energy losses ranging from a minimum of about 12.0 
volts, through a sharply defined most probable loss at 12.6 volts, to a complete loss 
of all the energy. It was found that as the energy of collision was increased the proba- 
bility of the larger energy losses increased with respect to the probability of the 
smaller energy losses. Angular distribution curves were steeper, the smaller the amount 
of energy lost at collision for a given collision energy. Also, for any particular energy 
loss, the steepness increased with the speed of the colliding electrons. The electrons 
which are torn off from atoms in ionizing collisions are termed ejected electrons and are 
presumably those which appear in our experiment with but little or no energy. The an- 
gular distribution of the ejected electrons, having various amounts of energy (1 to 8 
volts), were studied for different collision energies (35 to 340 volts). In every case there 
was a noteworthy absence of considerable scattering at small angles, so characteristic of 
the scattering of faster electrons. In many cases, pronounced peaks were found in the 
angular distribution curves at large angles (90° to 160°). These peaks change in pos- 
ition and size in a regular manner with the energy of collision and the energy of the 
ejected electron. 


HIS paper is an account of a continuation of our experimental investiga- 

tions on electron scattering in gases. The results of our investigations on 
argon have already been published in this journal.! As argon is monatomic, 
it seemed desirable to select a diatomic gas for the second investigation; 
hence hydrogen was chosen for study. The method of experimentation with 
hydrogen was identical with that used in the work on argon. 

Scattering of electrons by hydrogen has been studied by Harnwell,* 
McMillen,’ Arnot,‘ Bullard and Massey,’ and by Ramsauer and Kollath.* 
Harnwell and McMillen found that the number of electrons scattered 
elastically fell off steeply as the angle was increased (angle range 7° to 60°; 
energy range 40 to 180 volts). Arnot, working with electrons between 30 and 


* This work was made possible by assistance to the senior author from a grant made by 
the Rockefeller Foundation to Washington University for research in science. 

1 A. L. Hughes and J. H. McMillen, Phys. Rev. 39, 585 (1932). It is convenient to men- 
tion here that several paragraphs were accidentally omitted from the end of the paper referred 
to. The omission was rectified in a letter to the Editor of the Physical Review 40, 469 (1932). 
An error in Eq. (2) of the same paper was also corrected. 

2G. P. Harnwell, Phys. Rev. 33, 559 (1929); 34, 661 (1929); 35, 285 (1930). 

3 J. H. McMillen, Phys. Rev. 36, 1034 (1930). 

*F. L. Arnot, Proc. Roy. Soc. 133, 615 (1931). 

5 E. C. Bullard and H. S. W. Massey, Proc. Roy. Soc. 133, 637 (1931). 

6 C, Ramsauer and R. Kollath, Ann. d. Physik 12, 529 (1932). 
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800 volts energy, confirmed these results, but found that, at about 90°, the 
curve representing the number of elastically scattered electrons began to in- 
crease. (His angular range extended to 120°.) Similar results were obtained 
by Bullard and Massey who studied the elastic scattering of electrons having 
energies between 4 and 30 volts. Ramsauer and Kollath also found, for very 
slow electron (1.5 to 11.0 volts), a definite but rather slow increase in the 
elastic scattering as the angle was increased above 90° to as far as 165°. So 
far as we know there have been no investigations of the distributions of 
electrons scattered inelastically by hydrogen. 


EXPERIMENTAL METHOD 


The same apparatus was used in the study of electron scattering in hy- 
drogen as in the previous work on argon.’ The reader is therefore referred to 
the previous article for a description of the apparatus and method of experi- 
mentation. Hydrogen was prepared by the action of pure dilute sulphuric 
acid on pure zinc. The gas was dried by passing through a tube immersed in 
liquid air. It was then stored in a reservoir over mercury and admitted to the 
apparatus through a capillary tube and a liquid air trap to freeze out the 
mercury. A suitable constant pressure could be maintained in the apparatus 
by regulating the driving pressure of the hydrogen on the high pressure side 
of the capillary tube, the gas being pumped out of the apparatus by a con- 
densation pump. This constant flow method provides us with pure hydrogen 
in the scattering experiments. In experiments on scattering in hydrogen it is 
particularly important to avoid impurities, as the scattering by a hydrogen 
molecule is small compared with the scattering by heavier molecules. 

Our results on hydrogen are expressed in the same units as those used in 
the previous article; that is, we measure the number of electrons diverted 
from the main beam at an angle @ and then multiply by sin @ to compensate 
for the change with angle of the effective scattering path length. The final 
values are therefore proportional to the scattering coefficient as defined on 
page 589 of our previous paper. * 


ELASTIC SCATTERING 


In view of the recent comprehensive work of Bullard and Massey and of 
Arnot on elastic scattering of electrons in hydrogen, we did not make an ex- 
tensive study of the elastic scattering. It seemed worth while, however, to 
attempt to extend the studies to large angles, inasmuch as the investigators 
referred to did not go beyond about 120°. Our results are shown in Fig. 1. It 
will be seen that we confirm the discovery of Bullard and Massey and of 
Arnot that there is a minimum at about 90°. In addition, we find a definite 
maximum at large angles which is contrary to the expectations of Arnot who 
thought that “the curves probably continue to rise steadily as the angle is 
increased to 180°.” 

We have indicated the positions of Arnot’s 29 volt curve and Bullard 


7 A. L. Hughes and J. H. McMillen, Phys. Rev. 39, 585 (1932). 
8 A. L. Hughes and J. H. McMillen, Phys. Rev. 39, 585 (1932). 
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and Massey’s 30 volts curve in the neighborhood of our 35 volt curve. It is 
evident that although all three curves agree in giving a minimum at about 
90°, no two curves are in satisfactory accord. We can only say, therefore, that 
there is qualitative agreement insofar as the curves rise and fall together. 


INELASTIC SCATTERING 


The distributions of energy lost by 50 volt primary electrons, when 
scattered at various angles, are shown in Fig. 2. The various curves have been 
arbitrarily adjusted to the same height for a loss of energy amounting to 20 
volts. It will be seen that the “most probable” energy loss is approximately 
12.6 volts. It is evident that the larger energy losses are relatively more 
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Fig. 1. Elastic scattering of electrons by hydrogen molecules. The numbers attached to 


the curves represent the electron energies in volts. Circles with dots, Arnot’s curve for 29 volt 
electrons. Triangles, Bullard and Massey's curve for 30 volt electrons. 





























probable for the larger scattering angles. A small subsidiary maximum ap- 
pears at about 16 volts in the 30° and 45° curves. In Fig. 3 we have the energy 
distributions of electrons, all scattered at 10°, for various primary energies 
between 35 and 200 volts. (The curve for 412 volt electrons, scattered at 5°, 
is included). As before, the curves are arbitrarily fitted together at the 20 
volt loss. It is seen that, for a given scattering angle, here 10°, the 12.6 volt 
loss becomes relatively more probable the /ower the primary electron energy. 
It will be noticed that there is a distinct slowing up in the rate of the fall of 
the curve between 15 and 20 volts, and in one curve (the 412 volt, 5° curve) 
there is a small subsidiary maximum at 17.5 volts. 

Various energy losses, viz., 12.6, 16.4, and 25.0 volts, were selected for a 
study of their angular distributions. The 12.6 volt loss, the most probable 
loss, is an excitation loss. The 16.4 loss is just above the ionization potential 
which is 15.9 volts. (The 15.9 volt loss was not studied, as the resolution of 
the apparatus was imperfect and consequently the electrons selected at a 
15.9 volt setting would really include a number of electrons causing excita- 
tion). Finally the 25 volt loss was chosen because it was well above the min- 
imum ionization loss and, at the same time, the yield was not too small to 
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Fig. 2. Energy distribution of electrons scattered inelastically by hydrogen molecules, 


showing variation with angle. Primary energy in all cases =50 volts. All curves are adjusted 
arbitrarily to the same height at the 20 volt loss. On the 30° curve, the positions of the theo- 
retically possible types of effects are indicated by suitable vertical lines. The associated hori- 
zontal lines indicate the possible vibrational levels. 
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Fig. 3. Energy distributions of electrons scattered inelastically by hydrogen molecules, 
showing variation with primary energy. All curves are adjusted arbitrarily to the same height 
at the 20 volt loss. 
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make measurements unreliable. According to Fig. 4, the less the energy loss, 
the steeper the angular distribution curve, a result already implied in the 
previous paragraph. The curves are all adjusted to the same value at 10°. 
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Fig. 4. Angular distribution of electrons scattered inelastically. Primary electron energy 
in all cases = 50 volts. Energy losses in volts are indicated on the curves. All curves are adjusted 
arbitrarily to the same height at 10°. 
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Fig. 5. Angular distribution of electrons scattered with 12.6 volts loss of energy, showing 
dependence on primary electron energy. The primary energy, expressed in volts, is indicated 
for each curve. 


Finally, in Fig. 5, we show how the angular distribution curves for a 
constant loss, viz., 12.6 volt (corresponding to excitation of the molecule), de- 
pend on the energy of the electrons before collision. Up to about 90° or so, 
the curves are steeper the higher the primary energy. (Too much reliance 
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should not be placed on the exact shape of the curves above 90°, as the read- 
ings were extremely small). 


EJECTED ELECTRONS 


Following the terminology used in our previous paper on argon,° we call 
those electrons which appear with very low energies the ejected electrons. 
When ionization takes place, an electron is detached from the atom. Con- 
sequently, after an ionizing collision we have two electrons; the faster we as- 
sume to be the colliding electron and the slower the ejected electron. (We 
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Fig. 6. Angular distributions of ejected electrons. Each row represents a single ejected 
electron energy which takes the values 1, 3, 5.5, and 8 volts as we go from one row to the next. 
Each column corresponds to a different energy of collision. All curves are to the same scale. 


cannot, of course, demonstrate that the slower of the two electrons is the one 
ejected from the atom; it seems probable, however, and suggests a con- 
venient terminology.) For a study of the angular distribution of ejected 
electrons we selected the 1, 3, 5.5, and 8 volt ejected electrons resulting from 
a collision between hydrogen molecules and primary electrons having energies 
ranging from 35 to 340 volts. The results are shown in Fig. 6. Perhaps the 
most striking feature common to all the curves is the absence of the exces- 
sively large scattering at small angles which is so characteristic of electrons 
scattered inelastically with considerable amounts of energy (see Fig. 4). The 
next feature worthy of mention is the presence, in many cases, of certain 
highly preferred directions for the emission of ejected electrons. A certain 
regularity in the way in which these peaks shift can be traced. Thus for any 


* A. L. Hughes and J. H. McMillen, Phys. Rev. 39, 585 (1932). 
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one ejected electron energy, an increase in the energy of the primary electrons 
is accompanied by a progressive shift of the position of the peak to smaller 
angles. On the other hand, if we keep the primary energy fixed and study the 
various groups of ejected electrons, we find that the higher the energy of the 
ejected electrons the /arger the angle at which the peak is to be found. For 8 
volt ejected electrons the angular distribution curves have become practi- 
cally featureless. It is to be noticed, too, that the peaks in the angular dis- 
tribution curves for all the ejected electrons practicany disappear when the 
primary energy is low (35 and 50 volts). 
DISCUSSION 

Several theoretical physicists have amplified and extended the pioneer 
work of Born, who first applied wave mechanics to the scattering of electrons 
by atoms, with the result that certain features of electron scattering are now 
accounted for theoretically. Very little, however, has been done on the corre- 
sponding problem for molecules. Massey and Mohr"? have carried through a 
calculation of the scattering of electrons by diatomic molecules. It is assumed 
that the diatomic molecule may be represented by a suitable potential field 
having axial symmetry over which passes an electron wave. The amplitude of 
the wave scattered in any desired direction is calculated. An average is then 
taken for all possible orientations of the molecular axis with respect to the 
direction of the primary wave. Finally it is found that the scattering by such 
a diatomic molecule is obtained on multiplying what is substantially the 
scattering by the constituent atoms by a diffraction factor (1+sin x/x), in 
which x =47d sin (@/2)/A, where d is the distance between the atoms in the 
molecule and X the electron wave-length. This diffraction factor changes sign 
periodically as x increases. Because of the extremely steep character of the 
curve for atomic scattering, however, it is not easy to detect the superposed 
oscillation of the diffraction pattern.'' Thus Massey and Mohr were unable to 
show that Arnot’s experimental curves (for 80 volts and above) for the scat- 
tering of electrons by hydrogen were in better accord with their formula for 
scattering by molecular hydrogen than with the formula for scattering by 
atomic hydrogen. Massey and Mohr state that the scattering curves obtained 
by McMillen” fall off much more rapidly with angle of scattering than the 
calculated values. This is incorrect, for McMillen’s curves are almost super- 
posable on Arnot’s, and if anything, they are less steep. Arnot’s range of 
angles was nearly twice as great as McMillen’s, and so afforded a better test. 

The rise in the scattering curves with angle as the angle is increased be- 
yond about 90° as shown in Fig. 1 of this paper and in Arnot’s work is not 
accounted for by Massey and Mohr’s theory. This theory is probably too 
simplified to do more than describe the part of the scattering curves below 


10H. S. W. Massey, Proc. Roy. Soc. 129, 616 (1930); H. S. W. Massey and C. B. O. 
Mohr, Proc. Roy. Soc. 135, 258 (1932). 

The oscillations of intensity superposed on rapidly falling atomic scattering curves, con- 
tributed by factors such as (1+sin x/x) have been considered by Mark and Wierl, principally 
in connection with their scattering of high speed electrons. (See H. Mark and R. Wierl “Die 
experimentellen und theoretischen Grundlagen der Electronenbeugung.” Gebriider Born- 
traeger, Berlin, 1931). 

2 J. H. McMillen, Phys. Rev. 36, 1034 (1931). 
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about 90°. To account qualitatively for the presence of maxima and minima 
in the scattering of electrons by atoms, Bullard and Massey™ called attention 
to the necessity for considering the effects of distortion of the electron wave, 
electron exchange and polarization. These effects have not been included in 
Massey and Mohr’s treatment of scattering by diatomic molecules. 

We are accustomed to think of a clear cut distinction between an inelastic 
and an elastic collision between electrons and atoms. In the latter case the 
electron leaves the atom unchanged, while in the former case the electron 
excites (or ionizes) the atom. The lowest excitation state differs in energy so 
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Fig. 7. Franck-Condon diagram for hydrogen. The diagrams given by Smyth and by Mulli- 
ken agree substantially in the case of all curves except the 1°S curve. Mulliken’s 1° curve is 
shown as a broken line. 
considerably from the normal state that there is no difficulty in identifying an 
electron which has rebounded from an atom elastically and one which has 
rebounded inelastically. With a molecule, however, it is conceivable that an 
electron collision could change the vibrational state and the rotational state 
of the molecule without altering the electronic state. We should therefore be 
prepared, perhaps, to find evidences of energy losses of the order of a few 
tenths of a volt. Experimentally these might cause an asymmetrical broaden- 
ing of the so-called elastic peak towards the low energy side. We have made 
as yet no special search for such an effect but it seems possible that it could 
be detected in view of the experimental work of Ramien™ on the average 
energy lost by electrons in colliding with many hydrogen molecules, and in 
view of the theoretical discussion of Massey.” 

The possible molecular changes resulting in a hydrogen molecule after an 
electron impact may be visualized by the help of a Franck-Condon diagram. 
In Fig. 7 we reproduce a somewhat simplified copy of the diagram given by 
Smyth" for hydrogen. The curve A is for the normal molecule, and the stable 
state is the one corresponding to a nuclear separation varying between / and 
N (corresponding to the lowest vibrational state). Excitation by electron 
impact is represented by a vertical displacement, within the shaded area, to 
one of the other curves. If the transition is to the 1*S state, the atoms im- 
mediately fly apart and dissociation results. If, on the other hand, the transi- 


8 E. C. Bullard and H. S. W. Massey, Proc. Roy. Soc. A133, 637 (1931). 
4H. Ramien, Zeits. f. Physik 70, 353 (1931). 

1° H.S. W. Massey, Proc. Camb. Phil. Soc. 28, 99 (1932). 

1 H. D. Smyth, Rev. Mod. Physics 3, 347 (1931). 
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tion is to the B state (i.e., to the 2'S state), we get an excited molecule which 
does not dissociate so long as it is in this state. The energy necessary to pro- 
duce either state is not sharply defined because of the finite width of the 
shaded column resting on J/.N. Thus a transition to the 1° state requires 
anything between 10.9 volts and 13.5 volts, while a transition to the 2' state 
requires an amount of energy between 11.5 and 12.2 volts and a transition to 
the 2'II state calls for something between 12.2 and 13.6 volts. This overlapping 
makes it impossible to decide whether our most probable loss—12.6 volts— 
corresponds to a transition to the 1°2 state to the 2! state or to the 2'Il state. 
It should be stated that the Franck-Condon diagram as given by Mulliken"? 
differs slightly in scale from that given by Smyth, yet in such a way as to sug- 
gest a different interpretation of our results. Briefly, if the Mulliken represen- 
tation is correct, we cannot regard our 12.6 volt loss as indicating transitions 
to either the 1° or to a low vibrational level of the 2'S state; we must attrib- 
ute it to a transition to the 2'II state, or to a high vibrational level of the 
2'S state. Dr. E. U. Condon informs us that, on account of the approxima- 
tions made in the calculations of these curves, they cannot be regarded as 
exact; thus the position of the 152 curve above 1/N (Fig. 7) may be in error 
by 2 or 3 volts. In view of this situation it is evident that the identification of 
our 12.6 volt loss with one or other of the possible transitions cannot be made 
with certainty. It may well be that our 12.6 volt peak corresponds to transi- 
tions to all three states. In view of the sharpness of the energy losses to be 
expected in afomic excitation, as compared with this blurring to be expected 
in molecular excitation, it is possible that the most probable energy loss peak 
in hydrogen (e.g., in Fig. 2) would be appreciably less steep (especially on the 
side next to the elastic peak) than the 11.6 volt peak in argon taken under 
identical conditions. Tests, however, showed that if the 12.6 volt peak in 
hydrogen and the 11.6 volt peak in argon were superposed, then the hydrogen 
peak had no more spread toward the elastic peak than had the argon peak. 
Massey and Mohr have calculated by wave mechanics the angular distri- 
bution of electrons which collide inelastically with hydrogen molecules to 
give rise to the 2'S state. These distributions have maxima at certain angles, 
e.g., for 30 volt primary electrons, a sharp maximum appears at 20°. Fig. 5 
gives the angular distribution of the 12.6 volt loss. Unfortunately the curve 
does not extend to angles smaller than 25°, but it seems obvious that there is 
no indication of the presence of a maximum at 20°. Massey and Mohr do not 
give curves for the angular distributions of electrons which have collided with 
hydrogen molecules and put them into the 1°~ state. Consequently we are 
unable to compare them with our experimental curves. There is one point 
of significance, however. By integrating over all angles Massey and Mohr 
have determined the probabilities of excitation of the 2'S state and the 1° 
state of hydrogen as a function of the energy of electron impact. The scale 
of their curves implies that the probability of excitation of the 1° transition 
is about 10,000 times as large as the probability of excitation of the 2'S 
state. If this is correct, then our 12.6 volt peak must be associated almost 
wholly with transitions to the 1°X state and the number of transitions to the 
17 R. S. Mulliken, Rev. Mod. Phys. 4, 1 (1932); see also W. Finkenburg and W. Weizel, 
Zeits. f. Physik 68, 577 (1931). 
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2'S state must be neglibly small. (This ignores, of course, the possibility of 
transitions to the 2'II state). 

In the case of scattering of electrons by argon atoms, we saw that if the 
energy lost by the colliding electrons was greater than 15.7 volts, the ioniza- 
tion potential of argon, then a second electron, which we termed the ejected 
electron, must appear. Thus the parts of the scattering curves corresponding 
to energy losses greater than 15.7 volts give us the distributions of the 
colliding, or ionizing electrons, as well as the distributions of the ejected 
electron. It was found convenient to assume that the ejected electron was 
always the slower of the two electrons moving away from the atom after an 
ionizing collision. The situation is still more complicated when we have a 
diatomic molecule like hydrogen instead of an atom like argon. When the 
collision is of the type which leaves a molecular ion, we have, as in the case of 
argon, two electrons and an ion. However, as Tate, Lozier and Bleakney'* 
have shown, some types of ionization result in dissociation with the consti- 
tuent atoms flying apart with considerable kinetic energy. We thus have, after 
collision, two ions and two electrons moving away from the place of collision, 
each with its own angle distribution characteristic of its energy and of the 
energy of the electron before collision. So far as we know, no attempt has 
been made to solve such a problem theoretically. It is evident, from the curves 
shown in Fig. 6, that, when electrons are ejected with a certain amount of 
energy from certain types of collisions, strongly preferred directions of emis- 
sion are to be found. We can also trace certain regularities in the displace- 
ment of the peak representing preferred directions of emission as the energies 
concerned are progressively altered. Thus, for a given energy of ejection, the 
position of the peak moves to smaller angles as the energy of collision is in- 
creased. On the other hand, if we keep the energy of collision unchanged we 
find that, as the energy of the ejected electron is increased the position of the 
peak moves to larger angles. There is also a regularity in the way in which the 
magnitude of the peak changes as we go across a row or down a column. A 
comparison of these results with those for argon will show that the distribu- 
tions of electrons ejected from argon atoms are similar in certain respects to 
the distributions for hydrogen. 

The results for hydrogen were obtained with the same apparatus as that 
in our work on argon. In our previous article we pointed out that the method 
of speeding up the electrons by a suitable field into the analyser introduced a 
certain amount of blurring of the sharpness of the curves which was greater 
the smaller the energy of the electrons under investigation. The same com- 
ment applies also to this work on hydrogen. 

So far as we know, there has been no attempt to work out a theory for the 
distributions of the ejected electrons. The frequent occurrence of strongly 
preferred directions of emission for these electrons suggests a peculiar kind 
of diffraction effect. 

We take great pleasure in thanking Dr. E. U. Condon for his helpful com- 
ments on the situation relating to the interpretation of the Franck-Condon 
diagram for hydrogen. 


18 W. Bleakney, Phys. Rev. 35, 1180 (1930); W. W. Lozier, Phys, Rev. 36, 1285 (1930); 
J. T. Tate and W. W. Lozier, Phys. Rev. 39, 254 (1932). 
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A description of electronic structures of polyatomic molecules in terms of atomic 
and molecular orbitals (one-electron orbital wave functions) is suggested, as follows. 
Unskared electrons are considered as occupying atomic orbitals associated with the 
various nuclei to which the electrons belong. Shared electrons are described from sev- 
eral points of view, whose simultaneous consideration should give a better understand- 
ing of their condition and functions in the molecule. In the first place, they are 
described, in terms of atomic orbitals, from the point of view of each nucleus in the 
molecule. In the second place, they are described from a unitary molecular point of 
view, in terms of molecular orbitals. It is pointed out that the chemical evidence which 
led Lewis to his concept of the electron-pair bond can now all be explained by the 
quantum theory without the necessity of using such a concept. It is noted that the 
Heitler-London, Pauling-Slater quantum-mechanical electron-pair bond concept dif- 
fers markedly from Lewis's, also that it is of more restricted application to chemical 
data. It is pointed out, in agreement with Hund, that properties of the H. and L., 
S. and P. electron-pair bond which make it useful in dealing with chemical combina- 
tion are also possessed by the concept of molecular orbitals. For example, the Pauling- 
Slater criterion that bonds which correspond to a maximum overlapping of atomic or- 
bitals are the strongest is just as characteristic of bonding molecular orbitals as of 
electron-pair bonds. An electron-pair bond is here interpreted as being little other than 
two electrons occupying a bonding molecular orbital. Or in general one must say that a 
set of » electron-pair bonds is interpreted as a set of 2” electrons occupying n molec- 
ular orbitals, because molecular orbitals are not necessarily localized between two 
nuclei like electron-pair bonds. It is concluded that the essential facts of molecular 
electronic structure can be qualitatively understood in terms of the mode of descrip- 
tion stated in the first paragraph above and of a semi-empirical valence rule (not essen- 
tially new) which summarizes, in the light of quantum theory, the most important 
regularities in regard to the types of chemical compounds which are stable. 


ATOMIC AND MOLECULAR ORBITALS FOR DESCRIBING 
MOLECULAR ELECTRONIC STRUCTURES 


1. Introduction 
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HE series of which the present paper is the second! is a development of 
a program initiated by Hund? and the writer’ and carried forward 


principally by us, by Herzberg, and by Dunkel.**°:7:5-50-12 The object 


1 A brief preliminary paper outlining some of the principles used and some of the main 


results has already been published: R. S. Mulliken, Phys. Rev. 40, 55 (1932). In this paper, 
the following changes should be made in the last paragraph on p. 58. Line 8, delete “tetrahedral- 
ized but uncombined”; lines 9 and 12, delete brackets around o, and delete “uncombined” be- 
fore 2pa; lines 10-11, delete “tetrahedralized uncombined.” Also on p. 59, line 10, delete brack- 


ets. Elsewhere, do not delete brackets or “tetrahedralized.” 





2 F. Hund, Zeits. f. Physik 51, 759 (1928). 
3 R.S. Mulliken, Phys. Rev. 32, 186, 761 (1928); 33, 730 (1929). 
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of this program has been to describe and understand molecules in terms of 
one-electron orbital wave functions of distinctly molecular character. The 
first paper of the present series! will be referred to hereafter as I. The present 
paper is to a large extent a preliminary critical review rather than a presenta- 
tion of new results. 


2. Atomic and molecular orbitals 

From here on, one-electron orbital wave functions will be referred to for 
brevity as orbitals. The method followed here will be to describe unshared 
electrons always in terms of atomic orbitals but to use molecular orbitals for 
shared electrons. This procedure was first used by Lennard-Jones for dia- 
tomic molecules, and is a partial departure from the original formulation of 
the above program. Among the main objects of the present series are the 
determination of the forms of molecular orbitals and the application of a 
knowledge of these forms to an explanation of valence and related chemical 
phenomena. 

Shared electrons include two types, namely bonding and anti-bonding 
electrons, at least in diatomic molecules; unshared electrons belong to the 
class of non-bonding electrons.*:*'°"! Anti-bonding electrons occur in dia- 
tomic molecules only when accompanied by a larger number of bonding 
electrons. Often it is hard to draw the line between shared and unshared 
electrons; in such cases, the electrons in question consist of bonding and 
anti-bonding electrons in equal numbers, at least in diatomic molecules 
(cf. third following paragraph and references 8, 9). 

By an atomic orbital is meant an orbital corresponding to the motion of 
an electron in the field of a single nucleus plus other electrons, while a molec- 
ular orbital corresponds to the motion of an electron in the field of two or 
more nuclei plus other electrons. Both atomic and molecular orbitals may be 
thought of as defined in accordance with the Hartree method of the self- 
consistent field, in order to allow so far as possible for the effects of other 
electrons than the one whose orbital is under consideration. 

Every non-degenerate orbital can be occupied by at most two electrons, 





4 F. Hund, Zeits. f. Physik 63, 719 (1930). 

°F, Hund, Zeits. f. Physik 73, 1 (1931); 74, 429 (1932). 

6 F. Hund, Zeits. f. Physik 73, 565 (1932). 

7 F. Hund, Zeits. f. Physik 74, 1 (1932). 

8 R. S. Mulliken, Chem. Rev. 9, 347 (1931). On p. 351 (foot), p. 353 (middle) and p. 355 
(middle), the writer gives the incorrect impression that, omitting the energy of repulsion of the 
nuclei, the promoted wave-function 2po has a higher energy than the original 1s atomic wave- 
function. Actually the energy is a little lower (for R>O) than that of 1s in the H atom. This, 
however, does not in any essential way affect the validity of the arguments given. The fact that 
the energy of 20 goes up sufficiently rapidly with reference to that of 1se, in particular the fact 
that the energy of 2/0 plus the nuclear repulsion energy causes a net repulsion of H and H’*, is 
sufficient. 

* R.S. Mulliken, Rev. Mod. Phys. 4, 1 (1932). 

10 G, Herzberg, Zeits. f. Physik 57, 601 (1929). 

" G. Herzberg, “Molekiilstruktur,” Leipziger Vortriige, 1931, p. 167. S. Hirzel, Leipzig. 

® J. E. Lennard-Jones, Trans. Faraday Soc. 25, 668 (1929). 
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corresponding to the two possible orientations of the electron spin. Orbitals 
having n-fold degeneracy can be occupied by at most 2” electrons. Examples 
of non-degenerate orbitals are s orbitals in atoms, ¢ orbitals in diatomic or 
linear molecules’ and most orbitals in non-linear polyatomic molecules. 
Examples of degenerate orbitals are p, d,... orbitals in atoms (2/+1-fold 
degeneracy), 7, 6, . . . orbitals in diatomic or linear molecules (2-fold degen- 
eracy), [7] orbitals (2-fold degeneracy), in molecules like NH3, NOs;-, ete. 
(cf. 1), d, (2-fold degeneracy) d, and p (3-fold degeneracy) in molecules 
having tetrahedral or octahedral symmetry (cf. I). 

It should be noted that the use of atomic orbitals for describing the con- 
dition of unshared electrons in molecules often gives a false impression of the 
amount of degeneracy. For example if one writes 1s*1s*02s* for the electron 
configuration® of Lis, the description suggests only one ionization energy for 
the 1s electrons, since 1s orbitals on two Li atoms are identical in energy. If 
one writes o1s*o0*1s*o2s*, however, using molecular orbitals exclusively,® it is 
made explicit that there are two distinct energies corresponding to ols and 
o*1s. In this particular example the energy difference is negligible, but in 
many cases, especially where there are electrons which are near the border 
line between unshared and shared, considerable energy-splittings may exist. 

For example in Oz one may perhaps best write 1s*1s?2s*2s* . . . but must 
then grant that the 2s type, although essentially unshared, has a markedly 
double ionization energy.’ In Ne, the sharing of the 2s electrons is so strong 
that it is best to write 1s*1s°¢2s*¢*2s*, the types 2s and o*2s being very dif- 
ferent in energy.® In CF,, one might write { 1s?2s?2p[z ]}*} -*} 12} -a°b", where 
a and b represent the shared electrons. In so doing, it should be recognized 
that four distinct energy values are to be expected for the unshared 1s fluorine 
electrons, likewise for the 2s and for the 2p[z],—or possibly eight values for 
2p|z| because of a possible splitting up of the degeneracy implied by the 
symbol [7]. The four-fold splitting might be appreciable for the 2p[7] and 
possibly also for the 2s, although of course negligible for the 1s. The carbon 
1s is of course very different in energy from the fluorine 1s. These examples 
should be sufficient to show how to guard against a possible misinterpretation 
of the use of atomic orbitals in describing molecules. 

In the present method, molecular orbitals are conceived of as entities 
quite independent of atomic orbitals. Nevertheless in practise molecular 
orbitals can usually be conveniently approximated by building up linear 
combinations of orbitals of the atomic type. The present method of thinking 
in terms of the finished molecule, used already by Lewis in his valence theory, 
avoids the disputes and ambiguities, or the necessity of using complicated 
linear combinations, which arise if one thinks of molecules as composed of 
definite atoms or ions. 


3. One-nucleus and other viewpoints 


Understanding of the electronic structure of molecules is greatly aided 
by introducing a set of partially overlapping descriptions whereby the elec- 
trons immediately surrounding each nucleus in a molecule are described in 
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terms of a set of atomic orbitals corresponding to the viewpoint of that 
nucleus. (Sidgwick,—reference 13, p. 98,—has used essentially this view- 
point.) Such one-nucleus viewpoint atomic orbitals must be taken as consider- 
ably deformed compared with ordinary atomic orbitals, the deformation be- 
ing thought of as caused by the strong fields of the other nuclei. Unshared 
electrons, of course, are here assigned exclusively each to a particular nucleus, 
so that the one-nucleus viewpoints for them give the same description as that 
introduced in the first paragraph of section 2. Shared electrons, however, 
being considered as belonging to two or more nuclei, receive non-mutually- 
exclusive descriptions in terms of deformed atomic orbitals associated with 
each, or at least with several, of these nuclei (cf. descriptions of N.*, HCI, and 
PtCl,= in section 8 below and of CF, and other molecules in I as examples). 

An understanding of molecular electronic structures can not infrequently 
be advanced by using still other auxiliary viewpoints. The well-known 
united-atom viewpoint is sometimes useful (cf. reference 9, p. 19). In hydrides 
especially, this or, rather, a viewpoint in which each hydrogen nucleus is at 
first regarded as united with a neighboring nucleus of larger charge, is valu- 
able. This amounts to saying that in hydrogen compounds, the hydrogen 
one-nucleus viewpoints can often appropriately be treated as unimportant 
compared with the viewpoints of other nuclei in the molecule or that the 
hydrogen nuclei can be regarded merely as perturbing force-centers. The 
validity of this point of view is indicated by a number of facts concerning 
diatomic hydride band spectra,’ also by chemical data, notably the behavior 
of the boron hydrides.* 


The familiar chemical method of regarding many molecules as built up 
of “radicals” suggests the usefulness of a point of view making use of shared 
orbitals belonging to radicals instead of, or in addition to, the viewpoint of 
molecular shared orbitals. 


QUANTUM THEORY OF VALENCE 


4. Review of theories of valence and molecular structure 


A semi-historical survey of some of the theories dealing with valence and 
molecular structure will put the whole problem in better perspective. Inci- 
dentally, the writer hopes to show that there are no compelling reasons, either 
empirical or theoretical, for placing primary emphasis on electron pairs in 
constructing theories of valence. He hopes thereby to remove possible ob- 
jections to the present method based on its lack of such emphasis. 

The best chemical theory of valence cevering all types of compounds is 
generally agreed to be that developed principally by G. N. Lewis." To a 
rather large extent, the essential features of this theory still stand, although 
their meaning has been made clearer and more specific by interpreting them 
in the light of the quantum theory. The most important features of Lewis’s 
theory are perhaps the following: 


13 G. N. Lewis, Valence and The Structure of Atoms and Molecules. The Chemical Catalog 
Co., New York, 1923. Sidgwick’s excellent book should also be consulted: N. V. Sidgwick, The 
Electronic Theory of Valency. The Clarendon Press, Oxford, 1927. 
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(1)(a) Every atom tends so to give, take, or share electrons as to be 
surrounded by an outer “group of eight” electrons (reference 13, p. 79),— 
except that for atoms such as H, He, Li this is replaced by a group of two. 
(6) Electrons which are shared may be shared equally or unequally by 
atoms, so that one can readily account for all intermediate stages between 
homopolar molecules like Hz or Cl, and polar molecules such as NaCl. (c) 
Besides the cases having a stable group of two, already mentioned, Lewis 
noted the existence of other probable exceptions to the “rule of eight,” 
namely cases like those of PCl;, SF., PtCle= where a central atom appears 
to share more than eight electrons with other atoms (reference 13, pp. 102, 
114 et seg.). Lewis suggested that in such compounds, the shared electrons in 
excess of eight have “passed into a secondary valence shell” of higher energy. 

(2)(a) The single bond of the old valence theory of organic chemistry is 
interpreted as a pair of electrons held jointly in the outer shells of two atoms, 
and this concept of the chemical bond is extended to inorganic compounds. 
This electron-pair bond may be symmetrically shared as in He or Cle, but 
more often the electron-pair is nearer one atom than the other. (b) The 
valence of an atom in any molecule is defined as the number of electron pairs 
which it shares with other atoms (reference 13, p. 104). This is Langmuir’s 
“covalence.” Lewis's definition of valence makes the numerical valence the 
same as the “coordination number” of Werner (reference 13, p. 114). Using 
this definition, atoms which form a “group of eight” (Langmuir’s “octet”) by 
sharing have a valence of four, and “we may regard the maximum valence 
of four as a sort of norm,” although atoms sharing more than four pairs have 
a higher valence, e.g., S or Ptin SF, or PtCl,= has a valence of six. When there 
is an outer shell of eight electrons, there is a strong tendency for these to be 
all shared. In dealing with ionic or ionized molecules, Lewis uses the term 
“polar number” (reference 13, pp. 70, 104), e.g., in Cot*Cl>e, the Co atom 
is said to have a polar number +2, the Cl a polar number —1. This is 
Langmuir’s “electrovalence.” (c) Usually each atomic partner in an electron- 
pair bond furnishes one of the two electrons, as e.g. in He, CHy, CCl,, but 
very often one partner furnishes both: examples, BR;- NH; and other am- 
monia complexes, where the N atom furnishes both electrons of a pair; 
SO,-, ClOg-, etc., where the central atom furnishes all the shared electrons 
except the ionic charges. One-sided sharing as here and in (a) causes polarity 
in the molecule, but does not require a really ionic conception of valence, such 
as Kossel has used," except in extreme cases. (d) Double and triple bonds 
are considered to involve sharing of two and three pairs of electrons, but 
are considered to be something very different from merely two or three 
ordinary electron-pair bonds. Multiple bonds are apparently rarely formed 
by atoms other than C, N, and O. 

The fundamental ideas of Lewis’s theory are perhaps the three following: 
(A) each atom (better, each nucleus) in a molecule tends to become sur- 


4 W. Kossel, Ann. d. Physik 49, 229 (1916); Naturwiss. 7, 339, 360 (1919); cf. A. E. van 
Arkel and J. H. de Boer, Chemische Bindung als elektrostatische Erscheinung, S. Hirzel, 
Leipzig, 1931 (German edition). 











54 ROBERT S. MULLIKEN 


rounded by a set of closed shells of electrons; (B) shared electrons, forming 
chemical bonds of the homopolar type, are usually localized between two 


atoms (or nuclei) which they link together; (C) a chemical bond usually con- 
sists of a pair of electrons somehow rather closely united to each other. Of 
these three ideas, Lewis seems to have considered (C) as the most funda- 
mental (cf. reference 13, pp. 79-81). Ideas (A) and (B) were already in use 
before Lewis's 1916 paper, although Lewis’s work greatly developed them 
and their application. Lewis’ most important contribution to the theory con- 
sisted perhaps in his introduction of the idea of the completion of closed 
groups by sharings. Idea (A), limited, however, to the “group of eight,” was 
proposed by Abegg in 1904 for polar valence and was also used in the 
same way by Lewis and later by Parson (1915) and Kossel (1916). Idea (B) 
was adapted by Lewis from the concepts of organic chemistry. Langmuir," 
in further developing the subject in 1919, emphasized idea (A) and spoke 
of the “octet theory” of valence. Bury“ in 1921 also emphasized idea (A). 
With the development of the quantum theory, a more fundamental 
theoretical background has gradually been created for the originally largely 
empirical ideas of the Lewis theory. Bohr’s theory of the periodic system, as 
modified by Stoner and Main Smith, shows why idea (A) is important in 
valence theory. The group of two found in H~-, He, Li* is 1s? and presumably 
the shared pair of electrons in Hg is essentially the same, while the highly 


‘ 


stable “octet” is ns°np®. Bohr’s theory also provides for less stable groups 
of 18 and 32 electrons. Sidgwick” showed in a general way in 1923 on the 
basis of the original Bohr theory how a group of 18, containing ws, np and 
nd electrons,can be used to interpret complex ions like PtCl.= or Co(NH3)¢***. 
Revised interpretations which regard the shared electrons in such compounds 
as belonging to nd, (n+1)s, and (n+1)p, have been given later (cf. Pauling," 
also I). Knorr!’ has discussed the interpretation of Lewis's theory from the 
standpoint of the modified Bohr theory. 

Bohr’s theory was unable to give a satisfactory understanding of the 
sharing of electrons in molecules, but the new quantum mechanics is showing 
itself capable of doing this. 

London and Heitler, generalizing results obtained from a quantum- 
theoretical study of the formation of He from H+H, attempted to construct 
a valence theory which has often been supposed to be the quantum-mechani- 
cal equivalent of Lewis’s, and which emphasizes Lewis’s ideas (B) and es- 
pecially (C). This so-called spin theory of valence emphasizes the pairing 
of electrons and their spins, but deals primarily with the interactions of 
atoms as wholes. It has, however, not proved very successful.>:*"° London 
made the suggestion that the shared electrons in excess of eight (ms? np*) in 


6 T, Langmuir, J. Am. Chem. Soc. 41, 868 (1919); Science 53, 290; 54, 59 (1921). 

1 C,R. Bury, J. Am. Chem. Soc. 43, 1602 (1921). 

17 N. V. Sidgwick, Trans. Chem. Soc. 123, 725 (1923). 

18 L. Pauling, J. Am. Chem. Soc. 53, 1367, 3225 (1931); 54, 988 (1932). 

19 C, A. Knorr, Zeits. f. anorg. allgem. Chem. 129, 109 (1923). Also W. A. Noyes (1917), 
and others: cf. references given by W. A. Noyes, Zeits. angew. Chem. 44, 893 (1931). 




























































ELECTRONIC STRUCTURES OF MOLECULES 


PCl;, SF, and the like are nd electrons; this idea is more specific than the 
similar idea of Lewis (cf. item 1 (c) above), and seems to be a good one.'* 
One should distinguish carefully* between Heitler and London's valence 
theory and their valuable perturbation-method for calculating energies of 
molecule-formation. 

Pauling'* and Slater?® have generalized the Heitler and London results 
for He in a different way, considering more specifically the interactions of 
pairs of electrons, one electron from each of two atoms. In this way they have 
obtained important results on valence and valence angles in polyatomic 
molecules. Their work emphasizes Lewis’s idea (C) in modified form but 
focuses attention more on idea (B). Their theory seems more limited than 
Lewis’s in that it requires that the two electrons forming an electron-pair 
bond necessarily come from two different atoms. This makes it inapplicable 
to many chemical compounds (cf. summary of Lewis’s theory, item 2c) un- 
less one assumes them to be formed from ions, a procedure which in many 
cases seems decidedly artificial. The theories of Heitler and London, Pauling 
and Slater might be called electron-pairing theories if Lewis's is called an 
electron-pair theory. It should also be pointed out that the H.L.P.S. electron- 
pair differs very considerably from Lewis’s conception of the electron-pair 
bond in that the electrons are much less closely associated (see below); in 
this respect it approaches the truth much more closely than does Lewis’s 
conception.—Pauling and Slater consider a double bond to be merely two 
ordinary single bonds sticking out from each atom in different directions, 
and treat the triple bond in a similar way. In this they do not agree very well 
with Lewis (cf. summary, item 2d), nor do they agree with results obtained 
from the method of molecular orbitals (see below). 

Dunkel*! emphasizes Lewis's ideas (A) and (B), and describes the shared 
electrons in molecules in terms of the symbols [¢] and |], implying orbital 
wave functions having properties similar to those of o and 7 orbitals in 
diatomic molecules. The symbol |o], for instance, indicates a molecular 
orbital concentrated in the region between two nuclei and roughly sym- 
metrical around the line joining them. This is justified by the fact that the 
shared electron moves in a Hartree field which is roughly symmetrical around 
this line. Although shared electrons of any given kind nearly always occur 
in pairs in stable molecules, this seems to be incidental in Dunkel’s as in the 
present work. Unshared electrons are classified as they would be in free atoms. 
Dunkel’s viewpoint evidently resembles that of the present paper. Hund in 
his recent papers®“ has used Dunkel’s [co], [7] classification, but has used 
the simple symbols ¢, 7. 

Hund in his important recent papers’ *:? concludes that the results ob- 
tained by Slater and Pauling using the electron-pair bond method can also 
be obtained, sometimes more easily, by the method of molecular orbitals, and 
that the two methods are in many respects equivalent if one restricts the 
use of molecular orbitals to a type localized between two nuclei,—in agree- 


20 J.C. Slater, Phys. Rev. 37, 481; 38, 325, 1109 (1931). 
2 M. Dunkel, Zeits. f. phys. Chem. [ B) 7, 81; 10, 434 (1930). 








56 ROBERT S. MULLIKEN 


ment with Lewis’s idea (B). Reference should be made to Hund’s papers for 
a critical discussion of many points concerning the various quantum theories 
of valence. 


Hund has classified valences according to a number of types, depending 
in part on whether they are formed by s or p electrons or involve partial 
hybridization of s and p (Hund’s g type). Hund also classifies chemical bonds 
according to a number of types, and concludes among other things that single 
bonds are always of the type [a]?, double bonds of the type [a]? [z]?, triple 
bonds of the type |o|? [x]. It will be noted that this description of double 
and triple bonds agrees better with Lewis’s conclusions (cf. summary, item 
2d) than does the description given by Slater and Pauling. The absence of 
free rotation about double bonds can be well accounted for’:'’ by the form- 
ulation [o}? |z]*, which really goes back to Hiickel,” while the validity 
of the explanation given by Pauling and Slater seems doubtful. 

Hund shows’ that the familiar rules of organic chemistry can be explained 
in terms of the quantum-mechanically predicted behavior of some of the 
possible types of valence and of bonds. He shows that not al! of the theoreti- 
cally possible types are found in ordinary organic compounds, and shows 
how this can be understood in terms of energy relations; in particular, by 
assuming that of the types [a | and [7] corresponding to a p valence electron, 
|o | in practise nearly always has the lower energy. 

In all this work, Hund’s purpose has been to show how the ordinary 
rules of valence can be derived from the principles of the quantum theory. 
In doing this, he has found it necessary to specialize the concept of molecular 
orbitals to the case of orbitals localized between two nuclei. Such a localized 
molecular orbital occupied by two electrons he regards as corresponding to 
the valence bond of organic chemistry or to that of Lewis, and as being 
essentially equivalent in the case of ordinary stable molecules to the electron- 
pair bond of H.L.S.P. 

Hund remarks incidentally that the use of localized molecular orbitals 
is a somewhat poorer approximation than the use of the non-localized orbitals 
which are advocated in the present series, and illustrates this’ by discussion 
of a case which corresponds to CH, or H2Q. He also points out that in some 
cases localized orbitals cannot be used at all, and illustrates this by a con- 
sideration of crystal lattices and of the benzene molecule CgH¢ and related 
molecules. In his plausible and important interpretation of CsH«, Hund uses 
the same point of view that is emphasized in the present series. His non- 
localized “x” type in CsHe¢ is very similar to the central-nucleus-viewpoint 
type 2p{o] in COs" or NO3°~ if this type is assumed to be filled by two elec- 
trons shared by the O atoms in addition to their ordinarily assumed sharing 
of one pair each (cf. I, p. 61). 

In the present series, idea (A) of Lewis’s theory is emphasized, idea (B) 
is adopted in generalized form (non-localized molecular orbitals), while idea 
(C) is considered as corresponding more to an incidental than to a really 
essential characteristic of chemical combination. While the pairing of elec- 


# E. Hiickel, Zeits. f. Physik 60, 423 (1930). 
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trons, which is a consequence of the Pauli principle combined with the fact 
that the electron spin can have two and only two orientations, is very im- 
portant for the general theory of atomic and molecular structure, especially in 
connection with the formation of closed shells, the importance of the specific 
role of electron-pairs or of electron-pairing in the formation of molecules has, 
in the writer’s opinion, been somewhat over-emphasized, first in Lewis's 
theory, then in the work of Heitler and London, Pauling and Slater. Reasons 
for this opinion have been given in a previous paper dealing mainly with 
diatomic molecules.* It seems desirable, however, to take the matter up 
again in some detail here, approaching it historically. 


5. Analysis of the origin and development of the concept of the electron- 
pair bond 


At the outset it is important to notice that Lewis’s development of the 
concept of the electron-pair bond preceded our present knowledge of the 
electron spin and its properties and of the Pauli principle. 

In developing his theory, Lewis noticed that while nearly half of the 
known atoms contain an odd number of electrons, very nearly all stable 
molecules contain an even number. This and other information on atomic 
structure and on molecules suggested that electrons have a tendency to 
form stable groups of two, held close to each other (cf. reference 13, p. 82 
and Fig. 22) by strong forces of unknown character. These forces Lewis 
for some time considered to be magnetic, since the paramagnetism char- 
acteristic of atoms or molecules with an odd number of electrons usually 
disappears when they combine to give molecules with an even number. 
Lewis considered that in the formation of a molecule from two atoms con- 
taining unpaired electrons, the electrons attract each other in pairs, and that 
this attraction is, to a considerable extent at least, the cause of molecule 
formation (cf. e.g., reference 13, p. 149, second paragraph). 

The subsequent development of the quantum theory has provided ade- 
quate explanation for the above facts concerning electron-pairs. The fact 
that most of the electrons in atoms as well as in molecules are paired is now 
well understood in terms of Pauli principle and electron spin. The fact that 
paramagnetism goes with unpaired electrons, diamagnetism with paired 
electrons, is also well understood. The fact that electrons are more often 
unpaired in atoms than in molecules can be understood in terms of the Pauli 
principle, electron spin, and orbital degeneracy (see below). 

According to the quantum theory an electron-pair, in a molecule or in 
an atom, consists of two electrons, with spins antiparallel, occupying the 
same orbital and “symmetrically related” so far as their orbital motion is 
concerned.* The magnetic forces between the two electrons are very small 
while the electrostatic repulsion is large. It is true that the symmetrical rela- 
tion keeps the electrons nearer together than they would otherwise be, and 
in this respect shows a similarity to Lewis’s concept of a pair of electrons 
held together by a strong attraction. But, directly contrary to Lewis's idea, 
the increased nearness of the electrons usually only increases their energy 
of repulsion. 
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Thus the formation of a pair of electrons when a molecule is formed, 
i.e., the occupation of a shared molecular orbital by fwo electrons, far from 
being per se a cause of combination of two atoms initially having each an un- 
paired electron, should tend to prevent combination. According to the 
present point of view the two electrons concerned act as bonding electrons 
in spite of considerable repulsive forces acting between them. Further details 
will be found in a discussion of Het and Hein reference 8. 

Why then, are unpaired electrons so very much rarer in stable chemical 
molecules than in atoms? The answer involves several distinct considerations. 
In the first place, every atom having an odd number of electrons has neces- 
sarily at least one unpaired, while as soon as the formation of molecules is 
permitted, this necessity disappears. This leaves three questions. (1) Why 
do atoms often have more unpaired electrons than the minimum possible 
number, which is zero for an even total number of electrons, and one for an 
odd total number of electrons? (2) Why do stable molecules nearly always 
have an even number of electrons? (3) Why do such molecules nearly always 
have all their electrons paired? 

The answers to all these questions depend largely on the occurrence and 
behavior of non-degenerate and of degenerate orbitals. If every orbital 
in atoms and molecules were decidedly different from every other in energy, 
then no atom or molecule with an even number of electrons would contain 
unpaired electrons when in its normal state. The electrons would settle 
down into the orbitals of lowest energy, with two electrons in each such 
orbital. 

But when there is orbital degeneracy, i.e., when two or more orbitals 
form a group whose members are equal in energy (cf. section 2 above), then 
the energy of the electron system as a whole is a minimum if every electron 
occupies a different orbital belonging to the degenerate group. For example 
in the nitrogen atom the energy is a minimum when we have the electron 
configuration 1s? 2s? 2p* and the state 4S. The orbitals 1s and 2s are non- 
degenerate and each is occupied by two electrons, but the type 2p has a 
three-fold degeneracy, being capable for example of giving rise to three dis- 
tinct orbitals, 26,:, 20, and 2p_; in a strong magnetic field. The state 4S 
corresponds to one electron, unpaired, in each of the three orbitals 2p,:, 
2/0, and 2p_,. The energy is a minimum for this state because it permits the 
orbitals of the three electrons to be antisymmetrically related. The re- 
sultant spin S(S=13 here) is always equal, in molecules as well as in atoms, 
to 3 times the number of unpaired or antisymmetrically related electrons, 
here three. If two electrons occupied the same orbital, e.g., 2p_1, they would 
necessarily be symmetrically related; but a symmetrical relation, as already 
noted, forces the electrons together, and would raise the energy above that 
of the ‘S normal state. 

It remains to give reasons why stable molecules nearly always have an 
even number of electrons, in pairs. In the writers’ opinion, this can best be 
explained in terms of Lewis’s idea (A) according to which each nucleus in a 
molecule tends to become surrounded by an atomic closed shell of electrons. 
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Since such a closed shell, e.g., 2s? or ns*np*®, contains an even number of 
electrons,—for the reason that it is in the last analysis made up of pairs of 
electrons, although these often belong to degenerate groups, e.g., p*,—every 
molecule which satisfies Lewis’s idea A must automatically contain an even 
number of electrons, all paired. The quantum-mechanical background of idea 
A will be taken up below. 

A further reason why molecules with an even number of electrons much 
more rarely contain unpaired electrons than atoms with an even number 
of electrons is to be found in the fact that orbital degeneracy is less usual 
in molecules than in atoms (cf. section 2, above). The Os molecule is a 
good example of an “exception that proves the rule.” Here the electron 
configuration is of a type - - - 7’, and the normal state is *2~, with S=1 
and so two unpaired electrons. The 7 type of molecular orbital has a two- 
fold degeneracy, so that with two 7 electrons, one can be z+, the other z~, 
and the two can be antisymmetrically related and so give rise to the normal 
state with S=1. 

From the foregoing discussion it will be seen that the empirical evidence 
which led Lewis to emphasize his idea (C), according to which the formation 
of electron-pairs has a peculiar importance for chemical phenomena, can 
now all be explained satisfactorily on the basis of general quantum-theoretical 
considerations without any necessity of adopting idea (C). 

Nevertheless it must be admitted that most of the striking chemical facts 
are not inconsistent with an electron-pair bond theory of homopolar valence 
similar to Lewis’s. (His idea of a special attraction between electrons as the 
cause of pairing must of course be dropped.) A pair of electrons occupying 
equivalent localized molecular orbitals appears to be on the whole the near- 
est quantum analogue of Lewis’s bond. For the best-known stable molecules 
the H.L.S.P. electron-pair,—cf. Hund (section 4), and section 13,—is equally 
good or perhaps better, but for molecules where a pair of electrons is fur- 
nished by one atom (cf. Lewis’s theory, item 2c, and section 6), H.L.S.P. elec- 
tron-pairs are less suitable than pairs of bonding orbitals. On the other hand 
for describing molecules having unusually loose binding, or molecules in a 
high state of vibration or in process of dissociation, provided dissociation 
would cause the breaking up of electron-pairs into unpaired electrons, elec- 
tron-pair bonds of the H.L.S.P. type should often constitute a much better 
approximation than bonding orbitals. 


6. Advantages of the present method 

The concept of the bonding molecular orbital has the following ad- 
vantages over the Heitler-London, Slater-Pauling electron-pair bond. (1) 
It is not necessary to have two electrons to get a chemical bond; one electron 
occupying a bonding orbital has a bonding effect, although of course not 
as strong as if two were present.® Such “one-electron bonds” are rarely found 
in stable molecules, it is true, but the concept of bonding orbitals which 
explains them as a special case can lay claim to greater generality than the 
electron-pair bond concept which has to be replaced by a different special 
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concept in order to deal with them. The importance of one-electron bonds'® 
is greatly increased if one ventures into the realm of spectroscopically ob- 
servable molecular entities and of intermediate products in chemical re- 
actions. (2) Bonding molecular orbitals can be constructed for any degree of 
polarity or unequalness of sharing of electrons between nuclei. H.L.S.P. 
electron-pair bonds, however, cannot be used for molecules of the extreme 
polar type, and intermediate cases can be taken care of only by forming linear 
combinations of wave functions of the polar and electron-pair bond types 
(cf. section 13). In cases like Cu(NH3),** or CO, the H.L.P.S. method can 
be used only by starting with assumed states of ionization like Cu> and NH3;* 
or C~ and ©*, while the method of molecular orbitals can start with Cu*+ 
and NH; or C and ©. For example, the pair of electrons which any N atom 
in NHs shares with the Cu or Pt atom in Cu( NH3)4**+ or Pt(NH3),.**** ean 
be assigned to a symmetrically related pair of molecular orbitals which at 
first belong wholly to the NH; but take on more and more of the character of 
Cu or Pt atomic orbitals as the NH; approaches the metal atom. Here the 
method of molecular orbitals is superior even for loosely bound molecules. 
These various possibilities make the concept of molecular orbitals more gen- 
erally useful than that of H.L.S.P. electron-pair bonds. (3) Bonding mole- 
cular orbitals are not restricted, like electron-pair bonds, to holding just two 
nuclei together, but may be distributed between several nuclei. This makes 
the concept a much more flexible one than that of the electron-pair bond, 
and makes it possible to account rather directly for a number of phenomena, 
many of them more or less spectroscopic but others definitely chemical in the 
ordinary sense, which the electron-pair bond method does not touch. In par- 
ticular, molecular orbitals can be chosen in conformity to the actual sym- 
metry of the nuclear arrangement, while the electron-pair bond method often 
disregards (even though it often predicts) this. 

In general it may be said that there are many phenomena which can 
be interpreted in terms of electron-pair bonds only if after setting up these 
bonds, various linear combinations are formed, while the molecular orbital 
concept goes more directly at the solution, although often seemingly neglect- 
ing certain features expressed by the electron-pair concept. It appears prob- 
able that in practise one can expect most of the phenomena expressible by 
the special concept of electron-pair bonds to drop out of the application of 
the less specialized method using molecular orbitals. This should become 
clear in the detailed discussion of examples in later papers. It is quite possible, 
however, that the electron-pair bond method may for many problems be 
more adapted to quantitative calculations than the present method. 


7. Statement and justification of a valence rule 


Most of the ordinary numerical aspects of valence, alike for homopolar, 
heteropolar and for intermediate types of compounds, appear to be expressi- 
ble by a simple valence rule (cf. I): Every nucleus in a molecule tends to be 
surrounded, by means of sharing or transfer of electrons, by an electron dis- 
tribution corresponding to some stable configuration having a total charge ap- 
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proximately equal to or somewhat exceeding the charge of the nucleus. This is 
essentially idea (A) of Lewis’s theory together with a generalization of his 
idea (B). This valence rule, it will be noted, is expressed in terms of one- 
nucleus viewpoints, according to which the state of the electrons near each 
nucleus is described by means of a set of electron quantum numbers or atomic 
orbitals associated with the viewpoint of that nucleus. 

The important tendency of many atomic nuclei to have all their outer 
electrons shared (cf. section 4, item 2), as e.g. in the formation of NH,*, 
BR3;:-NHs3, or Cu(NHs3),** from NH; is not accounted for by this valence 
rule, but its explanation appears to follow incidentally from the application 
of the present method (cf. discussion of NH; and CH;in I). 

By “stable configuration” in the valence rule is usually meant a set of 
atomic orbitals completely occupied by electrons (i.e., a set of closed shells) 
and of such type that further electrons could go only into orbitals of dis- 
tinctly higher energy. Stable configurations are usually sharply marked off 
from other configurations with either more or fewer electrons, by virtue of 
large energy changes that go with a change in principal quantum number, 
or with the change from a penetrating to a non-penetrating type of orbit 
even without change in principal quantum number. 

Stable configurations which are sharply defined in respect to energy tend 
to give sharply-defined valence relations. In the transition groups where a d 
shell is in the process of being built up, there is at first no limiting sharply- 
defined stable configuration, and highly variable valence is the result. As the 
d shell approaches completion, however, the complete group d'’s*p® appar- 
ently begins to serve as a stable limit, which is approached or reached in a 
number of complex molecules. 

A quantum-theoretical justification of the above valence rule follows in 
part from the fact that, for the electrons in the neighborhood of every nucleus 
in a molecule, the Pauli exclusion principle makes essentially the same re- 
quirements in regard to quantum numbers and closed shells as for the elec- 
trons in an isolated atom. This is obviously true for the unshared, inner 
electrons, and for the ionic structures in definitely heteropolar valence, while 
its truth for shared electrons appears plausible but requires further investi- 
gation. The best evidence for its correctness for shared electrons seems to be 
the empirical evidence of the success of the rule in interpreting chemical facts. 
Dunkel* has justified his (more or less tacit) use of a similar rule by similar 
arguments (cf. also Sidgwick, reference 13, p. 98 et seq.). 

As a result of sharing, the numbers of electrons surrounding certain nuclei 
sometimes very considerably exceed the numbers in the corresponding neutral 
atoms. In such cases the sharing is usually relatively loose and presumably 
is strongly one-sided, although not necessarily ionic. Thus, for example, the 
electrons shared by the Pt or Cu nucleus in Pt(NHs3)¢**** or especially 
Cu(NHs3),** doubtless belong much more to the N nuclei than to the metal 
nucleus. The shared molecular orbitals would be approximated by linear 
combinations containing a relatively small proportion of an orbital,—nd, 
(n+1)s, or (n+1)p,—of the metal nucleus and large proportions of 2p¢ 
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orbitals of the N nucleus. Nevertheless the magnetic and other properties 
of such molecules'* indicate that even such one-sided sharing is commonly 
(but not always) strong enough so that, as far as the Pauli principle is con- 
cerned, part or all the metal-nucleus-viewpoint orbitals nd, (7+1)s, (u+1)p 
are effectively filled. 

In general, the ability of a nucleus to hold more electrons, but a rather 
definitely limited number more, when in a molecule than in an atom, appears 
to be conditioned by the fact that shared electrons are, so to speak, electro- 
statically on a part-time basis for each nucleus while with respect to the Pauli 
principle they are serving full time. 

Proceeding further with the quantum-mechanical justification of the 
valence rule, one needs next to show why an energy-decrease should occur 
when atoms or ions so combine that each nucleus becomes surrounded by a 
“stable configuration.” The heteropolar case is fairly well understood. A 
qualitative understanding also of the energy decrease which occurs in homo- 
polar sharing seems to be obtainable in terms of one-nucleus viewpoints. 


8. Cause of homopolar valence forces 

London and Heitler concluded that the previously mysterious homo- 
polar valence forces are explained by “resonance” or “exchange” integrals,— 
like Jz in Eqs. (2), (3) below and similar integrals occurring in the electron- 
pair case,—which give rise to sharp attractions or repulsions. As will be seen 
in section 10, these energy integrals correspond to the fact that bonding 
molecular orbitals give a higher, anti-bonding orbitals a lower, electron 
density in the regions between nuclei than if they were formed by mere over- 
lapping of the electron densities of the atomic orbitals from which they 
might be formed. In a certain sense the above simple explanation of the 
homopolar valence forces is an adequate one. But it is of interest to see if 
one can go farther in understanding the matter physically by seeking to find 
reasons why these changes in electron density and in energy should occur 
when molecular orbitals are formed. 

The simple case of H2* is instructive.* If we let E.E. represent the energy 
of the electron, taken as zero for R=, and N.E. the energy of repulsion of 
the nuclei, then the energy change as the two H nuclei come together is 


AE = N.E. + ELE. 


If E.E. decreases with decreasing R considerably faster than N.E. in- 
creases, until R reaches a fairly small equilibrium value, we have a stable 
molecule (negative AZ). If it increases, or decreases less rapidly than N.E. 
increases, we have repulsion. The two states o1s and o*1s of H2* correspond 
to these two cases. The exchange energies in the two cases are proportional 
respectively to —J2:and +s. 

The fact that AE is negative for o1s and positive for o*1s corresponds to 
a rapid decrease of E.E. with R for ¢1s and a very much less rapid decrease 
for o*1s. A good physical reason for this can be seen when we note that as 
R-0, the 1s atomic orbital of H+H+ must shrink through o1s to a much 
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more concentrated 1s orbital of the united-atom, while through o*1s one 
reaches 2po0 of the united-atom, which is of about the same degree of con- 
centration and has the same energy as 1s of H+H?*. 

These considerations suggest that the behavior of bonding molecular 
orbitals like o1s of Ht is, qualitatively at least, essentially the result of an 
increase in effective nuclear charge without change of quantum numbers (no 
promotion), as R decreases. This effect must of course be intense enough to 
make AE negative in the above equation over a considerable range of R 
values. The fact that it actually does so for o1s of H,+ could hardly have been 
predicted from our qualitative explanation, but this does not render the latter 
invalid, although it does show the superiority of the Heitler-London method 
for quantitative prediction. With anti-bonding orbitals like o*1s, the changed 
quantum numbers (promotion) more or less neutralize the effect of the 
increasing effective nuclear charge, thus permitting N.E. always to exceed 
the negative of E.E. in the above equation. 

A somewhat more complicated case than H.* is that of the N2* molecule 
(cf. also section 10). The formation of the molecule in its normal state may be 
expressed as follows: 


N(1s22s22p3) + N+(1s22s22p2)  Ngt(1s?1s2e2s2o*2s?x2p'o2p). 


If we let R= x—R=R,.—R=0 (N+Nt—-N,t-Si*), we should probably 
have? 


N [2s] — N2*[o2s] — Sit[2s], N[2s] > Ny*+[o*2s] — Sit[3p(o) J, 
N[2p(c) ] > Net [o2p] — Sit [3s], N[2p(r)] — Not [2p] — Sit[2p(x) J. 


Here the bonding molecular orbitals ¢2s and 72p behave like ols of H2* in 
that they are unpromoted as R-0, while the anti-bonding orbital behaves 
like o*1s of H2* in that it is promoted as R-0. But ¢2), which is promoted as 
R-0, is a bonding orbital, like 2p which is unpromoted. These results in- 
dicate that the difference between bonding and anti-bonding orbitals is not 
simply a difference between an unpromoted and a promoted condition if 
promotion is defined according to what happens in the united-atom (R=0). 

If, however, we redefine promotion in terms of the one-nucleus viewpoints 
of the nuclei in the molecule, introducing the new word “premotion” in order 
to avoid misunderstanding, bonding electrons are essentially unpremoted 
electrons or (from an energy standpoint) sometimes slightly premoted elec- 
trons, while anti-bonding electrons are (strongly) premoted electrons. A few 
examples will make clear the exact sense in which premotion is here defined. 
According to the one-nucleus viewpoint of either of the two N nuclei in 
N.* (or No), the seven (or eight) electrons occupying the bonding molecular 
orbitals ¢2s, ¢2p, and 72p are reckoned as 2s, 2po, and 2pz electrons, hence 
all unpremoted; in Ne they form, from the point of view of each nucleus, a 
complete, although completely shared, L shell. The anti-bonding o*2s molec- 
ular electrons, however, must be reckoned as 3-quantum, hence premoted, 
electrons. In HCI, two bonding electrons may be considered as shared by the 
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H and Cl nuclei. From the H nucleus viewpoint, these represent 1s*, from 
the viewpoint of the Cl nucleus they are 3po*. But from both viewpoints 
they are unpremoted in the sense that their quantum numbers are the same 
as those of electrons already present in the H and Cl atoms. In PtCl«> (cf. I), 
two bonding electrons may be considered as shared between each Cl and the 
Pt nucleus. From the point of view of each Cl nucleus, the electrons it shares 
are 3pla|?. From the point of view of the Pt nucleus, the twelve electrons 
it shares are 5d, ‘6s* 6p", while 5d,° is also present but unshared. The shared 
5d and 6s electrons are unpremoted, the 6/ electrons are premoted in the 
sense that the Pt atom alone, in its normal state, has a configuration 5d°6s 
without 6p electrons. But the fact that 6p electrons are only a little less 
easily ionized than 5d or 6s in the Pt atom makes it possible for the Pt 
nucleus in PtCly= to use them as bonding electrons of a slightly premoted 
type. 

The results of the preceding paragraphs can be formulated in part as 
follows. Chemical combination of the homopolar type is a result of the shrinkage 
and consequent energy-decrease of atomic orbitals in the fields of neighboring 
nuclei, when such orbitals are shared with little or no premotion. 

A condition which is necessary for successful sharing is that the binding 
energy shall not be too different in the atomic orbitals involved. Otherwise 
the sharing becomes a one-sided affair: either an electron is almost completely 
transferred from one nucleus to another, or else almost no sharing or transfer 
at all takes place (cf. Hund?® for further discussion and details). Polar mole- 
cules like Cs*F~ approximate the former and loosely bound molecules like 
Cu(NH;)4** the latter case; unstable molecules like HeH are extreme ex- 
amples of the latter case. 


9. Valence saturation 


The important phenomenon of valence saturation can now be somewhat 
understood. It should first be pointed out that saturation is usually a re- 
lative rather than an absolute matter. Whether an atom or molecule acts 
in a saturated manner usually depends on the other atoms or molecules 
with which it is placed in contact, and furthermore it depends very much on 
circumstances such as temperature, pressure, and presence or absence of 
light. In practise, however, a molecule is generally considered saturated if it 
is stable at room temperature in the presence of others of its own kind and 
of air and diffused sunlight and perhaps of certain common chemicals. With 
this definition, a given atom can occur in many different saturated mole- 
cules. 


A state of saturation usually exists when each nucleus in a molecule 
is surrounded by a set of closed shells which is stable with respect both to 
gain and to loss of electrons, including gain by sharing. A set of closed shells 
is stable toward loss of electrons if it has a high ionization potential; toward 
gain of electrons if a fairly large energy would be required to transfer an 
electron to an orbital outside the set. Closed shells ending with ns*np*, also 
the shell 1s°, usually fulfill both criteria, while closed shells ending with nd‘ 
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fulfill neither since the ionization potential is relatively low and the energy 
of excitation required to get md® (n+1)s is always small, or even negative. 
Even “stable” closed shells do, however, often take part in further sharing, 
usually by allowing two of their electrons to be shared by another nucleus 
(cf. section 4, item 2) and the examples of NH; combinations cited in section 
7). Closed shells ending with ms’, n>1, are relatively stable toward gain of 
electrons, but not toward loss. Atoms which are surrounded by stable closed 
shells when alone do not form molecules, because their electrons are too 
strongly held for either transfer or sharing and at the same time their outer 
shells contain no vacant places of low energy for receiving transferred or 
shared electrons. From the preceding discussion it will be seen that the 
property of saturation is not a property of closed shells as such, but depends 
on the occurrence of marked energy discontinuities, which often but not al- 
ways are associated with closed shells. 

Electrons are not accepted beyond the point of saturation simply because 
the unfavorable energy effect corresponding to premotion to a vacant position 
is greater than the normal favorable energy effect corresponding to sharing 
without premotion. The quantum-mechanical forces which cause saturated 
molecules to repel other molecules if they come too close may be described as 
premotion forces. The action of the premotion forces in polyatomic molecules 
is of the same nature as the anti-bonding action °° of anti-bonding electrons 
in diatomic molecules. 

In diatomic molecules there is sometimes a forced sharing of anti-bonding 
electrons, incidental to the sharing of bonding electrons, which may be con- 
sidered as a sort of supersaturation. That is, electrons are included in the 
molecule in promoted orbits which would not be expected from the valence 
rule. An example is the Ne molecule with the electron-configuration 1s? 1s? 
o2s* o*2s* r2p' o2p*. From the point of view of each nucleus, that nucleus is 
surrounded by 1s? unshared, then by 2s? 2pm‘ 2po* (shared, cf. section 8), 
forming a complete K and L shell, and in addition by two (partially) pre- 
moted electrons probably classifiable as 3-quantum electrons, corresponding 
to the shared o*2s*. The molecule NO is even more supersaturated, by the 
inclusion of a m*2p electron in its configuration.® This functions very de- 
cidedly as a premoted 3-quantum electron from the point of view of one or 
both nuclei. 

More commonly one finds an opposite condition in which nuclei in mole- 
cules are surrounded by less than a complete stable shell. This often results 
from a competition of different nuclei for electrons, or an inability of certain 
nuclei to share enough electrons to complete the shells of all the other nuclei. 
Probable or possible examples of nuclei with incomplete outer shells are B 
in BCls, C in CO3>, N in NO37, Sin SO; (cf. I). 


COMPARISON OF PROPERTIES OF MOLECULAR ORBITALS 
AND ELECTRON-PAIR BONDS 


10. Properties of molecular orbitals 
It has been stated in previous sections that the concept of molecular 
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orbitals possesses all the important properties which make the concept of 
electron-pair bonds valuable in explaining and predicting molecular struc- 
tures, and further, is less specialized and so adapted to explain a wider variety 
of phenomena. The characteristics of the two concepts can be well under- 
stood by developing them for the simple molecules H+ and He, then general- 
izing.—The reader should also refer to Hund’s valuable discussion®* of molec- 
ular orbitals (Hund’s concept c in reference 5) and their relation to electron- 
pair bonds (Hund’s concept 0). 

Pauling™ was the first to apply the Heitler-London method, developed 
originally for Hz and there yielding among other things the Heitler-London 
electron-pair bond, to the one-electron molecule H.+. Applied to He, it 
helps one to understand molecular orbitals. One finds that when H(1s) and 
H* approach, either attraction with formation of normal He*, or repulsion, 
may occur. The two modes of interaction correspond to the two energetically 
distinct states one would have for H(1s) +H? if the two H nuclei were slightly 
different in charge. 

Let the two H nuclei be designated A and B and the corresponding 1s 
atomic orbitals @, and @z. The molecular orbitals resulting when H and 
H* are brought together can be expressed in zeroth approximation as follows: 


go = cola + op), gd; = loa —_ or), (1) 


where the normalizing factors c depend on the distance R between the nuclei. 
The energy changes AE are in first approximation as follows :* 


go:AE = e/R + [e’/(1+S)\(—h—- kh) (2) 
@:AE=e/R + [e'/(1—S)|(-h+h). (3) 


Here 


= J lon: ra|dr, 5 = [ esdndr, Ip = J lon raldr, (4) 


where r4 is the distance of the electron from nucleus A. The “resonance” 
term + /2 turns the scales in favor of attraction in the case of @o, of repulsion 
in that of ¢). 

The orbital ¢o, whose presence leads to formation of a stable molecule, 
belongs to the class of bonding orbitals and may in fact be taken as the 
prototype of these, at least for diatomic molecules with equal nuclei. Simi- 
larly @; may be taken as a prototype for diatomic anti-bonding orbitals. It 
will be noted that ¢o is symmetrical in the nuclei while ¢; is antisymmetrical, 
and that ¢o is nodeless while ¢; has a nodal plane half-way between the 
nuclei. 

The behavior of the AE’s for do and ¢, (Eqs. 2, 3) is easily understood 
when one considers the mean charge density e*¢* corresponding to each. 
One finds 


go? = Co°(ba* + bn? + 2hadn), 2 = 7 (ba? + GB? — 2b4gGz). (5) 


3 L. Pauling, Chem. Rev. 5, 173 (1928). 
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This shows that e*@¢? is relatively more concentrated in the region between 
the nuclei than if one had a mere overlapping of the densities e’,? and 
e*,*, which would give 3e?(¢17+@27) for e*d?. On the other hand ¢; is much 
less concentrated between the nuclei than for a mere overlapping of e*,° 
and e*¢,”. With mere overlapping of densities, one would have a slight re- 
pulsion between H and H?*, corresponding to e®/R—e®J; in Eqs. (2) and (3). 
It is the added concentration between the nuclei or withdrawal from this 
region, in @y and ¢;, which expresses itself in + J: in Eqs. (2) and (3) and is 
decisive in making ¢» bonding and ¢; anti-bonding. 

The above results for molecular orbitals in He* can be generalized to 
other molecules. N2*, obtained by the union of N*(1s* 2s? 2p?) + N(1s? 2s? 
2p*), will serve as an example. Strictly, we should use also a certain admix- 
ture of N**++N-, and should consider various excited states of N*+N, but 
it is unlikely that this would seriously affect the results. Also, it should be 
noted that in this discussion we are treating the group of three states of N* 
(@P,1D, and 'S) and of N(4S,2D,?P), which arise from the electron configura- 
tions mentioned, as if each group were a single state. This is, however, 
legitimate because of the large energy of formation of N2*.Noe* in its normal 
state contains one o2p electron, which may be considered, approximately 
at least, as going over on dissociation to a 2p atomic orbital of the neutral 
N atom. This may be on either nucleus A or nucleus B, so one gets approxi- 
mately 


-2p = C(bA2p + ¢ B2p) ; 


where of course it is understood that the 2p atomic orbitals are based on 
a suitable atomic Hartree field. Instead of the bonding orbital ¢2p, one 
could also get the anti-bonding orbital ¢*2, which should occur in some 
excited, perhaps unstable, state of No*:¢,+2,=c*(@12»—On2,). The bonding 
orbitals ¢2s and 72p in N2* could similarly be shown to be capable of being 
respectively approximated by sums ¢;1+@, made up of atomic 2s or 2p 
orbitals, and the anti-bonding orbital o*2s could be approximated by a 
difference d,4 —@pz. 

Also in unsymmetrical diatomic molecules, molecular orbitals can be ap- 
proximated by linear combinations of atomic orbitals: bonding type, ¢=ad@, 
+b@». In strongly ionic molecules, a>6 or vice versa. The results stated in 
this and the preceding paragraph have already been given by Hund®* and 
need not be further discussed here. 

It seems clear that molecular orbitals in polyatomic molecules also can 
be approximated as linear combinations of atomic orbitals (cf. I for ex- 
amples), and that bonding orbitals are additive linear combinations. In 
general, polyatomic bonding orbitals may be expected to surround several 
nuclei, although in special cases such orbitals may be largely confined to the 
neighborhood of two nuclei. 


11. Maximum overlapping as a criterion of bonding power 


It is now of interest to learn how the bonding power of a molecular orbital 
depends on its form. This can be seen from expressions like that for @» in 
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Eq. (1). Since the lowered energy of bonding orbitals results from the fact 
that they are relatively concentrated between the nuclei, one readily sees 
that in order to have the bonding energy as large as possible, products such 
as dagen in Eq. (4) should be as large as possible, which is true if the atomic 
orbitals $4 and ¢z overlap as much as possible. This can also be seen directly 
from the energy equations (cf. Eqs. 2, 4). This is Pauling’s and Slater's cri- 
terion of maximum overlapping of atomic orbitals. The process of applying this 
criterion to the problem of determining the zeroth approximation atomic 
orbitals which give strongest bonding (Pauling’s “best bonding eigenfunc- 
tions”), when @,4 and (or) ¢z belong originally to degenerate orbitals, seems 
to be the essence of the method used by Slater and Pauling (cf. reference 20, 
p. 1141). 

Slater’s and Pauling’s original derivations of this criterion from the equa- 
tions for electron-pair bonds (cf. Eq. (8) and related AE equations) give the 
impression that the existence of the criterion depends on the overlapping of 
the orbits of two actual electrons. The present approach shows, however, that 
the criterion for finding a best bonding molecular orbital for a pair of atoms 
depends only on the overlapping of two atomic orbitals, regardless of whether 
the resulting molecular orbital is occupied by two electrons, by one elec- 
tron, or by none; the validity of the procedure, used in the present series, of 
regarding such an orbital occupied by two electrons as essentially equivalent 
to an electron-pair bond will be considered shortly. 

Furthermore, the present approach shows that the criterion of maximum 
overlapping is just as applicable to molecular orbitals which are approximate 
linear combinations of atomic orbitals of several atoms as to those which are 
formed from orbitals of just two atoms. The present use of molecular orbitals 
which connect several nuclei is perhaps the most essential difference between 
the present method and that of Slater and Pauling. 


12. Doubly- or multiply-occupied molecular orbitals 


In nearly all chemically stable molecules, every bonding molecular orbital 
which is occupied at all is occupied by the full quota of electrons allowed by 
the Pauli principle, namely two if the orbital is non-degenerate, 2m if it is 
degenerate and n is the degree of degeneracy. The most important properties 
of such multiply-occupied bonding orbitals can easily be generalized from 
a consideration of the doubly-occupied o1s orbital in normal He. Omitting 
constant factors and spins, the wave function x for normal Hz may be ap- 
proximated by the following expression : 


X = $0(1)¢0(2) (6) 


Eq. (6) implies that the two electrons move entirely independently, in 
the o1s orbital, here denoted by ¢o; the numbers 1 and 2 refer to the coérdin- 
ates of the two electrons. This, however, is true only as a rough approxima- 
tion. Although the agreement of Eq. (6) with the truth can be made rather 
good by making the best possible choice of the form of ¢9,—the simple form 
given in Eq. (1) can be greatly improved on by adjusting it to correspond 
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to a one-quantum orbit in a suitable Hartree field,—it can never be made 
exactly right without adding a correction term. This is necessary in order to 
take care of detailed effects of the two electrons on each other’s motions. The 
same kind of problem is met with in atomic spectra, for example in the He 
atom. Eq. (5) is a fairly good approximation for this atom if @o is taken to 
refer to the 1s atomic orbital in a suitable Hartree field. 

Slater* has found, in a paper preceding the work already referred to 
above, that if @o is approximated by ¢co(@4+¢z) in accordance with Eq. (1), 
the corrected expression for x can be closely approximated by: 


X = adgo(1)bo(2) — bgr(1)¢r(2), (7) 


where ¢; is as given in Eq. (1). For the molecule in its equilibrium condition 
(R=R,.), Slater finds a/b =8, approximately; but for a molecule in the process 
of dissociation (R->~ ), a/b approaches 1. This indicates that Eq. (6), even 
with @» taken as co(6,4 +z), and of course all the more so if a better @»o is 
used, is a good approximation for He in its ordinary stable state. 

Reflection indicates that it is safe to generalize the type of approximation 
contained in Eq. (6) to unsymmetrical and to polyatomic molecules (cf. 
also reference 5). It even appears that the approximation may often be rela- 
tively better for polyatomic molecules than for He. Hence the general use of 
molecular orbitals in describing shared electrons in stable molecules seems to 
be justified as a good approximation. In the case of chemical bonds with large 
R and small dissociation energy, however, cases may occur where the use of 
molecular orbitals as in Eq. (6), although still formally possible, does not 
constitute a good approximation. 

The nature of the correction required in Eq. (6) is somewhat different 
for the case of He than for that of He, because of the non-centralness of 
the Hartree field in He. The nature of this difference can be seen if we con- 
sider what would happen if we increased R to a large value in He. One would 
then find that the molecule would usually be very much like two H atoms 
with one electron on each, which means that x would be essentially of the 
form $4(1)¢2(2) or O4(2)62(1). There would, however, be a small possibility 
of finding both electrons attached to one H nucleus, in which case x would 
be essentially $4(1)¢4(2) or d2(1)¢2(2), corresponding to H~+H* or H* 
+H-. Returning to R=R,, it would seem reasonable that the true x should 
be approximately a linear combination of the four expressions just given. 


13. Comparison with Heitler-London electron-pair bond 


Heitler and London’s original approximation™ for normal He, omitting 
constant factors, was 


x = ¢a(1)¢a(2) + a(2)oa(1). (8) 


* J.C. Slater, Phys. Rev. 35, 514-5 (1930). 
2 W. Heitler and F. London, Zeits. f. Physik 44, 455 (1927); Y. Sugiura, Zeits. f. Physik 
45, 484 (1927). 
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This expresses a state of affairs in which if one electron is in @4, the other 
is necessarily in @y. For large R, this would keep the electrons separated, 
with one near each nucleus, but for R= R,, where @, and @» overlap strongly, 
the effect is much less pronounced. 

Slater®' has found that the true state of He for R= Ry is better expressed by 


x = aloa(1)on(2) + o4(2)o2(1)] + Bloa(1)o4(2) + on(Doa(2)], () 


with @/8 roughly equal to 8. [That the ratio a/8 here is nearly the same as 
a/b in Eq. (7) is accidental.] Eq. (9) allows the two electrons to be close 
together near the same nucleus oftener than Eq. (8) would permit. 

Eq. (9) is a step from Eq. (8) toward Eq. (6); Eq. (6) implies that the 
probability of any stated position for one electron is independent of the 
position of the other. On the other hand, Eq. (7) isa step from Eq. (6) toward 
Eq. (8), if dy and ¢; are as given in Eq. (1). In fact, Eqs. (7) and (9) become 
identical if Eq. (1) is used and if the coefficients a, b,a, 8, are properly related. 
Slater (private communication concerning details of Ref. 25) found that for 
R=R, in He, the ratio of the coefficients co and c,; in Eq. (1) is about that 
given by (¢:/¢o)?=6.21. Substituting this, also a/b=8, in Eq. (7), the latter 
reduces to the form of Eq. (9) with a 8=14.21/1.79, which happens to be 
very nearly the same as a/b. 

Eq. (8) is the expression, in terms of wave functions, of the Heitler-Lon- 
don, Pauling-Slater concept of the electron-pair bond, for the case of He. 
Eq. (6) is the expression in terms of wave functions, for He, of the concept 
of a bonding orbital occupied by two electrons. Slater’s Eqs. (7) and (9), in 
view of his result that a/b and a@/8 are nearly equal, suggest that the two 
concepts, while departing from the truth in opposite ways, are about equally 
good approximations to it in their ability to describe the positions of two 
equivalent bonding electrons relative to each other and to two nuclei. Other 
considerations, pro and con, may enter besides the values of the ratios a/b 
and @/8, and of course such ratios would be different for other examples than 
He, but it seems fairly safe to assume that the molecular orbital concept 
ordinarily does not give too bad an approximation in its description of the 
electron positions. We have already seen that in two other important proper- 
ties,—the fact of concentration between the nuclei and the existence of the 
criterion of maximum overlapping,—the concept of bonding molecular 
orbitals gives the same results as that of electron-pair bonds. 

For electron-pair bonds, the first of these properties follows if one de- 
termines the electron density for a single electron; this is given by /x*d72 
and comes out proportional to ’ 


ga> + on* + 2Soads, (10) 


where S is as in Eq. (4). The concentration term here is less by a factor S 
than in Eq. (5). When one takes /x°dr2 using Eq. (6), however, one gets the 
same result as is expressed in Eq. (5). The existence of the criterion of maxi- 
mum overlapping follows in the case of the electron-pair bond by a considera- 
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tion of the above electron density expression, or of AE expressions,” in much 
the same way that it does for bonding orbitals. 

It is instructive to compare the application of the method of H.L.P.S. 
electron-pairs and that of molecular orbitals to molecules which have un- 
equal nuclei and are partly polar, such as LiH, HF or HCI. If one uses elec- 
tron-pair bonds, one forms a linear combination of polar and electron-pair 
wave functions :'> 


x = alos(1)br(2) + $1(2)xn(1) | + Bdr(lidpz(2), (11) 


where B is the more negative atom, Using molecular orbitals, Eq. (6) applies, 
with 
Po = ag, + bop, b> a. (12) 


Using Eq. (12), Eq. (6) becomes identical with Eq. (11) if we put ab=a, 
6? =3, except that Eq. (6) gives an extra term a°@,(1)@,(2). This term is re- 
latively unimportant if the molecule is strongly polar: for example if b= 2a, 
the coefficient of @4(1)¢4(2) is only one fourth that of @x(1)dx(2). As in the 
case of equal nuclei discussed above, however, the truth must lie between Eq. 
(11) and Eqs. (6), (12). If both forms represent about equally good approxi- 
mations, as seems likely, Eq. (6) may well be preferred for many purposes be- 
cause of its simpler conception and formulation. 
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Evidence of Space Quantization of Atoms upon Impact 


By Hi. Keun anpd O, OLDENBERG 
Zweites Physikalisches Institut, Goettingen, and Physical Research Laboratory, 
Harvard University 
(Received May 16, 1932) 

The present paper otters an interpretation of certain diffuse bands observed in the 
spectra of metal vapors mixed with rare gases. Two separate maxima on the short 
wave-length side of the resonance line of the metal are not due to vibrational quanta 
as in the case of other similar diffuse spectra. Instead they are interpreted as separate 
electronic levels of the pair of colliding atoms described by space quantization taking 
place in the mutual approach. The diffuse structure on the long wave-length side is 
due to vibrational quanta of molecules. 


I. EXPERIMENTAL FACTS 


SEVERAL vears ago, some fluorescence experiments were performed,' the 
aim of which was to investigate the interaction in a single elementary 
process of the three kinds of energy: light, excitation energy, and kinetic en- 
ergy of atoms. It was expected that an excited atom, when colliding with a 
rare gas atom, might react in two ways: either it might split wp its energy of 
excitation into two parts, radiation and kinetic energy of the two atoms, thus 
giving rise to a continuous spectrum spreading from the resonance line to the 
long wave-length side; or instead the kinetic energy of impact might join 
with the excitation energy thus producing radiation that belongs to a con- 
tinuum on the short wave-length side of the resonance line.” 


He+A 





Hy 2576 





Mercury rare gas bands. 


In order to check this type of interaction by experiment, the fluorescence 
radiation of mercury vapor with a large amount of a rare gas added has been 
investigated. In this experiment, the excited mercury atom has a good chance 

1 (. Oldenberg, Zeits. f. Physik 47, 184 (1928); 51, 605 (1928); 55, 1 (1929), 


2 A theoretical discussion of this process has recently been given by V. Weisskopf, Zeits. f. 
Physik 75, 287 (1932) and H. Margenau, Phys. Rev. 40, 387 (1932). 
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EVIDENCE OF SPACE QUANTIZATION 


to collide during its lifetime with a rare gas atom. Actually, continuous spec- 
tra have been observed on both sides of the resonance line of mercury. Under 
most conditions, however, instead of having a uniform distribution of intens- 
ity, they showed an unexpected diffuse structure as seen in Fig. 1. (The in- 
tense Hg line 2537 itself was taken out of the spectrum by an absorption 
tube with Hg vapor, in order to prevent scattering of light at the photographic 
plate.) When A or Kr is added to Hg vapor, two broad diffuse maxima appear 
at the short wave-length side of the resonance line and a finer diffuse struc- 
ture on the long wave-length side. Only a trace of structure has been ob- 
served in Hg+Ne and no structure in Hg+He or Hg+ Xe. 

The same type of structure has recently been observed by Krefft and 
Rompe* in mixtures of thallium and other metalic vapors with rare gases. 
It is of interest that also with He added they obtained two well-defined sepa- 
rate maxima, even broader than with the other rare gases. 


II]. CRITICISM OF THE INTERPRETATION AS A BAND SPECTRUM 


The fine structure on the long wave-length side has been explained by the 
formation of molecules from mercury and rare gas atoms. In the fluorescence 
experiments only the mercury, not the rare gases, was excited. It must be 
assumed that the excited mercury atom attracts the rare gas atom by the 
polarization force (cf. section IV). The conspicuous broad maxima on the short 
wave-length side, however, which present quite a different aspect, have not 
been explained satisfactorily as yet. In order to interpret them in the usual 
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Fig. 2. Hypothetical band structure. AC, molecular band; FG, atomic line. 


way as vibrational quanta of a band, potential curves must be assumed as 
in Fig. 2. The short wave-length position of the maxima must be represented 
by a longer vertical distance AC than FG representing the line emission. It 
must be assumed that excited molecules with strong vibration combine with 
normal molecules with no vibration or with just one quantum. It is not evi. 
dent, however, why the vibration of the excited molecules is restricted to such 
a small range of high energy near A that the vibrational structure of the 
normal state comes out clearly. 

It has been suggested* that the diagram of Fig. 2 does not apply because 
it fails to represent the rotational energy of the molecule, which can exceed 


3-H. Krefft and R. Rompe, Zeits. f. Physik 73, 681 (1932). 
* Oldenberg, Zeits. f. Physik 55, 10 (1929). 
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considerably the dissociation energy. This argument cannot be maintained 
because of the selection rule. In the emission process the angular momentum 
of the molecule changes only by the fundamental unit 4/27. For heavy mole- 
cules as HgA or HgkKr this corresponds to such a small change of rotational 
energy that the single band cannot spread far from its zero line. Hence the 
rotational energy will not affect the position of the bands to any considerable 
extent. 

These criticisms based on a rather vague argument are strongly supported 
by the new observations of Krefft and Rompe, who observed the same type 
of structure with a particularly broad extension in the mixture of Tl vapor 
and He. The frequency difference between the two maxima is about 500 
cm™!, that is of the order of magnitude of the vibrational quantum of a regu- 
lar molecule.* Now it cannot be assumed that the molecule TlHe exists with 
a considerable energy of dissociation (to be observed as a large quantum of 
vibration) and that it plays a part in the discharge in 800°C, for the He atom 
has by far the smallest polarizability® of all rare gases (by far the lowest con- 
densation point, the smallest heat of vaporization, the least adsorption). 
Therefore the type of structure as observed in Tl+He on the short wave- 
length side of the Tl line must be interpreted by an entirely different process. 


II]. RADIATION DURING COLLISION; SPACE QUANTIZATION 


It seems that this process can be understood on the basis of the theory of 
band spectra without an arbitrary hypothesis. It must be assumed that the 
excess energy of the short wave-length diffuse structure, not explained as 
dissociation energy of molecules, is due to the kinetic energy of thermal im- 
pacts.’ The mercury atom emits this radiation instead of the sharp resonance 
line at the moment of a collision with a rare gas atom. This is just the process 
mentioned above, for the study of which the experimental work was under- 
taken. This process has been described on the basis of potential curves by 
Winans and Jablonski.’ By the same mechanism Weizel® recently explained 
the asymmetric short wave-length broadening of some helium lines discovered 
by Hopfield.!° Weizel has succeeded in connecting this broadening with the 


® The correction discussed by H. Kuhn (Zeits. f. Physik 63, 458 (1930) and 72, 462 (1931)) 
cannot change the order of magnitude of the observed difference, while still discrete maxima 
of intensity are observed. 

6 J. C. Slater and J. G. Kirkwood, Phys. Rev. 37, 682 (1931); H. Margenau, Phys. Rev. 
37, 1425 (1931). 

7 From this assumption it follows that Condon’s wave-mechanical interpretation of diffuse 
maxima (Phys. Rev. 32, 866 (1928)) does not apply to the spectrum under discussion. Condon 
considered a sharply defined initial state (lowest vibrational level of the upper electronic state) 
and a final state which belongs to a continuous range. He concluded that diffuse maxima might 
appear closely connected with the waves representing the motion in the final state. This theory 
does not apply to the conditions described by Fig. 2, because here the initial state, too, is not 
sharply defined but is itself diffuse due to the velocity distribution. Moreover the obseved 
maxima seem to be broader than those described by Condon’s theory. 

8 T. G. Winans, Phil. Mag. 7, 558 (1929). A. Jablonski, Zeits. f. Physik 70, 723 (1931). 
® W. Weizel, Phys. Rev. 38, 642 (1931). 
10 J. J. Hopfield, Astrophys. J. 72, 133 (1930). 
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molecular spectrum. Discussing the broadening of the He line 585 on the 
basis of the theory of band spectra, he finds one potential curve that de- 
scribes the interaction of two normal He atoms colliding with each other, 
and four separate curves for an excited 2'P atom colliding with a normal one. 
Due to a selection rule only two of these four curves, both representing re- 
pulsion, combine with the normal state (Fig. 3). 

It can be assumed that in the normal state the Hg atom has a smaller 
effective cross section than in the excited state. Hence in the diagram the 
potential curve of the normal state remains flat even for a small internuclear 
distance. Thus the kinetic energy of the collision—at the very moment in 
which it is turned into potential energy (A in Fig. 3)—takes part in radiating 
a larger quantum AC than the excited atom would radiate by itself (FG). The 
superimposed effect of the two potential curves of the upper state is observed 
as pressure broadening of the atomic line. 

In the present problem—mercury with an added rare gas—the interaction 
of the two atoms is simpler because they are different from one another. In 
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Fig. 3. Radiation during collision. AC, radiation during collision; FG, atomic line. 

Fig. 4. Potential curves for the mercury and rare gas collision. AC, BD, diffuse maxima 
on the short wave-length side; AQ, BP, limits of the spectrum; KL, MN band structure on 
the long wave-length side; FG, atomic line. 


the normal state of the two atoms again only one potential curve exists (repre- 
senting an ' term), because both the mercury and the rare gas atom have 
spherical symmetry in their normal 'S states. In the upper state (excited 
atom-+normal rare gas atom) the method of counting the molecular elec- 
tronic levels, given by Hund," depends upon the coupling. In Weizel’s prob- 
lem (He’+He) the multiplet structure (due to the interaction between the 
orbital momentum and the spin) is narrow as compared with the extension 
of the continuous spectrum (due to the interaction between the colliding 
atoms); this interaction between the atoms is described as a quantized orien- 
tation of the orbital momentum with regard to the internuclear axis. In the 
present problem, however, it must be taken into account that the splitting 
up of the multiplet structure in the mercury or thallium atom is much larger 
than the extension of the diffuse molecular spectrum under consideration. 
Hence the interaction between the orbital momentum and the spin is large 


" F, Hund, Zeits. f. Physik 36, 659 (1926). 
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as compared with the interaction between the two colliding atoms and a 
quantized orientation of the total angular momentum of the Hg or Tl atom 
is to be expected with regard to the internuclear axis, just as in the case of a 
Stark effect in a weak electric field. The number of separate energy levels of 
the Hg’ or TI’ colliding with a normal rare gas atom is given by the number of 
terms in the series j, j—1.--.-, 0 or 3. (Negative values of j describe an- 
other orientation but give the same energy levels as the positive values of j).” 
This argument leads to the conclusion that in the upper state Hg’(2°P,) +rare 
gas (1S) two molecular levels exist, given by j/=1 and j7=0, represented in the 
diagram by two separate potential curves. From both excited levels the 
transitions to the normal level are allowed. (As in general the collision will 
not be head on, a rotation will be superimposed with a constant angular 
momentum. ) 

Along one or the other of these curves the rare gas atoms and the mercury 
atoms will approach each other to a certain range of potential energy A or B 
determined by the thermal energy. 

The effective cross section of the Hg atoms in the normal state will be 
supposed to be smaller than in the excited state. Hence the lowest potential 
curve of Fig. 4 representing the normal state will just barely begin to rise for 
the values of nuclear distance that belong to A and B. Therefore both transi- 
tions AC and BD, that are vertical in the diagram—following the Franck- 
Condon rule—represent radiations of larger energy or shorter wave-length 
than the atomic line FG, BD shifted by a larger amount than AC. This rea- 
soning leads to two separate diffuse maxima on the short wave-length side 
of the atomic line, both of them shifted by an amount smaller than the mean 
kinetic energy of the collisions. This just describes the observed spectra. The 
two separate maxima are interpreted not by the vibration of a molecule but 
by two separate electronic levels of a pair of atoms to be described by two 
separate directions of mutual approach, both of them belonging to the same 
electronic levels of the atoms. Thus the separate maxima give evidence of 
space quantization in the collision between atoms. This interpretation is 
strongly supported by the fact that the two separate maxima persist at high 
temperature (cf. section IV). 

The diffuse bands show a rather sharp limit on the short wave-length side 
and in some cases a much longer diffuse extension to long wave-lengths. 
This may be understood, too, on the basis of the potential curve diagram 
(Fig. 4). The dotted lines represent transitions in the case of a small change 
of internuclear distance during radiation. The long distance BP, representing 
radiation of short wave-length, cannot be much longer than BD, representing 
the most probable transition from B. AQ, however, representing the long 
wave-length extension, can be considerably shorter than AC due to the steep 
branch of the lower potential curve.” 

® Cf. W. Pauli, Handbuch d. Physik XXIII, p. 246. 

‘8 The extension to long waves plays the major part in the absorption experiments. at 
high pressures of Fiichtbauer, Joos and Dinckelacker (Ann. d. Physik 71, 204, (1923)) and 
Bonhoeffer’s and Reichardt’s experiments with Hg dissolved in liquids (Zeits. f. Physik 67, 
780 (1931)). 
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The appearance of the same diffuse structure in the absorption spectrum 
will be explained by the inverse process, the absorption being produced by a 
mercury atom at the moment of a collision with a rare gas atom. 

That these separate maxima were not observed by Hopfield in the spec- 
trum of helium, was probably due to their position in the extreme ultraviolet, 
where the resolution of the same differences in frequency is much more diffi- 
cult. 

The diffuse maxima, reported by Krefft and Rompe in thallium vapor 
with an added rare gas near the green thallium line 2*.S,,.—2?P3,2., are ex- 
plained by separate molecular levels (j= +3/2 and +1/2) belonging to the 
lower atomic level instead of the upper one. All other observations, too, can 
be interpreted on the same basis. One cannot expect, however, to observe 
separate maxima in every case, for which they are derived from the theory, 
because they may overlap too much. This is evidently the case in Hg+He, 
where the existence of two maxima can only be deduced by analogy with the 
other rare gases. 


IV. QUANTIZED MOLECULES FORMED OF EXCITED 
MERCURY AND RARE GAs ATOMS 


This explanation does not take into account the finer diffuse structure on 
the long wave-length side of the atomic line (Fig. 1). Its previous interpreta- 
tion as oscillation quanta of a molecule will be maintained.“ The simple 
structure and its position on the long wave-length side is most easily ex- 
plained—contrary to the explanation published previously—by the forma- 
tion of molecules of Hg atoms in their excited state with rare gas atoms. One 
expects the excited Hg atom to be able to induce an electric moment in the 
rare gas atom and thus by polarization form a loosely bound molecule. The 
normal Hg will have this power to a smaller extent because it has a high sym- 
metry similar to that of a rare gas atom. Therefore the upper potential curves 
of Fig. 4 ought to show shallow minima for large values of the internuclear 
distance, for which the lower curve is approximately flat. This type of excited 
molecule would emit a diffuse band on the long wave-length side of the atomic 
line, in which the vibrational quanta still might be distinguished. 

Hence this band system is explained by a group of vibrational levels in the 
upper electronic state combining with a lower state in which vibrational levels 
can not be resolved. It has been described as an apparently simple set of 
bands with a trace of convergence to the long wave-length side; it has been 
mentioned, however, that possibly the observed set of bands is due to the 
superposition of several sets. According to the present interpretation, two 
simple sets of bands (KL and WN, Fig. 4) belonging to the two upper poten- 
tial curves should be superimposed. The smaller slope of the potential curve 
of the normal state might modify the position of the observed bands. There- 
fore no conclusion can be drawn concerning the convergence from the struc- 
ture as shown in Fig. 1. 





4 Cf. F, London, Zeits. f. Physik 11, 248 (1931). 
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This interpretation of the diffuse structure on the long wave-length side 
as Vibrational quanta of a band system is supported by observations. This 
structure does not appear for the light atoms He and Ne, which—due to 
their smaller polarizability—might not form molecules. (It remains uncertain 
if this type of structure exists for Hg+X. The pictures with Xenon were 
taken with a quartz bulb of a few millimeter diameter, from which the scat- 
tered light was so strong that the structure might have been concealed.) At 
high temperature the fine structure on the long wave-length side disappears, 
while the broad maxima on the short wave-length side remain (observed in 
Hg+A). This is understood as the complete dissociation by heat; in these 
circumstances there will remain only the broad maxima on the short wave- 
length side that are explained not by vibrations of molecules but by separate 
electronic levels of a pair of colliding atoms. 

V. CONCLUSION 

We come to the conclusion that three different processes are to be distin- 
guished in the collision of Hg’ with a rare gas atom. In most cases the atoms 
will separate without any effect, leaving the Hg’ such that it radiates its 
sharp atomic line (FG in Fig. 4). This may be the most probable process. In 
some other cases the radiation will take place in the moment of the collision 
giving rise to the broad maxima on the short wave-length side of the atomic 
line (AC and BD in Fig. 4). Finally in a few cases, if the temperature is not 
too high, some energy may be taken away in a triple collision; real excited 
molecules with vibrational quanta will be produced emitting a band spectrum 
(KL and JN in Fig. 4). 

The diffuse maxima in the fluorescence radiation of mercury vapor with a 
rare gas added are surprisingly intense. The corresponding spectra in the 
neighborhood of some other lines, for instance the D lines of Na, are much 
fainter or missing entirely. This leads to the question, whether in some 
atomic terms, but not in all of them, the collision process might produce an 
instantaneous increase in the probability of emission. This problem will be 
further investigated. . 

According to the present interpretation, the main conclusion is that the 
two separate maxima of intensity, observed in the spectra of certain metal 
vapors mixed with rare gases, furnish evidence that the mutual approach of 
these atoms can be described by space quantization. 


1% M. Born and J. Franck, Ann. d. Physik 76, 225 (1925). 
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The Equation of State of a Non-ideal Einstein-Bose or 
Fermi-Dirac Gas 


By G. E. UHLENBECK AND L. GROPPER 
University of Michigan 
(Received May 25, 1932) 

With regard to the question if from isotherm measurements one can obtain an 
experimental test for the existence of Bose statistics in real gases, as is required by 
theory, we prove the following general theorem. The “Zustandsumme” of a non-ideal 
Bose or Fermi gas is given by the classical integral provided one replaces the Boltz- 
mann exp(— ¢4;, kT) factor by: 

e~*ii/*T (1 + exp | — 4x2mkTr;,2, h?| (1) 
for each pair of molecules (ij). For the second virial coefficient B, this has, i.e., in a 
Bose gas, as a consequence that: 


B = B non-ideal class + Bideal Bose + B’ (2) 
where: 


B’ =24nN f drr*(1 — eP/4T) exp [— 4x2mkTr?/h?]. (3) 
0 


Only at very low temperatures do the last two terms in (2) become appreciable. They 
are then of the same order of magnitude, but have opposite signs. Due to this fact, due 
to the lack of precise knowledge of the molecular forces, and due to the absence of 
accurate measurements of B at very low temperatures, one can as yet not decide from 
isotherm measurements alone whether or not real gases obey the Bose statistics. 


I. INTRODUCTION 
1. 
SIS well known, it has been proved from the fact that electrons and pro- 
tons obey the Fermi-Dirac (F.D.) statistics, that a gas consisting of par- 

ticles of even charge obeys the Einstein-Bose (E.B.) statistics.! A neutral 
gas, for example helium, is a special case of this. The question therefore be- 
comes very important, if one could for example verify this theoretical predic- 
tion from equation of state measurements, in particular from the experi- 
mentally determined values of the second virial coefficient B. 

The theoretical expression By, for a non-ideal gas, obeying Boltzmann sta- 
tistics is, as is well known? 


Ba = ann f drr2(1 — e~ 6) / kT) (1) 
0 


N is the number of molecules per mol. It is here assumed that the intermole- 
cular forces are radial and have the potential ¢(r). 

The expression Bz.g. or Br.p. for an ideal gas obeying Einstein-Bose or 
Fermi-Dirac statistics, is: 

1 E. Wigner, Sitzungsber. der Ungarischen Acad. 1928, p. 1; P. Ehrenfest and J. R. Oppen- 


heimer, Phys. Rev. 37, 333 (1931). 
2 See e.g., R. H. Fowler, Statistical Mechanics, Ch. VIII. 
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Ber. \ 1 h? sti 
> = F- —l[ — ) V (2) 
Brn. f 24\amkT 


In the recent work on the quantum theory of the equation of state of real 
gases (London, Slater, Kirkwood, Keyes, Margenau) one has made only use 
of (1), introducing for ¢(r) the theoretical expression obtained from a gen- 
eralized Heitler-London calculation.’ The reason that (2) was neglected is, 
that only at very low temperatures (say below 20°) it becomes appreciable, 
and because it was always assumed that the B for a non-ideal Bose gas was the 
sum of (1) and (2). 

We wish to show in this paper that this additivity is not strictly true, but 
that in addition there occurs a term due to the interaction between the Bose- 
ness and the non-ideality of the gas. We get for a non-ideal Bose gas: 


-B = Bat Bez. + Ben. (. 


we 
— 


where: 
B’r an. = anv f drr2(1 — e%!*T)exp [— 422mkTr?h?| (4) 
0 


2. The proof of (3) and (4) rests upon a kind of generalization of the 
exp (—V/k7T) theorem of Boltzmann for an E.B. or F.D. gas. Recently 
Slater? has pointed out the analogue of this fundamental theorem in the 
quantum theory. He remarked that quantum mechanically, due to the ortho- 
gonality and normalization of the wave functions yy, belonging to the energy 
values E,,, the Zustandssumme Sg may be written: 


Se = >G,e En/kT — J 7s fax = dxy doen Fn sty 2. (5) 
n n 


In Boltsmann-statistics for a given E,, all the G,, linearly independent wave 
functions y, belonging to it, are allowed, and must therefore occur in the 
sum,* so that the G, does not explicitly appear any more. In Bose statistics the 
G, =1, because only those £, are allowed, to which belong the symmetric 
wave functions, and only these have therefore to occur in the sum. Similarly 
in the F.D. statistics we have the E, to which the antisymmetric functions 
belong. Hereafter we will distinguish between these three cases by writing 
Sos: SOE.B.» Sor.p.+ 

Slater compares (5) now with the classical “Zustandsintegral” after in- 
tegration over the momenta: 


3 F. London, Zeits. f. Physik 63, 245 (1930); Zeits. f. phys. Chem. B. 11, 222 (1930). 
J. C. Slater, Phys. Rev. 32, 349 (1928); J. C. Slater and J. G. Kirkwood, Phys. Rev. 37, 682 
(1931). H. Margenau, Phys. Rev. 36, 1782 (1930); 37, 1014, 1425 (1931); 38, 747, 1785 (1931); 
Proc. Nat. Acad. 18, 56, 230 (1932). J. G. Kirkwood and F. G. Keyes, Phys. Rev. 37, 832; 38, 
576 (1931). J. G. Kirkwood, Phys. Zeits. 33, 39 (1932). One has made small corrections in (1), 
though, to allow for the existence of discrete states (formation of polarization molecules). 
‘J.C. Slater, Phys. Rev. 38, 237 (1931). 
5 Comp. P. Ehrenfest and G. E. Uhlenbeck, Zeits. £. Physik 41, 24 (1927). 
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2armkT\3* 2 ; 
Sa = ~— f te fe ‘Tdx,-- + dty (6) 
2 


where (x1 - - - gy) is the total potential energy due to the intramolecular 
forces and eventually also due to external forces. 

One might expect from the correspondence principle that for sufficient 
high temperatures where the influence of the discreteness of the translational 
energy levels can be neglected with respect to the influence of the intermolec- 
ular forces (which is already the case, excluding the electron gas, say for 


7 >1° absolute) °® 
sstea te 2rmkT\3X"2 _ 
done En Ty, 2 —— eV IAT. (7) 


(T>1°) 
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Slater makes this also plausible from the meaning of y,? as a probability 
density. In Chapter III we will give an analytical proof. This gives the 
justification even in the quantum theory for the use of (6) in the case of a 
Boltzmann gas. 

For the case of a E.B. or F.D. gas we will prove a theorem analogous to 
(7): 


»—E,/ kT 2 
doz.n.c~*!*Ty, seeieen 
: ] 1 /2amkT\3*!2 
> — — Cc V/AT x 
> F.p.e En/kTy 2 \ (T>1°).V! h? 
n - 


I](1 + exp [— 42°mkTr;2/h?] 


‘) 


(8) 


where the product has to be taken over all pairs of molecules, and r;; is the 
distance between the (7j) pair. This again holds as indicated for temperatures 
so high that the influence of the discreteness of the translational energy levels 
is negligible, but not yet so high that the deviations due to the E.B. or F.D. 
statistics from the classical gas laws can be neglected. These are given by (8) 
in first approximation; in addition therefore the volume must not be too 
small. One sees from (8) that there is an apparent attractive force in the E.B. 
gas between each pair of molecules with the potential: 


@’ = — pT log }1 + e~ det mk h2r? } (9) 


Of course for temperatures so high that the deviations due to E.B. or F.D- 
statistics can be neglected compared to the non-ideality of the gas (7, say for 
He, >20° absolute) (8) goes over into (7). 

3. Because the total potential energy V (in absence of external fields) can 
be written in the form: 


V= Yo(ris) 


(8) can by using (9) be put in the form, 
1 /2xmkT 


3N/2 
Deus. = (= ) []e~ at ou? /47, (8a) 
2 


N! h? 


6 We will distinguish the integrands of Sgz, Sgz.g by subscripts B, E.B. etc.and sometimes 
write for the sums simply Lg, Lz. g. etc. 
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Exactly as (7) is carried through by way of (6) to give (1), so (8a) is carried 
through and will give again (1) when we replace there ¢ by ¢+¢’. Introduc- 
ing then the value of ¢’ we immediately obtain (3) and (4). 

To obtain an idea of the contribution of Bg.z.’ to the second virial coeffi- 
cient, we have made in the case of helium a graphical integration, introducing 
for o(r) the theoretical expression used by Kirkwood and Keyes, and the 
slightly different one proposed by Margenau, as well as the half empirical 
one calculated by Lennard-Jones.’ We got: 


a. r= 20° ° Byrn. = — 0-8 B’ ror. = 


b. T=5°: Berw. = —6-4 Bex. S 


5 
+ + 
- 


We see therefore that Bg.z. and Bg.z.’ are of opposite sign and of the same 
order of magnitude at these low temperatures. At higher temperatures of 
course both go to zero rapidly. It is clear therefore that because of this fact, 
it will be still more difficult to decide from the experimentally determined B 
values at low temperatures, between the classical and the E.B. statistics. In 
fact, with the uncertainties still existing at present, both in the theory of the 
intermolecular potential ¢(7) and in the experimental material, in spite of the 
rather great accuracy and extent of the latter,’ it seems to us rather hopeless 
to try to make this distinction with this method. 

Quite analogous results are obtained for the F.D. statistics, the signs be- 
ing merely reversed, 


Brio. = — Ben Brn. = — Ben. 
so that also the existence of the F.D. statistics may not be ruled out. 


Il. A New TREATMENT OF THE IDEAL GAs 


4. In order to prove (7) and (8), it will be convenient to show it first in the 
case of the ideal gas, because as we will see in Chapter III, the proof for the 
general case can be made to rest upon the evaluation of the sum for the ideal 
gas. Consider first as the simplest case two atoms in a one-dimensional box of 
length L. 

a. With Boltzmann statistics we have to take ail linearly independent 
wave functions belonging to each energy value, and we get therefore for the 
sum: 


4 Nywxy Not Xe 
22 L? Do Dd erat? L2(m2+ ma?) (2)8nn, 2 sin? —— sin? — 
. ny 


4 4 





2no 


| Wyk. | Net xX, 
+ sin? ——— sin* ——— 
. 10 
& = = Nywx, Nom X9 (10) 
= — Dien art/L2(n 2+n 2) sll cadena catioiae 
L* aint ngut 4 L 


7 Comp. R. H. Fowler, Statistical Mechanics, Ch. X. 
’ Comp. G. P. Nijhoff, Dissertation Leiden, 1928. 
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) Nywx 2 P 
v4 ashes a wx, 2 sti ; Net Xe 
= E > e7 ar? Len 3 sin? we { > ener? Leng? cin? 


t <= 4 2 L 


where: 
a = h?/8x*mkT 


In the neighborhood of zero values of x; and x2 we may replace the sums by 
integrals and obtain: 


or = 


2armkT 


h? 





(1 — e7,?/*)(1 — e-7,7/*), (11) 


From the symmetry of the problem it is clear that for x; and x, in the neigh- 
borhood of L we have an analogous expression replacing x; and x2 in (11) by 
L—x, and L—x2. We see therefore that =, is zero at the boundary of the 
square OA BC and rises rapidly from all sides to the value 2rmk7/h?. In the 
center the wave functions in (10) will average out, because of their different 














4 C 
L 
O L A X 


Fig. 1. 


phases. Replacing them by the average value 1/4 and the rest of the sum by 
the integral, (10) again becomes 2mmkT/h®, so that Xg, remains constant 
throughout in accordance with (7). As we will see in §5, the exponentials in 
(11) account for the correction due to the discreteness of the energy levels. 


b. With E.B. or F.D. statistics we get: 


Zz 2 mya xX Not X2 
E.B. 2 _ MTX, | Net Xe 

\ =— )) eremstsin sea) teod sin —— sin —— 
ee.) L? nya ng L 


4 4 





We see that OC is for each term of the sum a symmetry line; in the F.D. case 
it is a zero line, in the E.B. case a maximum line. This will therefore also be 
the case for the total sum. We get working the bracket out and going over 
from the double sum to products of single sums as in @: 


Dens. 
> r.p. 


utr 


1 c) 
\ =} > 2 | Dewits| cos > — Xe) 
n=0 





2L? 
ur 2 
— cos 7 + x2) . 


-) 


(13) 
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All the terms in the first sum of the bracket are in phase along the line x, = x2, 
whereas the second sum has a value only in the corners O and C; in the rest 


of the square the second cos averages out to zero. Therefore we can neglect 
the second sum with respect to the first sum. Neglecting also the exponentials 
in Lz (this can really only be justified for N particles, see $5) and replacing 


the sum by an integral, we get: 
2rmk y. 


This is in accordance with (8). We have therefore practically everywhere half 
the classical value; only along OC it rises in the E.B. statistics sharply to the 
classical value, and in the F.D. statistics it decreases equally sharp to zero. 

5. 

Going over now to NV particles in a one-dimensional box of length L, we 
have: 

a. for Boltzmann statistics: 


2rmkT\*%'? 
Le 2 (—— IQ -— e*?*). (15) 
i=1 


h? 


This is clearly (11) extended to .V particles. From this follows the free energy 


y, according to: 
L 
core f ~ fare dey a (16) 
0 


Of course one must remember that (15) holds only for every x; between 0 and 
3L, and that for the other half we have an analogous expression replacing x; 
by L—.x;. Because of this symmetry we can write Eq. (16): 


2rmkT\*!2 etl? . 
evi kT = | ——_—— - 2% f - fan dis dxy Ila — Sere) 
h? 0 i 


Of all the terms in the product we need only those which after integration give 
terms proportional to L* and L*~'; because we are only interested in the 
second virial coefficient. We then obtain: 


guar (=) ‘yy E N ( h? ) ] (17) 
a = cnepuemenpananatenns ‘ oe sm ee / 
: h? V \SrmkT 


writing V for L now. This gives for the pressure: 


oy [14 1 1 ( h* ) | (18) 
a aves OV V 23/2? \amkT 


This agrees with the result obtained by the direct calculation of the Zustands- 
summe.’ The second term gives the correction due to the discreteness of the 





® Fora gas of N particles in a n-dimensional vessel this gives: 


_NETC 1 1 ( h? "] 
i ri 23/2 Yilm \aemkT 
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energy levels. This correction is independent of N, in contrast to the correction 
due to the E.B. or F.D. statistics which are proportional to N. For this reason 
the above correction is completely negligible. Because it arises from the ex- 
ponentials in (15) we will from now on neglect these. 

b. For E.B. or F.D. statistics we get for the sum: 


Lie.w.y a i 
= — b > =e Diem art L2(n 24+ ++ my?) g , 

" is ‘ Sipe hy 
>r - f L* N ! an 


. ny2ngz°-- 
(19) 
‘ns MTX, _ Uytxy |? 
> sin ——— - - - sin ——— 
P I 


P (-1) < < 


where the sum goes over all permutations of the x; keeping the quantum 
numbers fixed. The factor n,...my is one, when all the »; are different; is “1 
when two are the same and the other different; {; when three are the same, 
etc. When we work the square out we get first of all the square terms, which 
are just the Boltzmann terms evaluated in a. From the rest of the terms we 
must select only those which after integration over x; will give terms pro- 
portional to L*~'. One can easily convince oneself that these terms are of the 
type: 


Mgt NX, | | NGM Xe 
sin? —-— sin? ——— - - 
Us een Ue ee oe Pe Ln (20) 
ia —— da —— sin —— te {ere TT 
L ‘ ; L | sin? ——— sin? 
L L 
| 
\ 
As in $4 we write for the first four sin: 
nr ny 7] 
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cos —(x; — 2;) — cos —(x; + Xj) | X 
he F ( ce 1’ j 


~ _ 


Net Not 7 
cos —(x,; — x;) — cos —(x; + Xj) 
. 4 I 


4 “ = 








For the same reason as in §4 we may neglect again the cos of the sum x;+.7;. 
Going over from the multiple sum to products of single sums, and replacing 
the sums by integrals, we get: 


Dees.) 1 (we) + 1 h? 
+0) ~N! zi 


h* 


x} J g astiEint eos —(%; - rant | 
ij 0 L ) 


(21) 
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The sum and product are to be taken over all pairs. The last two lines differ 
only in terms which after integration over x; give terms proportional to 
powers of LZ lower than N—1, and therefore they have been neglected. 
Formula (21) agrees with (8), when we put there V =0. 

To check (21) further we may calculate the free energy. One obtains: 


ewkT = f- fas tee avy} Les. 
Drv. 
ae : | 1 V(V - —(- h? ) | 
Smt yal RY 
N! h? v 2 darmkT 


where again we write V for L. Neglecting the 1 with respect to .V, we get for 


the pressure: 
ay ml f _N1 ( h? ) 1 ” 
| — + —_— —_—_-- Z2 
, OV V V2? \armkT ' 


which agrees with well known results.'!° Comparing (15) and (21), and their 
consequences (18) and (22) resp. one sees that the reason why the E.B. or 
F.D. correction has the factor NV, in contrast to the discreteness anes is 
simply that the number of pairs, w hich one can form with V objects (x,), 
3N(N-1). 

These considerations can immediately be extended to a gas in a three- 
dimensional rectangular vessel. One obtains then instead of (21): 


Den. 1 /2rmkT\3*'? | 
VT Hd = e~ Cla) { (zi aj) 2+ (yi uj) 2+ i355 )?} (21a) 
>r.p. N! h? ‘j 


so that one sees that the apparent potential due to the Bose-ness is radial 
(see Eq. (9)). 











Ill. Tue Non-IpDEAL E.B. or F.D. Gas 


§6. 
The wave equation for the whole gas we write as: 
(H — Ey = 
where: 


H = Yo —(bet + bet + det) $V 
_ pa + py? + pzi »)+ V(m--- Zn) (23) 


Li + Dr 


where ¢;; is the interaction potential between the particles and 7; an eventual 
external potential. We shall abbreviate always all the coordinates by g and 
all the quantum numbers by 2. We wish to show in this paragraph that: 


10 For an E.B. or F.D. gas in an n-dimensional vessel one finds according to the usual 


method: 
a Nk — 1 ( hh? \n/2 “] 
VL om =a) V 
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Loe Batty? = Didae*! My. (24) 


In the right-hand side // is to be considered as an operator working on @,, 
which are the normalized eigenfunctions for the ideal gas.'' Formula (24) is 
true regardless of the statistics, that is to say, it holds equally well when we 
take on both sides only the symmetric wave functions (E.B. statistics), or 
only the antisymmetric wave functions (F.D. statistics), or all the wave 
functions (Boltzmann statistics). 

For the proof one remarks first that for any function (power series) f(x), 
we have: 


S(EWD Wn = [UD 
in which // again is the operator, in particular: 
eal kTY, = e~H/ktTy, , (25) 


We develop now y, according to the normalized eigenfunctions of the ideal 
gas. In the case of a cubicle vessel this is simply a development in a multiple 
Fourier-series 





¥n(q) _ >ocnmbm(q) (26) 
m 
where one must then remember that @¢,,(q) is an abbreviation for :” 
2\8%? myx, . Mery |, Msn TEn 
Pm,-+-mgN (21 as? ty) - (-) sin L sin L -*-Sin > . 


Substituting (25) and (26) in the left-hand side of (24) we get: 


DivneFnltty, = DED Lenicndie "7; = Didie#/*%6, 
n n a 7 i 


since: LaCniCnj = 5i;. 
§7. 
Introducing for simplification: 


V’ = 8x°mV/h? E’ = 82?mE/h? 
the wave equation becomes: 

d*y /dq? + (E' — Vw = 0 
and for the sum, using (24) we get: 


a . ¢ -_ , 2 —V’) 
Loe Fal tty,? = Dierabny,? = Didaer(eide Vg, 
n n n 


Comp. F. Bloch, Zeits. f. Physik 74, 295 (1932). We came to this result independently 
and because our point of view is somewhat different, we thought it not superfluous to repeat the 
proof. 

2 This is in the case of the Boltzmann statistics. For the E.B. or F.D. statistics one takes 
the symmetric or antisymmetric combinations, and because these form with respect to sym- 
metric resp. antisymmetric functions a complete orthogonal set, the argument remains the 
same. 
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where a has again the meaning h*, 827°mk7 and where g abbreviates again the 
3.N coordinates of the gas. 

We will prove now that when a is very small, we have in first approxima- 
tion: 


2 ino? , E - , ° ~ 
> o,,€° d2/aq?—\ >», S¢ al > en ae 2] ) 
” 


where /*e,,', 82° are the energy values of the ideal gas. In particular for a cu- 
bicle vessel 
€,) = w7(ny? + nto? + - ++ mgy?)/L?. 

In this way the evaluation of the sum for the non-ideal gas reduces to that 
for the ideal gas. Using then the results of Chapter II, we will then have 
proved our main theorems (7) and (8). 

Formula (27) is not obvious, because d*/dq? and V’ are not commutable, 
so that in general: 
ya (d2/dq2-V" = en V’+d2/ dq?) + c aV, ead? dq? 

= ert? dq? . @ av’ 
To prove (27), write: 
fF = ea (a? dq? Vd, 
then F will fulfill the differential equation :" 
OF, 0a = (a 0q") —_ VF 
which we rewrite: 
OF /da + (e,’ + V’)F = (6°/dg? + €,)F (28) 
Fora=0, F=¢, and the right member is zero. In first approximation there- 
fore: 
Fm gyrate (29) 
which gives (27). In Note I we will prove that in higher approximations: 


Lig neriae—vg, = eel" Diem atng,2- [1 + afi(q) + afe(q) +--+ | (30) 


which justifies more strictly that (27) is true for small a; fi(q), fe(q) are cer- 
tain functions of the coordinates, which contain V and its derivatives. 

As a special example, which is of interest because the summations can be 
carried out exactly, we will consider in Note II the case of the harmonic 
oscillator. 


Note I. 

To find the next approximation we substitute (29) in the right-hand side of (28). Considered 
as a differential equation for F as a function of a, we can immediately integrate. Making use of 
the Schrédinger equation: 


13 Comp. F. Bloch, reference 11. 
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(d?/dq? + €x')d. = 0 


: a§( dy’ do, dv") a /dV'\? 4 
F = [ +. = soe + 2—: oa - -6.(<-) Jew ) (31) 


we obtain: 


dq* dg dq ig J 3 dq 


Substituting in the left-hand side of (27) and calling the sum for the ideal gas: 
= > en Gn? 

we get: ; 

( ay’ dv’ S) 


pert ldeV dy, = eV" I 5 — S+— 
X i” a ( dq? dq dg 


4 “(<) 5! 

3 \ dq - 
Substituting the value of S from (15) and (21), only the third term of the right-hand member will 
give a term proportional to a. Neglecting the terms which give a? and a’, we obtain: 

a. for Boltzmann statistics: 


” 1 ” 2 IV aye t lo 
Se= (Keele - [1-0 Sat] ae 
Ta i ; ve : 


b. for E.B. or F.D. statistics: 
DYER 1 sf. 
Srp) N 


L 


We must remark, that the terms proportional to @ in (32) and (33) contribute to the second 
virial coefficient, but as one easily verifies they will give corrections to (3) which are propor- 
tional to a’ (for a gas in a three dimensional vessel), whereas (2) is only proportional to a®”. 

Further one must remark that the second term in the bracket of (32) and (33) are not the 
fi(g) of (30), for if one proceeds to the next approximation, one gets another term proportional 
to a, which contains 6?V'’/dx,;*. One can convince oneself though that in the second virial 
coefficient these terms give corrections proportional to powers of a higher than a’, 

One can approach the proof of (27) from a somewhat more general standpoint. When f and 
g are any two non-commutable quantities, one can show that for small values of a:" 


etito) = {1 + (a? 2)(gf — fg) 4... \ eS. @%9 


"Tas t essai) 


x (xi— xe Se (ee av’ -) oni J (33) 


Dte@ 7-222 \ax; ax; 


rn + 


a 
4 


When we identify now g with d?/dg? and f with V’, and let exp[a(f+g) ] work on ¢, one obtains 
again (31).% 
Note II. 
For the harmonic oscillator it is possible to carry out the sums exactly. Introducing the 
abbreviations: 
= 2rx(mv/h)'? 6 = hv/kT 


we can write: 


. He 
(2emkT)—2 F e-EnltTy,2 = (29) /%¢-8/2-v? Feo no! 0) (34) 


" 2"-n! 


4% Comp. J. E. Campbell, Theory of continuous Groups, pp. 54-57, for analogous questions. 
15 One uses here the special case m = 2 of a general formula for 


f ad” d" y, 
dg” dg" "i 


derived by Hylleras (Zeits. f. Physik 74, 216 (1932)). 
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where » is the frequency and //, the n Hermite polynomial. Introducing Weber's function 
D,,.(y):"8 


H,(2-"2y) = 2”/%e 4D, (y) 
the right-hand side becomes: 


(20)'/2¢-0/2 > e"*D,,2(2! 2y)/n!. 


We distinguish again two cases; 
a. Boltzmann statistics. We have to sum over all values of ». From the general result:!7 








D,.()Dx(n) ne £2 + ny? — 2Ene*) ss 
———_——_e "= i— e= | 2e(e-+n-))/4 2xXD ne ee => 35 
~ n’ ( ) _ a1—e) J (35) 
putting x= 0; £==2"/?y, we obtain: 
. 20 6 6 
(2rmkT)—V2 Ye FnthTy,2 -( 5) exp ( —— — y* tgh- -). (36) 
= l-—-e¢? 2 2 


One sees that for @ very small the right-hand side goes over into: 


Putt m= g-tetme%e2/AT = gV IAT, 


b. E.B. or F.D. statistics. We must suppose then that we have to do with two particles 
attracting each other with a force proportional to the distance between them, moving in a one- 
dimensional vessel. If x stands for the distance between the particles, in E.B. statistics we must 
take only the eigenfunctions even in x, and in F.D. statistics we must take the odd ones. Separat- 
ing into center of gravity and relative motion coordinates,'® we have then to consider for the 
relative motion simply the sum (34) where in E.B. statistics we have to sum over all even values 
of , and in F.D. statistics over the odd values’? and where one must remember to replace m by 
3m. Using now the general result ;?° 


s} é “2s 
cosh { ne =} : (37) 


sinh 1 —e 





-nz = (1 — ee) texp ( ion +9 cotgh r) 


n 


a D,(#)Dz(n) rn? 
ose eene 


From this follows then immediately: 


kT)722> EB. re 0 ° a 
(xmkT) 2EP 8 i(, 20 ) exp (— > - s*tsh +): 41 + exp (-=—) (38) 


(xmkT)-/2 > Fp. — ¢%# sinh @ 


For small @ the right-hand side goes over into: 
Rene my?z? kT (4 + ent mkT 2? jh?) 


in accordance with (8). 


16 Whittaker-Watson, Modern Analysis, p. 347. 

17 See G. E. Uhlenbeck and L. S. Ornstein, Phys. Rev. 36, 823 (1930). F. Bloch reference 11 
also gives the result, using a transformation function of Kennard (Zeits. f. Physik 44, 326, 1927.) 

18 Strictly this is only possible for an infinite vessel. Because we are not interested here in 
the discreteness effect of the total translational motion, this is quite allowed. 

19 The center of gravity motion will simply give an extra factor (4amkT/h?)"?. 

20 The proof is quite analogous as given for (35) by G. E. Uhlenbeck and L. S. Ornstein, 
reference 17. 
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On Ferromagnetism and Related Problems of the Theory 
of Electrons* 


By Paut S. Epstein 


California Institute of Technology 
(Received May 20, 1932) 


In this paper is presented a study of those electrons which hold together a 
chemically elementary crystal and are responsible for the homoeopolar bonds be- 
tween neighboring atoms. The starting point of the theory is Slater’s form of the 
secular equations for first order perturbations arising from the interactions of iden- 
tical atoms composing a system. A rigorous solution of these equations is given in 
section 4. The investigation is carried through in a unified way for magnetic and 
non-magnetic materials, with respect to internal energies and to magnetic proper- 
ties. As anticipated by Heisenberg, the deciding factor is the sign of the Heitler and 
London interchange integral J,. Substances with large negative J; are non-magnetic, 
those with large positive J;, potentially ferromagnetic. For non-magnetic bodies, 
the theory gives a confirmation of Bloch’s conclusions with slight differences of inter- 
pretation. On the other hand, the results for magnetic materials are new and entirely 
different from Bloch’s. As to specific heats, it is found that at very low temperatures 
they have the expression c =0.208 sR(7/6)*/?, where R is the gas constant, s the num- 
ber of valency electrons per atom, and @ has a close relation to the Curie point. With re- 
spect to ferromagnetism, the result is that a crystal satisfying Slater’s equations is 
spontaneously magnetized almost to saturation but that the polarity of this magnetiz- 
ation changes its sense at irregular intervals. This fact suggests that ferromagnetic 
crystals must have a block structure and that they are coherent, in the sense of the 
validity of Slater s equations, only within the blocks (compare sections 10 and 11). 
With this hypothesis the theory accounts for the fundamental facts of ferromagnet- 
ism. 


1. INTRODUCTION 


HE discoveries of quantum dynamics, especially Pauli’s exclusion prin- 

ciple and Fermi’s statistics based on it, gave a new stimulus to the theory 
of electrons in metals. While Sommerfeld! treated in an exhaustive way the 
subject of free electrons, it was pointed out by Heisenberg? that the exclusion 
principle must exercise a powerful influence on the orientation of the spins 
of valency electrons and that it may play the role of the molecular field postu- 
lated by Weiss for ferromagnetic materials. The theory of the more tightly 
bound electrons was investigated by Bloch in a series of important papers. 
While his work is fundamental in many respects, it treats the non-magnetic 
properties’ and the magnetic‘ separately and by different methods. Moreover, 
we believe that his results relating to ferromagnetism are open to objection. 


* An abridgment of this paper was read before the National Academy of Sciences on April 
26, 1932 at Washington, D. C. 

1 A. Sommerfeld, Zeits. f. Physik 47, 1, 43 (1928). 

2 W. Heisenberg, Zeits. f. Physik 49, 619 (1928). 

3 F. Bloch, Zeits. f. Physik 52, 555 (1928); 59, 208 (1930). 

4 F. Bloch, Zeits. f. Physik 57, 545 (1929); 61, 206 (1930); 74, 295 (1932). 

5 This assumption is made for simplicity but is not necessary. We shall see in section 5 
that our theory applies, almost without change, to more general cases. 
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2. FORMULATION OF THE MATHEMATICAL PROBLEM 


We consider with Heisenberg and Bloch a system of V electrons each per- 
taining, in the unperturbed state, to a different atom of a crystal.’ Let the 
orbital motion of the electron be that corresponding to the lowest quantum 
state. Besides, with every electron there is associated a spin the projection of 
which on a given axis can be either positive or negative. Let 7 electrons have 
negative spins and V—r positive. The projection of the total spin moment is 
then equal to mh, 27 ,where 2m = N—2r. If all the atoms are successively la- 
beled by the numbers 1, 2, - - - NV, let the positions of those containing an 
electron with a negative spin be fi, fo, - - «fr i<fe< --- <f,). We denote 
the antisymmetric wave function for such a distribution by (fi, fo, - - - f,). 
As the magnetic interactions are neglected, all the functions with the same 
number r, correspond to the same energy and the total wave function per- 
taining to the magnetic moment m can be written as the linear expression 

Wir) = DV alf,-- fh, +: f,). (1) 
Phot 

The coefficients a(fi, - - - f,) are determined, in the usual way, by a secu- 
lar system of equations when we introduce the mutual influence of the elec- 
trons and protons as a perturbation. The form of the secular equations was 
explicitly derived by Slater * 


ea(fi,-- fi) + Dy lath’, >> +f’) — af, --f)] = 0. (2) 


The summation is extended over all pairs of electrons having opposite 
spins and f,’, -- -f,’ and fi, - - - f, differ only in that this pair of spins is in- 
terchanged. The coefficient J;;, denotes the London and Heitler interchange 
integral for the pair and € is the energy measured from an arbitrary zero point. 
Usually, it is permissible to neglect the mutual influence of the more distant 
atoms and to take into consideration only the interactions between contigu- 
ous electrons. In this case the Eq. (2) is for most crystals reduced to 


ea(fi, fy ta dD lan’, -+ +f) — afi, > --f)] = 0. (3) 


The sum includes now all pairs of contiguous electrons having opposite 
spins. J, represents the interchange integral for two adjacent atoms. 

The main problem is to find a solution of the Eqs. (2) or (3), satisfying the 
conditions at the surface of the crystal and other subsidiary conditions which 
we shall state in the next section. Bloch has given a very elegant method of 
transforming these equations. However, we shall not use it because the solu- 
tion can be derived directly from the form given here. 


3. SOLUTION IN A SPECIAL CASE 


In order to show how the solution is arrived at, it is best to start from an 
example. We take, therefore, the case of a linear chain of atoms with r=2. 
The form which Eq. (3) takes in this case is 


6 J. C. Slater, Phys. Rev. 34, 1293 (1929); 35, 509 (1930). 
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ea(fi, fo) + Jila(fi — 1, fe) — a(fi, fe) + a(fi + 1, fe) — alfa, fo) 
+ a(fi, fe — 1) — alfa, fe) + a(fi, fo + 1) — a(fr, fe)] = 0, 


which holds for fe>fi+1. In the case fo=f;+1 two of the interchanges are 
precluded and the equation becomes 


ea(fi, fo) + Jila(fi -i,f)- a(fi, fe) + a(fi, fo + | a(fi, fo) | = 0). (4b) 


The difference equation (4a) has constant coefficients and its particular 
solutions can be written in the form exp ia(kifitefe). The corresponding 
energy level is 


(4a) 


€ ko = 4J,(sin® Sak; + sin? ake). (5) 


This expression is symmetrical in k; and k, so that the solution exp ia(kof; 
+kife) belongs to the same energy. We can, therefore, build up the slightly 
more general solution a(ki, ke) =c, expia(kifitkefe) +c expta(kefit+kife) and 
try to choose the coefficients c; and ¢ in such a way as to satisfy the subsidiary 
condition (4b). It is important to note that while (4a) is supposed to apply 
only in the case fo>f/;+1, our solution is of such a type that it satisfies the 
Eq. (4a) for all values of f; and fe, including fe=f;+1. In this last case, both 
Eqs. (4a) and (4b) are true and we can replace the subsidiary condition by the 
difference 

a(fi + 1, fe) + a(fi, fe — 1) — 2a(fi, fe) = 9, 
fo = fi +1. 


It is easy to see that this condition is satisfied if we put 


ciice = — exp ila(ki — ke) + 26(hki, he), 


(6) 


where 
sin ak; — sin ake 7 
o(ki, k2) = arctg — — * (7) 
2 — cos ak; — cos ak, 





The solution itself can then be written in the form 


a(ky, ke) = exp ila(kift + kof) + ak; + (hi, ke) | 


(8) 
— exp ila(kofi + kife) + ake + o(he, ki). 


If we substitute for the surface conditions the periodicity requirement of 
Bloch’s, all we have to do is to define ki, k, as integers and a=27/N. The 
Eqs. (8) and (5) will then represent a rigorous solution of the problem. It 
differs from Bloch’s result only inasmuch as a(ki, ke) vanishes for k; = ke. This 
means that the case of two equal quantic numbers does not occur and must 
be excluded. 

To obviate possible objections, we prefer, however, to show how a solu- 
tion satisfying the actual border conditions can be found. These conditions 
result from the fact that when an electron is at one of the ends of the chain 
(f1:=1 or fe=N) a further decrease (or increase) of the number f labelling its 
position is no longer possible. The corresponding transitions must be, there- 
fore, excluded from the conditions (4a) and (4b). This means 
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a(fi — 1, fe) — a(fi, fo) = O, when f; = 1, (9a) 
a(fi, fe + 1) — a(fi, fe) = O, when f2 = V. (9b) 


It is easy to see that in addition to a(ki,k2) the expressions a(—hi, ke), 
a(ki, —k2), a(—ki, —k2) are also solutions of the system (4a), (4b) pertaining 
to the same energy (5), because this energy does not change if we reverse the 
sign of k; or ko. Combining these four expressions with suitable coefficients 
we can build up a solution satisfying the condition (9a). In fact, the Eq. (9a) 
is fulfilled if the expression is arranged as a sum of terms each containing 
the factor cosak;(f;— 3), where j is 1 or 2. Such an arrangement is possible for 
any value of the parameters ak;. On the other hand, the border condition 
(9b) can be satisfied by the same expression only for definite and discreet 
values of these parameters. It is obvious from the symmetry with respect to 
the two ends of the chain that a must be taken as 

T 


a=—. 10) 
N (10) 





The condition (9b) is equivalent to the requirement that the solution 
can bearranged as asum of terms each containing the factor cosak; (fe—N—}). 


The arguments of the cos functions can be reduced to this form only if we 
impose upon the numbers &; the conditions 

k,’, 

“ (11) 


Where the numbers k,’, ke’ are integers (1, 2, - - - V) and represent the 
quantic numbers of the problem. The rigorous border conditions lead, there- 
fore, to a different result from Bloch’s periodicity relation: The numbers k; 
are no longer quantic numbers but non-integral auxiliary parameters. 

The coefficients of the secular problem acquire in this case the simple form 


a(fi, fo) = cos ak\(fi — 4) cos ako(fe — N — 3) 


F cos ako(fi — 3) cos aki(fe — N — 3). 


(1 — 1/N)ki — (1/m)[$(hi, ke) + (kr, — he) ] 
(1 — 1/N)ke — (1/m)[G(ke, hi) + (he, — fi) ] 


4. GENERAL SOLUTION 


The results of the preceding section apply with slight changes to the 
general case. For a system with r=}N—m negative spins the energy expres- 
sion is quite generally 

4N—|m| 


cn = >» €j- (12) 
j=1 


In the case of a linear chain, we have 


ak; 
, (13) 


¢ = 4], sin? 





T 


4 
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where the parameters k; are determined in relation to the quantic numbers 
k;’(=1,2, - - - N) by the equations 





1 1 
(1 - ak — — Dilek, ki) + o(ki, — kd] = &;’, 
i WT ixzj 
So, (14) 
i,j =1,2,---,3N—| m|. 
The expression for a(fi, - - - f-) becomes rather cumbersome. It is built 
up of terms derived from the product 


cos ak, (fy = 3) -cos [ako( fo we 3 2) a o(ke, ky) aad (ke, at ki) | 
-cos [aks(fs — 3/2) — (ks, ki) — o(ks, — ki) — O(ks, he) — (ks, — he) ] + + - 


by all possible permutations of ky, ke, - - « k,. 

If the number N is large, the difference between k; and k;’ does not ma- 
terially affect the energy expressions for smaller values of r, i.e., in the neigh- 
borhood of magnetic saturation. As to the opposite case r= N/2, m=0, the 
terms of the energy expression in which k; is small and those in which it is 
large (near to V) remain practically unchanged. But the middle range of 
values of k;, about V/2, is considerably influenced in the sense of an increase 
of this number. In the more important of our applications the states of small 
magnetic moment will play but a negligible role. Therefore, we shall be satis- 
fied with the approximation which is obtained by identifying k with k’, as it 
is amply sufficient for our purpose. 

For a two dimensional quadratic lattice 


_ . wh; al; ; : 7 
e; = 4J,{ sin? — + sin* — }; G? = N. (15) 
2G 2G 
In the case of a three-dimensional cubical crystal 


ml 





_ . wh; ; j oN; ; ; 
€; = 4/,| sin? —— + sin* —— + sin? —}; G* = N. (16) 
2G 2G 2G 


The parameters k;, /;, 2; are in both cases connected with a set of quantic 
numbers k;’, 1;’, n;’=1, 2, - - - G by the relations (14) (with G in place of V) 
and by two similar sets of equations which are obtained from (14) by substi- 
tuting / or instead of k. 

In all cases, the Eq. (13) can be divided into a main condition, analogous 
to (4a), which holds when no two negative spins are in contiguous atoms, and 
into subsidiary conditions like (4b) which are true in the opposite case. The 
main condition is satisfied by an exponential function (analogous to expa(ki/fi 
+hkofs) of the preceding section) and by a large number of other exponentials 
derived from this. The form of the energy terms in the sum (12) is here also 
completely defined by the main condition and by the exponential function 
satisfying it. The subsidiary conditions are satisfied by a proper combination 
of the exponentials and do not affect the form of the energy expressions. 
Finally, the border conditions serve to define the parameters entering into 
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the energy terms. Bloch succeeded in solving the main condition, therefore, 
his energy expressions are formally identical with ours, although the symbols 
have a somewhat different meaning. To this extent his results for other types 
of two and three dimensional lattices can also be used. However, if we inquire 
into the substance instead of the form, there is a profound difference in two 
other respects which are far more important than the changed definition of 
the parameters. 

(1) In the first place, no two terms of the energy sum (12) may be equal. In 
the example of section (3) we see from the formula (8) that a(fi, fe) vanishes 
when ki =k. This means that the corresponding energy level does not exist. 
The conditions are the same in the general case. Whenever two triples of 
numbers k;, /;, n; are equal the function a(fi, - - - f,) vanishes identically. 
This is by no means surprising because the exclusion principle is contained in 
the formulation of Slater’s Eqs. (2) or (3) and makes itself manifest in all the 
conclusions drawn from it. 

(2). In the second place, the number of terms in the sum (12) depends only 
on the absolute value of the magnetic moment. It is equal to }.V— | m|, and 
not to }N—m as Bloch has it. The mathematical reason for this is that the 
case of no two negative spins being in contiguous atoms can be realized only 
as long as the number of these negative spins is smaller than (or equal to) 
half the number of atoms (r=}.N—m S3.N, or m20). If this requirement is 
not satisfied, we have no main condition but only subsidiary conditions, which 
completely changes the aspect of the problem. However, we can reduce the 
new problem to the old one by reversing the direction of our system of coor- 
dinates, and by counting as negative those spins which were considered as 
positive. 70 positive and negative magnetic moments of the same absolute mag- 
nitude corres pond the same energy levels. This fact is an immediate consequence 
of the fundamental assumptions of this theory (section 2) which neglects 
magnetic interactions. Both directions in space are equivalent and the energy 
cannot depend on sense of the magnetic moment. 

These two features of our energy expressions have an important influence 
on the results and lead to conclusions which are considerably different from 
Bloch’s. Before closing this section, we should like to point out that the theory 
can be carried through also in the case of Slater’s general Eq. (2) which takes 
into account interactions between non-contiguous atoms. For instance, in the 
case of a cubical crystal the energy expressions become 

€5 = 4 DS So oog{sin? Jagik; + sin? dagel; + sin? dagsn;|. (18) 
Y 19203 
where J,,,,,, is the interchange integral for two atoms with a distance be- 
tween them having the projections gid, ged, g;3d if d denotes the lattice con- 
stant. 
5. THE SUM OF STATES 


To the energy expressions (12) there must be added the magnetic energy 
of the crystal equal to —2um//, where uw denotes the Bohr magneton and // 
the strength of the magnetic field. The total energy is, therefore, 
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€ = €m — 2umH. (19) 


The part €,, depends only on the absolute value of the moment so that 
€_m=€m. For every energy level (19) there exists, therefore, another of the 


form 
€= €, + 2umH, (19’) 


which corresponds to the magnetic moment —ym. 
Accordingly, the sum of states Q = Lexp(—€/kT) can be divided into two 


parts 
Q= 2H) + 2(- 4), (2)) 


the first containing all the terms with exponents of the type (19), the second 

all those of the type (19’). The two halfs Z(/7) and Z(—J//) differ only in the 

sign of /7 so that the sum of states Q is an even function of this quantity. The 

average magnetic moment of the system 

_d log Q ap 2) — Z'(— A) 
dH Z(H) + Z(-— H) 


M = kl 








(21) 


becomes an odd function of H which vanishes for J7=0. This fact seems 
to make this theory inadequate for the explanation of permanent magnetiza- 
tion. We believe, however, that this average is obtained by a type of statistics 
which has a limited physical reality. In the case of ferromagnetic substances 
the averaging must be done in a different way. We shall return to this ques- 
tion in section 10 and shall first develop a few equations for the evaluation of 
the function Z(77). With the help of the formula (12), we can represent each 
term of the sum of states as [5., exp( —€)/kT+2mpll/kT), where 2m = N—2r. 
But according to (13) or (16), €; is a function of the parameters &;, /;, ; and, 
for a given m, there will be as many different terms as there are possible 
choices of these numbers. The totality of them is exp (2mul1/kT)Z,, 


Z,= DIT exp (- «,/k7). (22) 
kj; j= 

We use here the subscript k; as an abbreviation for the triple of param- 
eters k;, /;, n;. As we found in the preceding section, all the €,; in the sum 
(12), or all those in our product, must be different. In other words, our prob- 
lem is subject to Fermi statistics. Our formula can, therefore, be treated in 
the usual way: If we consider the expression 


R =|J[1 + 2 exp (— «,/k7)], (23) 
kj 


then, Z, is obviously the coefficient of the term of the degree r of the expan- 
sion of R in powers of z: 
1 R 


Zz,=— dz, 
QriJ ott! (24) 





the path of integration being a circuit round the point z=0. 
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The half Z(//7) of the sum of states is now 
N/2 
Z(M) = > exp [(N — 2r)ull/kT]Z,. (25) 


Strictly speaking, the expressions R are not the same in the different 
functions Z,. But within the approximation of putting k; equal to k,’, ete. 
(section 4) we can regard them as identical. Substituting (24) in (25) and 
summing under the sign of the integral, we find 


1 R (se7)i4+! — J 
Z(H) = =| - — -———— — ds, (26) 


Qri giNtl ze? — 1 
where 7 is an abbreviation for 
7 = Qull/kT. (27) 


Although we spoke of one electron in every atom our equations apply 
as well to the case of several valency electrons of different energy. In fact, 
each of them will be in resonance only with its own kind in the other atoms 
and will satisfy a separate secular system. The total sum of states is then the 
product of the sums of state (26). If there are in an atom two or more electrons 
in nearly equivalent orbital states and with unlinked or weakly linked spins, 
the conditions are somewhat different; they are being now investigated by 
Mr. D. Weinstein. It is probable that our expressions will continue to be 
true, also in this case, with a change of the physical meaning of the sym- 
bol J; which will now include both the external and the inner interchange 
energies. 


PARTIAL EVALUATION OF THE SUM OF STATES 


We consider the case that €;/k7 is, in the average, large in absolute value 
compared with the quantity 7. In fact the characteristic feature of our theory 
is Heisenberg’s assumption that there is a strong inner directing force due to 
the energies €;. The case of small €; is that when this assumption is not justi- 
fied; it is, therefore, uninteresting in this connection. We have two terms un- 
der the integral (26), let us focus our attention on the second and evaluate it 
by the method of the steepest descents. In the usual way, we put the inte- 
grand equal to e® and obtain the position of the saddle point z=A from the 
condition d¢/dz=0 

{ exp (— i/ kT) A 
ia ——_———— — (}N + 1) - ———_ = O. (28) 
1+. narmags kT) A-e¢* 

The second derivative we denote by 


exp (— €; kT) > N+1 1 
ee oe ie 


1+ A exp (— €; kT) 








as 29 
ear 


The position of the saddle point depends on the sign of the energy e€; or 
of the interchange integral J;. If the sign is negative, A is small compared 
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with exp (—7): In fact, the left side of the Eq. (28) is negative when A =0, 
but since the exponentials are large, most terms of the sum are, practically, 
equal to one already for moderate values of A. The sum becomes equal to V 
and the expression positive. There is a root in between, which does not lie 
close to the point z=e~*. Therefore, the last term in (28) is negligible com- 
pared with the first two and must be omitted. Neglecting, at the same time, 1 
beside }.V, we write it in the form 


> = 


1+ A(— €)/kT) 





_ A exp (— €/kT) cd (28’) 


Turning to the first term of the integrand of (26), we see that it is regular 
within the path of integration. It has no pole in s=0 and the point s=e™ lies 
outside the circuit of integration. The second term represents, therefore, the 
whole integral which takes the form 


1 
Z(H) = ——— [][1 + A exp (— ¢€)/kT)]A~¥?-(1 — Aer)". (30) 
(2rqg)'/? x, 


It must be noted that A does not depend on the magnetic field because 
the term containing 7 is neglected in the Eq. (28’) Z(/7) depends on //, only 
inasmuch as it enters into the last factor of the expression. 

The results are quite different in the case of a positive energy €;. A be- 
comes now much larger than e~*. In fact, for A =~ the left side of the Eq. 
(28’) is positive. But for comparatively large values of A the terms of the 
sum become negligible so that the other terms dominate and make the ex- 
pression negative. The integral retains exactly the same form (30) but it does 
not represent now the whole function Z(//). The circuit of integration is 
widened to pass through the distant point s=A and the point z=e” is now 
within this circuit. Consequently, the first term of the integrand (26) has 
now a pole in this point and the integral does not vanish but has the value 
Z(H): 


ZA) = ZN) + ZOU), 
where 


N 
Z(H) = exp a TI [1 + exp (— + — ¢/k7)] (31) 
kj 


and Z® is identical with the expression (30). It is obvious , for physical and 
for mathematical reasons, that at low temperatures Z™ is negligible com- 
pared with Z™. In fact, apart from the unimportant factor (1—Ae’)", 
Z (IT) is identical with the expression Zy,, of the last section which repre- 
sents the totality of all terms of the sum of states corresponding to m=0. It 
follows from our expression of the energy (12) that the exponents of these 
terms are very large, because all the terms of the sum (12) must be differ- 
ent and, consequently, only few of them can have very small values. For this 
reason Zy,2 is negligible, if J,/T is sufficiently large. 
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7. EVALUATION CONTINUED 


We now substitute into our general expressions the formula (16) for the 
energy term. It will be sufficient to treat the case of the three dimensional 
lattice as the other cases have no physical interest. We use the approximation 


k;, 1;, n;=1, 2, - - - G which is good enough as was shown in section 4. If we 
define 
tk; ml; TN; 32) 
i om sone on > = — 
— aes Th 5 = oi (32 
2G 2G 2G 


the summation can be replaced by an integration with respect to £&, n, ¢ be- 
tween the limits 0 and 7/2. Moreover, we use the abbreviations 


B = 2/,/kT, cos 2 + cos 2n + cos 2¢ = p. (33) 
The terms of the formula (16) become then 
€; = 2BkT(sin* — + sin? n + sin? ¢) = BRT(3 — p). (3!) 


Turning to the condition (29’) we can include the first term 33k7' into the 
parameter A writing A’=A exp (—38). We obtain, then, 


8 r/2 A’ exp Bp 
saa ff J didndg = 3. (35) 
r 0 1+ A’ exp Bp 


The constant A’ can be determined from this relation in the following way. 
Let us subtract both sides of the equation from 1: The right side remains 
equal to 3 while the left side has an integrand with the same denominator as 
(35) but with the numerator 1. The equation can, therefore, be represented as 


8 #2 1/A’) exp (— Bp) 
— ff : tel det. te — didndé = 3. 
rJJo  J°1 + 1/4") exp (— Bp) 


We substitute now &=}32r—£&’ and similarly for n, ¢, so that p= — (cos 2’ 
+cos 2n’+cos 2¢’) = —p’, obtaining an equation which has quite the same 
form as (35) except that A’ is replaced by 1/A’. We conclude from this that 
A’'=1/A’' or 








A’ =1, A = exp 38. (36) 


If we omit in (30) the unimportant factor 1/(27q)!, Z(/7) takes in the case 
of negative interchange integrals J; <@ the form 


8\ 
log Z(H) = — (3/2)8N + — 
- 


w/2 
: ff fe [1+ exp B (cos 2 + cos 2n + cos 2¢) |dtdndeé 
0 


37 
— log [1 — exp (38/2 + 7)]. (37) 


Equally simple are the expressions in the case of a positive interchange 
integral (J; >0). Substituting the notations (32) into the expression (31) we 
find 
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8.V 


log Z(H) = 3Nr + — 
= 


tl 


. (38) 
ff flog [1+ exp (— r — 28(sin? + sin? » + sin*f)) |dédnd¢, 
0 


to this must be added negative term Z®? of the form (37) which is negligible 
for low temperatures. 


8. ELECTRONIC SPECIFIC HEATS 


If the crystal is not subject to magnetic forces (J7 =r =0), the two halfs 
of the sum of state Q become identical and Q=2Z(0). As the free energy is 
connected with Q by the relation y = —kT log Q, it is easy to find the expres- 
sions for the internal energy and the specific heat. 

We shall see in the next section that the case J; <0 is that of non-magnetic 
(diamagnetic) bodies. It we denote, in the usual way, kN =R and, moreover, 


B=-YT, @= — 2),/k, (39) 


and the number of valency electrons in an atom by s, we find for such materials 
the temperature dependent part of the internal energy 


8 w/2 lidndt 
u = sro, ff i; 2a (40) 
id 0 1 + exp (— Bp) 


This expression is but slightly different from Bloch’s and agrees with it 
in every important respect giving results analogous to those obtained by 
Sommerfeld for free electrons. For high temperatures we can expand with 
respect to 8p taking only terms of zero and first order. This gives 





log Z(0) = Ns(log 2 — 38/2 + 387/16), (41) 
the temperature dependent part of the internal energy 
U = — 3sR0?/8T, 


and the specific heat 
c= dU/dT = 3sR(0/T)*/8. (42) 


The evaluation for low temperatures can be carried out in a similar way 
as with the integrals of Fermi and Sommerfeld. The zero order approximation 
is obtained by taking the integral /ffpd&dndt over the part of the cube 
0<é<7/2, etc., limited by the surface p=@. This term which, in every 
case, constant and of little interest happens to vanish. For the next ap- 
proximation it is well to remember that the normal distance between 
neighboring surfaces p=cos 2&+ cos 2n+ cos 2{=const. is equal to 
|(dp/d£)?+(dp/dn)? + (0p/d¢)?|—Ap =} [sin? 2E+ sin? 2n+ sin? 2¢]-!Ap. Near 
p=0 it can be expressed by }|3—2(cos 2n+ cos 2¢)?+2 cos 2 cos 2¢]~!Ap. 
The factor of Ap varies between a minimum of 1/2-3! and a maximum of }. 
As the lower values have a larger weight, we take as a rough average § of 








102 PAUL S. EPSTEIN 


the minimum plus } of the maximum which is very close to }. On the other 


hand, we must know the area of that part of the surface p=0 which lies 
inside the fundamental cube. This surface goes through the middle of six 
edges of the cube which lie at the corners of a plane regular hexagon and also 
through the center of the hexagon. The surface itself is not quite plane and 
its area is, therefore, slightly larger than that of the hexagon. We take as 
its rough approximation the area of the circle circumscribed about the hexa- 
gon which is equal to 7*/8. Multiplying the area by the thickness we find 
that in the neighborhood of the element of volume element takes the form 
dr =r'°dp/32 and the correction of the integral becomes 


1 [— as (43) 
o i+e% 246? 
The internal energy for low temperatures is, therefore, 
U = (x*/24)sRT?/0, (44) 
and the specific heat 
c = (nr*/12)sRT/0. (45) 


The expressions are different in the case of magnetic materials (J; >0) and 
they are derived here for the first time. We have to use the expression (38) 
with t=O and find 


8 7/2 28(sin? € + sin? n + sin? ¢)dédnd ¢ 
re J Jo 1 + exp 28(sin? & + sin* y + sin*f) 





For high temperatures (81) we can again expand with respect to the 
exponent and obtain 


21 
log Z(0) = Vs( tog 2 — 36/2 + ft), (47) 


while —Z® (0) is given by the expression (41). We see that for high tempera- 
tures Z® is not negligible compared with Z“ but becomes of a similar order 
of magnitude, especially, as the factor 1/(27q)! is no longer small in this case. 
We introduce the slightly changed definitions 


B=0/T, 6 = 2J,/k, (48) 
If we determine U and c for 7>@ from Z™ alone, 


U = — (21/8)sRe°/T, ¢ = (21/8)sR(0/T)?, (49) 


we get the correct dependence on the temperature but coefficients which are 
slightly too large. 

The evaluation for low temperatures 7<6@ is facilitated by the fact that 
the exponent retains its positive sign in the whole interval of integration. 
If 6 is large, the integrand is appreciable only for small values of the vari- 
ables &, n, ¢. It is permissible to replace the signs by their arguments which 
lead to 
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3 41 1 1 
U = -——_—_|1 - —__+—.... |sar, 
2 (2x8)3!2 25/2 © 35/2 


U = 0.0830sRT®!2/93!2, (50) 


or 


and 
c = 0.208sR(T/6)°”?. (51) 


The coefficient in this formula is rigorous and not approximate as in (45). 

The measurements of Eucken and Werth’ on the specific heats of iron 
and nickel show for temperatures between 16° and 22°K considerable devia- 
tions from the third power law in the sense of an increase of heat capacity. 
Eucken himself points out that his observations did not go to sufficiently deep 
temperatures to establish the law which this excess follows. All we can say is 
this: If the deviation is attributable to the cause discussed in the preceding 
paragraphs, its order of magnitude leads to a value of @ of a few hundred de- 
grees and a value of J; of a few one hundredths of a volt. If the formula (51), 
for magnetic materials, were confirmed, it could give us very valuable in- 
formation. The knowledge of the characteristic temperature @ would be very 
helpful in elucidating the mechanism of ferromagnetism (compare the next 
two sections). 

It should be pointed out that our theory neglects the possible variability 
of J; and @ with temperature. This is not a rigorous procedure because J, 
must depend to some extent on the distance between the atoms. We pre- 
sume, however, that this source of error is very small at low temperatures be- 
cause of the insignificance of the thermal expansion in this region. Neither 
is the energy of magnetization contained in our expressions but this cause 
makes an appreciable contribution only in the vicinity of the Curie point. 


9. MAGNETIZATION 


We turn now to the behavior of our crystals in a magnetic field. The mag- 
netic moment is obtained by differentiating the logarithm of the sum of 


states Q with respect to // 
d log Q 


i oo iT 52 

M = k7 os (52) 

Let us first discuss materials with negative interchange integrals (J; <0, 
8<0) for which we have to apply the expression (37). The dependence of 
Z(H) on the strength of the field is all contained in the last factor [1—exp 
(—38/2+17)]— which is very small compared with the other factors. In 
fact, it does not contain the number of atoms N and this makes its logarithm 
entirely negligible beside the logarithm of the rest of the expression. Sub- 
stances with a negative J; are, therefore, magnetically inert as far as the 
electronic interchanges are concerned. \We must remember, of course, that our 
formulas are valid only when 8 is considerably larger than 7, or —J, larger 
than 2uJ7. For fields of about 10,000 gauss this would mean that —J; must 


7 A. Eucken and H. Werth, Zeits. anorg. Chemie 188, 152 (1930). 
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be larger than 10-4 volt. We have pointed out in section 6 that the assump- 
tion of very small values for J; has little interest as it corresponds to the 
case when Heisenberg’s hypothesis does not apply and the older theory of 
paramagnetism enters into its rights. The fact that materials with negative 
interchange integrals of appreciable absolute value must be nonmagnetic is 
evident directly from the energy expressions (12), (16). The state of vanishing 
magnetic moment (m =0) is, then, that of the lowest energy which dominates, 
in an overwhelming way, over all other states. We arrive, in this way to the 
following classification of crystals with respect to their interchange integrals 
J,. Materials with large negative values of J; are nonmagnetic as far as elec- 
tronic interchanges are concerned. Materials with small values of J; (negative 
or positive) are paramagnetic. Those with large positive values of J; are, 
potentially, ferromagnetic (compare with next section). 

In the next section we shall give reasons for the view that the formula (52) 
should be applied with restrictions to the case of magnetic materials (J; >0). 
Here we wish only to investigate what sort of a magnetization curve would 
result if this equation were always valid. We know from formula (21) that 
this curve is antisymmetrical of the general character of those pertaining to 
paramagnetic bodies. What is its slope in the point /7=0? It is easy to see 
that in this point the slope is given by the expression 


en ee eae 
di/ dH? dil H=0 


We shall be satisfied with low temperatures when 8 becomes large. The 
Eq. (38) gives us then the approximation 


T \3/2 ete 3? ; 
log Z(H) = 3Nsr + vs(5—) Je pct avo ee |. (53) 


270 


(dM/dH)) = (Nsu)*/kT, (54) 


neglecting a term of the order N. If we refer / to one/cm* of the material, 
the expression must be multiplied by n/N (n being the number of atoms 
per/cm*) and becomes nNs*u?/k7. The values of the constants are as follows 
uw=0.9-10-°°, r is of the order 10”, kT is about 4-10-'4 for room tempera- 
ture. (d.M//dI/), is, therefore, of the order of magnitude 2s? -10~*. This repre- 
sents a steep rise of the curve even for very small crystals. 


10. EXPLANATION OF PRIMARY FERROMAGNETISM 


We shall see in the next section that, in a sense, our Eqs. (52) and (54) 
describe the magnetization curve of certain systems. However, their physical 
reality is derived in an indirect way: We believe that the method of the sum 
of states is inapplicable to a coherent crystal in the conditions which prevail 
at low temperatures. In fact, what is the physical meaning of the sum of 
states Q= > exp (—e€/kT)? It is derived from Gibbs’ concept of the ensemble 
of systems which is equivalent to the succession, in time, of the different 
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states of one given physical system. The mean value of any quantity derived 
by the use of the sum Q represents, therefore, the time average of this quantity. 
In the cases in which this method has been heretofore used this time average 
is what we really observe. In the theory of gases states which occur with an 
appreciable probability deviate but little from the mean state. The time aver- 
age is very close to the most probable state of the system and very close to, 
or identical with, the readings of our instruments which always integrate 
over a finite, though sometimes small, time. The situation is entirely differ- 
ent in our theory of magnetization. It is true that positive and negative mag- 
netic moments will occur equally often in the course of long periods and the 
time average will vanish. But this time average is not at all close to the most 
probable state. On the contrary, under the conditions mentioned, it coincides 
with one that is so improbable that it practically never occurs. To show this, 
let us compute the time average of the absolute value of the magnetic moment 
in the absence of a field. 





= kT(d log Z(H)/dH)o = Nsu[1 — 0.058(7/6)*/?]. (55) 


) 
| <um | 


This means that, on the average, the crystal is in a state close to magnetic 
saturation. When we further compute the relative mean quadratic deviation 
from this time average 





(| m | —|m)\*/| m\*? = (0.37/Ns)(T/6)*?, (56) 


we find that the spread around the mean value is extremely small. These 
formulas show us that the crystal is, practically, always in a state approach- 
ing saturation. The time average zero is explained by the fact that the two 
senses of the magnetic moment are equally probable so that the crystal, at 
irregular intervals, spontaneously changes its magnetic polarity. 

These results refer only to the case of low temperatures. When the tem- 
perature becomes sufficiently high, the state of vanishing magnetic moment 
is no longer one of low probability. This is apparent from our expression (41) 
valid for high temperatures. This expression is, practically, identical with 
the logarithm of Zy,s of section 6 representing the totality of all the terms of 
the sum of states pertaining to the magnetic moment m=0. The Eq. (41) 
shows us that this term has a negative exponent and is small as long as 
38/2 —38?/16>log 2. But when the temperature rises the inequality is re- 
versed into 38/2—36?/16<log 2, the exponent becomes positive, and the 
term large. It is true that, in this form, the argument is not quite convincing: 
The formula is derived for small values of 8 and it is not certain that it still 
holds when 38/2 — 36/16 becomes equal to log 2 or smaller. However, we get 
qualitatively the same result when we go back to the original definition 
Zyj2= > exp (—e/k7) the summation being extended over all possible energy 
levels compatible with the condition =0 or r=N/2. Now, let us look at 
our energy expressions (12) and (16). The number of terms in the sum (12) 
depends on the magnetic moment 2m; it has its minimum 0 for saturation 
(2m=N) and its maximum N/2 for the nonmagnetic state (2m =0). As all 
the terms of the sum must be different, the mean energy also increases from 
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e =(), for saturation, to a huge value, for the conditions of vanishing magnetic 
moment. It is easy to see that the energy excess of the state 2 =0 over the 
state 2m =N cannot be smaller than 1.45NJ, (or 0.72k.N@) and not larger 
than 2.42.\VJ, (or 1.21kN@). On the other hand, the number of different en- 
ergy levels pertaining to the state 2m: =0 is N!/(N/2)!(.N/2)! or, within the 
accuracy of Stirling’s formula 2*. We find, therefore, for Zy,2 the limits 


N(log 2—1.210/T) < log Zyp, < N(log 2 —0.720/T). 


Owing to the large factor N, the change of Zy,;2 with 7 is very rapid. In 
the region log 2>1.210/7';(or T>1.86) it becomes a very large number; in 
the region log 2 <0.720/T (or T <1.046) it is a very small number. The transi- 
tion from a high probability to a low probability of the unmagnetized state 
occurs somewhere between these two limits at a point which is roughly of 
the order 

To = 1.50. (57 


The conditions are the same for the other states of magnetization of our 
crystal (2m = N—2r). Asr increases from 0 to V/2 the individual terms of the 
groups Z, defined by the Eq. (22) become smaller and smaller because their 
exponents grow in absolute value. On the other hand, the number of terms 
increases with r being equal to V!/(N—r)!r! For high temperatures this sec- 
ond effect overweighs the first: The probability of a state increases with 
r=}N—|m| and has its maximum for the unmagnetized state m=0. The 
characteristic temperature above which such a distribution exists is, nearly, 
independent of r and approximately given by the Eq. (57). For temperatures 
7’ > 7 , the most probable state and the time average derived from the sum of 
states coincide and the sum of states has a good physical meaning. The crys- 
tal will be unmagnetized or very little magnetized. The opposite is true for 
low temperatures. The negative exponents are so large that all Z; are neg- 
ligible but those pertaining to very small r. We have here the conditions 
analysed in the beginning of this section. The point Ty of Eq. (57) ts, there- 
fore, the primary Curie point of a coherent crystal: Below it the crystal is a 
state of spontaneous magnetic saturation, above it, in a, practically, unmag- 
netic state. 

Let us now return to the discussion of conditions below the Curie point. 
As we have seen, a coherent crystal is, then, in a state very near to complete 
magnetic alignment but its magnetic polarity spontaneously reverses itself 
at irregular intervals. What is the length of these intervals? Any change of 
magnetization is due to perturbations coming from the irregular thermo- 
kinetic state of the medium around the crystal and of the crystal itself. If 
there were no mutual magnetic interactions between the electronic spins, every 
electron would respond to these perturbations independently. In order to 
change its polarity, the crystal would have to pass through the extremely 
improbable state of vanishing magnetic moment and such an event could 
occur only once in an eon. In reality, there exists some measure of magnetic 
linking between the spins and they can turn over in groups so that the crystal 
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need not pass through the most improbable state. However, if the crystal is 
not small in size, the periods involved are still very long. The estimate can 
be made in a way analogous to that used by Bloch,* in a somewhat different 
case, from the magnetic linking of the spins and from the moment of momen- 
tum which the crystal has to pick up in order to reverse its magnetization. 
We call a crystal coherent if Slater’s equation applies to it as to a unit, the 
numbers V and G determining its size. If there existed coherent crystals of a 
considerable size they would be always in a state of complete magnetization 
in some direction and a field could only change this direction but not demag- 
netize them. This is contrary to observation and we are compelled to conclude 
that the coherence extends only over microscopic elements. Even a well de- 
veloped monocrystal must have a Ddlock structure. Slater’s equation and all 
the consequences contained in it which we have developed in the preceding 
sections applies only to the separate blocks. It is unimportant for us whether 
these discontinuities of a crystal are irregular (in the sense of the mosaic 
structure surmized by Darwin and Smekal) or form a regular secondary struc- 
ture, as advocated by Zwicky. The size of a few hundred atoms in each direc- 
tion G, postulated by these authors, would insure a sufficiently frequent spon- 
taneous change of magnetization. 

We arrive, in this way, at the following picture of the constitution of 
ferromagnetic materials. They are built up of microscopical blocks which are 
in a state of permanent spontaneous magnetization. The polarity of this 
magnetization is subject to frequent spontaneous changes. As far as our sim- 
plified theory goes the electronic spins might have any direction and the mag- 
netic axis could rotate freely within the block. In reality there must exist cer- 
tain directions of preference related to the main crystallographic directions be- 
cause of the interactions between spin and orbital momentum.’ The magnetic 
polarity will alternate between these preferred directions. The permanent 
magnetization of these micro-crystalline blocks is what we designate as the 
primary ferrogmagnetism. The secondary magnetization of the bulk ofa ferro- 
magnetic substance is produced by the alignment of the magnetic axes of the 
primary elements. While this picture is, by no means, new we claim to have 
demonstrated, for the first time, that the permanent magnetization of the 
blocks is not a hypothesis but a necessary consequence of the exclusion prin- 
ciple. 

11. REMARKs ON SECONDARY MAGNETIZATION 


Already in Maxwell's Treatise we find the opinion that ferromagnetic ma- 
terials are built up of permanently magnetized microcrystalline blocks. He 
believed that a magnetic field produces an alignment of the magnetic axes by 
rotating the whole block. Although this seemed little probable because of the 
enormous viscosity of these materials, I thought Maxwell’s hypothesis suffh- 
ciently interesting to have it tested experimentally. A few years ago I caused 
Mr. Yensen'® to look for a possible rotation using the Hull-Debye x-ray 


5 F. Bloch, Zeits. f. Physik 74, 295 (1932). 
® Compare F. Bloch and G. Gentile, Zeits. f. Physik 70, 395 (1931). 
0 T. D. Yensen, Phys. Rev. 32, 114 (1928). 
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method. It was not very surprising that the result was negative. According 
to our results of the preceding sections, Maxwell's assumption is unnecessary 
because the magnetic polarity of a block alternates while the block itself re- 
mains in place. In modern times the idea of permanently magnetized micro- 
crystals was used with great success by Weiss" in the case of polycrystalline 
materials. Recently, Bitter” gave a considerable amount of experimental evi- 
dence which tends to show that it applies also to monocrystals. 

The objective of this paper and its main result is to establish the primary 
ferromagnetism of the microcrystalline blocks. There exists a considerable 
literature about secondary magnetization and we limit ourselves, in this con- 
nection, to a few remarks. The agglomeration of blocks represents a large num- 
ber of similar systems side by side. We can interpret, therefore, the statistical 
averages of section 9 in a new way. They represent as well the time averages 
for a single block as the actual mean values for the agglomerate. The simpli- 
fying assumption is implied that the blocks are all equal and that the mag- 
netization has one of two opposite directions. It would be easy to drop these 
restrictions but hardly worth our while. On the contrary, we simplify still 
further by dropping the small second term of (53), i.e., by assuming that our 
blocks are completely magnetized. Eq. (52) gives then for the magnetization 
(per unit volume of the agglomerate) 

NspH 
kT 








M = nsytgh 


There are two effects left out of account in this formula, which can modify 
it sufficiently to produce a spontaneous and permanent secondary magnetiza- 
tion. The first is the field produced at a given point by the rest of the system. 
According to Lorentz," this cause can be taken care of by substituting in- 
stead of the outer field // the actual strength of field /7’=J7+aJM. In Lorentz’ 
cases a was of the order } but it is hard to say what it is equal to in our case. 
The condition for secondary permanent magnetization would then be 
(dM dll)y<0 or nNs*u? kT >1, a. The fact that monocrystals of iron and 
cobalt have little or no remanence indicates that this cause is not quite suffi- 
cient to produce permanent magnetization and points to a rather small size 
of the blocks or to a small value of a. The other effect which distorts the 
magnetization curve is the influence of mechanical stresses put forward by 
Akulov" and discussed in detail by Gans and by Becker.'® These two causes 
seem entirely sufficient to account for all the phenomena of secondary ferro- 
magnetism. We wish, only, to point out that they may give rise to a sec- 
ondary Curie point which would be observable by magnetic methods if it 
happened to be lower than the primary Curie point determined by the Eq. 


" P. Weiss, Phys. Zeits. 9, 361 (1908). 

2 F. Bitter, Phys. Rev. 38, 528 (1931); 39, 337 (1932). 

8H. A. Lorentz, Enzyklopadie der math. Wiss. V 14, section 37. 

4 N.S. Akulov, Zeits. f. Physik 64, 817 (1930). 

% R. Gans, Schriften d. Kénigsberger gelehrten Gesellschaft 8, 33 (1931). 
16 R. Becker, Zeits. f. Physik 62, 256 (1930). 
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(57). The problem of relating these theoretical Curie points to those found by 
observation is not simple and we reserve it for another communication. 


Note added in proof. It was overlooked by the writer that, in another paper 
(Leipziger Vortraege 1930, p. 67), Bloch gives the proportionality of the 
specific heat of ferromagnetic materials with 7*?. However, the agreement 
with our results is more or less accidental, as Bloch’s treatment is quite differ- 
ent and makes use of the two assumptions, pointed out at the end of our sec- 
tion 4, which we do not consider as permissible.—The experimental test of 
this law is complicated by the fact that the elements of the iron group contain 
two kinds of valency electrons, the first kind having, presumably, a negative 
interchange integral, the second a positive one. (Compare: Slater, note 6). 
The complete electronic specific heat is, then, to be represented by the sum 
of our formulas (45) and (51): c=RlaiT, ‘0, +.a250(7°/82)' 2], a, =m" /12, 
a2 = 0.208. Only 62 is related to the primary Curie point. 
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ERRATA 


The Shape of an X-ray Line 


By ARCHER Hoyt 


Department of Physics, Cornell University 
(Phys. Rev. 40, 477, 1932) 


It has been called to the author's attention that there are two errors in 
Table I of the article by the above title. 

Spencer'has been misquoted. The values published for half-maximum /a/f- 
widths are in reality half-maximum full-widths. The values for Mo Ka, and 
Cu Ka; should be reduced by a factor }. The values which shou/d appear in 


the table of half-widths at half-maximum ordinate are 


for Mo Ka, 4.8 sec. or 0.14. X.U. 
for Cu Ka, 10.8 sec. or 0.305 X.U. 


1 Spencer Phys. Rev. 38, 630 (1931). 


110 








JULY 1, 1932 


PHYSICAL REVIEW VOLUME 41 


LETTERS TO THE EDITOR 


Prompt publication of brief reports of important discoveries in physics may 
be secured by addressing them to this department. Closing dates for this depart- 
ment are, for the first issue of the month, the twenty-eighth of the preceding 
month; for the second issue, the thirteenth of the month. The Board of Editors 
does not hold itself responsible for the opinions expressed by the corres pondents. 


The Vapour Pressure of Sulphur at 50°C 


Earlier workers! have obtained widely for some time (24-48 hours) at a rate of about 








differing results for the vapour pressure of 
sulphur at 50° vis., 0.00008 and 0.00036 mm 
of mercury. We made several determinations 
using a new method. 

The apparatus is shown in Fig. 1. A is a 
drying tower containing soda-lime and cal- 
cium chloride. B is a thermostat heated and 
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Fig. 


maintained at 
50°C. C is a glass spiral containing sulphur 
(purified by distillation in vacuo) while D 
is a Pyrex tube heated electrically to 300°C. 
The water jacket E cools the issuing gases 
which are then absorbed in N/100 iodine 
solution in wash bottle F. Bottles G, G con- 
tain potassium iodide solution to trap any 


controlled electrically and 


iodine carried over. 

The procedure is as follows. Air is drawn 
through the apparatus by means of a filter 
pump fitted with constant leak. The air passes 
first through a gas meter and drying tower 
and then over the sulphur in the spiral which 
is at 50°C. The sulphur is oxidised to sulphur 
dioxide in the Pyrex furnace (owing to the 
very high partial pressure of oxygen this takes 
place alinost quantitatively) and the issuing 
sulphur dioxide after cooling is absorbed in 
the N/100 iodine solution. A\iter passing air 


10 litres/hour the wash bottles are detached 
and the iodine in them titrated with a sodium 
thiosulphate solution previously standardized 
against the iodine. 

Thus we find the loss in iodine and hence 
the weight of sulphur dioxide formed. Know- 
ing from the gas meter the volume of air 
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which has passed and also the barometric 
pressure and the temperature at the meter, we 
can readily calculate the vapour pressure of 
sulphur at 50°C. 

We obtained as the mean value 0.000155 
mm Ig assuming all the sulphur present as 
Ss. 

In conclusion we wish to offer our best 
thanks to Mr. A. Ritchie, B.Sc for his in- 
valuable help and advice. 

R. R. H. Brown 
J. J. Muir 
Edinburgh University, 
Chemistry Department, 
April, 1932. 


1 Ruff and Graf, Z. Anorg. Chem. 58 210 
(1908); Greuner, J.A.C.S. 1399 (1907); West 
and Menzirs., J. Phys. Chem. 33, 1887 
(1929), 


Variation of the Cosmic Rays with Latitude 


Definite differences in the intensity of the 
cosmic rays at different latitudes are shown 
by our measurements, which have ranged 


from 47° north to 46° south. As far as they 
have gone, these measurements indicate a 
uniform variation with latitude, showing a 
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minimum at or near the equator, and increas- 
ing intensity toward the north and south 
poles. At sea level, the difference between the 
intensity at latitude 45° and 0° is roughly 16 
percent, whereas at an elevation of 9000 feet 
the difference is about 23 percent. This would 
indicate that it is the least penetrating part 
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readings being made in such a way that the 
effect of any insulation leak is eliminated. 
To shield from the effects of the local radia- 
tion, a spherical shell of 2.5 cm copper and 
two spherical shells each of 2.5 cm lead, mak- 
ing a total of 7.5 cm of lead and copper, are 
used. The effect of the local gamma-rays 


TABLE I. Cosmic ray intensity at different localities 
(Ions per cc per sec. through 5 cm Pb, 2.5 em Cu and 0.5 cm Fe) 


Location 


Lat. Long. Elev. Barom. Ice Th Date 
1 Mt. Evans 40°N 106°W 14,200ft 17.6lin 6.88 ions 0.57 9/31 
2 Summit Lake 40 N 106 W 12,700 18.70 5.84 0.34 9 31 
3 Denver 40 N 105 W 5300 24.8 2.93 —_ 9/31 
4 Jungfraujoch 47 N 6E 11,400 19.70 5.08 0.51 10/31 
5 Haleakala 21 N 156 W 9300 21.47 3.35+0.05 0.60 4/32 
6 Idlewild 21 N 156 W 4200 25.99 2.40+0.05 0.37 4/32 
7 Honolulu 21N 158 W 70 30.09 1.89+0.02 0.11 4/32 
8 S. S. Aorangi 45S 173 W 60 29.65 1.83+0.05 0.32 4/32 
9 Southern Alps 44S 170 E 6700 23.69 3.39+0.05 0.22 4/32 
10 Southern Alps 44S 170 E 3900 26.10 2.70+0.04 0.21 4/32 
11 Dunedin 46S 170 E 80 30.08 2.16+0.03 0.11 4/32 

41S 175 E 400 29.85 2 0.12 5 32 





12 Wellington 


.16+0.03 





of the cosmic rays which varies most rapidly 
with latitude. No significant variations with 
longitude have been noted. 

Our present experiments are being made 
with a 10 cm steel sphere filled with argon at 
30 atmospheres. The ionization current is 
measured by a Lindemann electrometer, the 


traversing this shield is determined by com- 
paring the ionization with and without the 
outer shell of 2.5 cm lead. All the readings are 
referred to the ionization by the gamma-rays 
from a capsule of 1.3 mg of radium placed at 
1 meter from the center of the ionization 
chamber as a standard. 
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In Table I are listed the results of our vari- 
ous measurements. The column marked Jc 
represents the ionization in air at atmospheric 
pressure expressed in ions per cc per second 
due to the cosmic rays traversing our 7.5 cm 
filter. J, represents the ionization due to the 
local rays passing through this filter expressed 
in the same units. Measurements 1 to 4 were 
made with different apparatus, using com- 
pressed air instead of compressed argon in the 
chamber, and lead and copper shields of 
different shape. These results are thus not 
strictly comparable with the later ones, 
though the difference is probably not more 
than 2 percent. 

In Fig. 1 these results are shown graphi- 
cally, and are compared with those of Millikan 
and his collaborators (Phys. Rev. 37, 235 
1931) as shown in the dotted line. It will be 
noted that for corresponding barometric 
pressures our measurements show consistently 
more intense radiation at the higher latitudes. 
The high altitude measurements in Hawaii 
are of especial interest as indicating that the 
difference with latitude is greater for the 
less penetrating rays. 

Millikan’s measurements were made at 
latitudes varying from 34°N to 59°N, most 
of them being at about 40°N. Our measure- 
ments near this latitude are in satisfactory 
agreement with his considering the slight 
differences in experimental arrangement. We 
do not however confirm the independence of 
the intensity of latitude which he reports. 

The variation of the cosmic rays with 
atitude which these measurements show is of 
_ust the kind to be expected if the rays con- 
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sist of electrically charged particles which are 
deflected by the earth’s magnetic field, the 
less penetrating rays being the more strongly 
affected. 


This letter is the first report of an extensive 
program involving similar measurements by 
many physicists in widely distributed parts 
of the world. The design and construction 
of the instruments are due chiefly to Professor 
R. D. Bennett of Massachusetts Institute of 
Technology, with L. N. Ridenour, Dr. J. A. 
Hopfield and Dr. E. O. Wollan. Professor J. C. 
Stearns of the University of Denver with J. A. 
Longman, L. N. Ridenour and W. Overbeck 
made the measurements at Denver and Mt. 
Evans. Dr. Marcel Schein and Dr. Bernhard 
Frey of the University of Ziirich cooperated 
in the measurements at the Jungfraujoch. 
Professor Harry Kirkpatrick of the University 
of Hawaii made the local arrangements in 
Hawaii, and with Professor W. H. Eller and 
Mr. I. Miyake assisted in the measurements. 
In New Zealand, Professor P. W. Burbidge 
of Auckland University College organized the 
various expeditions, and with Dr. C. M. 
Focken of the University of Otago helped with 
the observations. Valuable assistance has 
also been given by A. A. Compton and others. 
Without the cordial cooperation of the vari- 
ous universities where the work has taken us, 
the measurements would have been much 
more difficult. 

ARTHUR H. Compton 

University of Chicago, 

The Tasman Sea, 
May 7, 1932. 


On the Production of Ultra-Short Electromagnetic Waves 


In a note published in this journal, Mr. G. 
Potapenko gave an account of his interesting 
investigations in the field of the ultra-short 
electromagnetic waves. 

The tubes that the author used in his ex- 
periments had pure tungsten filaments, as it 
was impossible to observe ultra-short waves 
with tubes having other kinds of filaments.! 

As a result of some experiments I carried 
out recently, I found that, under particular 
conditions, oscillations of very high frequency 
may be obtained by using tubes with a heater 
element and a cathode coated with oxides. 

Different tubes have been tested, namely, 
tubes of the type: —27 R.C.A. (American); 


E. 424 Philips (Dutch); R.E.N.S. 1104 Tele- 
funken (German). 

Although the common Barkhausen-Kurz 
arrangement may be used, that of Tank and 
Schiltneck? was found to be more suitable, 
the only modification being to connect the 
cathode to the negative terminal of the fila- 
ment. The heating current must have a value 
which is only 80 percent —90 percent of its 
normal rated value. 

Although unnecessary if the latter experi- 
mental arrangement is used, a slight positive 


1 Phys. Rev. 39, 630 (1932). 
2 Helv. Phys. Acta. 1, 100 (1928). 
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plate potential (0.5 to 2 volts) increases the 
amplitude of oscillations. Under these con- 
ditions a plate current is present even when 
the tube is not oscillating; this current being 
obviously increased when oscillations set up. 
This increase of current may vary (in ac- 
cordance with the particular tube, etc.) be- 
tween 10 and 200 microamperes. As the grid 
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potential runs between 20 and 100 volts, a 
band of wave-lengths between 180 and 50 
cm can be obtained. 
AMEDEO GIACOMINI 
Istituto di Fisica della Regia Universita, 
Milano, Italy, 
June 1, 1932. 


Predissociation in Nitrogen 


Predissociation in nitrogen is best illus- 
trated by the first positive bands B,;—Ag, 
By Ajo, Bys Ay and by the bands By;—As, 
By-Ag and Bys—Ajo.! On small dispersion 
spectograms most of the bands of the first 
positive group possess three fairly strong 
heads, which correspond to transitions from 
the three components of the *z initial state 
of these bands to the *S state. In the above- 
mentioned bands however, two of the heads 
are almost completely missing. About a year 
ago, I reported the excitation of these missing 
heads in mixtures of nitrogen and oxygen,? 
in which the concentration of nitrogen was 
relatively small compared with that of the 
oxygen. At that time, the explanation of this 
effect was proposed to be the prevention of 
the normally occurring predissociation by 
electric fields, which were introduced during 
collisions. There is an obvious alternative ex- 
planation in which electric fields play no 
direct role, namely, the effect of increasing 
the number of molecules in the normally pre- 
dissociating levels by collisions which give 
rise to transition from the non-predissociat ing 
component of the *z state to the two pre- 
dissociating ones. The energy differences be- 
tween the components of the *w state are 
very small and such transfers should take 
place without difficulty. The effect of pressure 
on predissociation has been discussed for 
AIH by Stenvinkel® and others, and this effect 
in nitrogen is in some respects, similar to the 
one in AIH, There is a difference, however, 
for in nitrogen the collision transfers are be- 
tween components of an electronic state, 
whereas in AIH they are between rotational 
levels. 

The occurrence of this effect in nitrogen 
suggests a similar explanation for the en- 


hancement of predissociation in iodine by 
argon.’ Mulliken’s assignments of *xo*+ to the 
upper state of the visible halogen bands, 
makes this state a triplet state and hence pro- 
vides the two other components, *x; and *z2, 
which are necessary in order that the iodine 
problem be similar to the one in Ne. In 
iodine, the two levels *x; and *z: may them- 
selves be unstable or they may combine with 
a third unstable level as in nitrogen. 

Before the appearance of Van Vleck’s® 
paper on the theory of the magnetic quench- 
ing of iodine fluorescence, I had proposed to 
myself the problem of studying the effect of 
magnetic fields on the predissociation in nitro- 
gen, in order to observe the effect of magnetic 
fields on collisions. It was thought at the 
time that the magnetic field might be able to 
enhance collision transitions from the non- 
predissociating *7 state to the two predissoci- 
ating ones, and thus be able to quench the 
fluorescence. I am not aware of any theoretical 
reason for believing that this explanation has 
any value, but in view of the effect of col- 
lisions on predissociation both in iodine and 
in nitrogen, it seems worth while to suggest 
that the effect of magnetic fields on col- 
lisions may be of some significance. 

JoserH KAPLAN 

University of California 

at Los Angeles, 
June 3, 1932. 


1 J. Kaplan, Phys. Rev. 38, 373 (1931). 

2 J. Kaplan, Phys. Rev. 38, 1079 (1931). 

3 Stenvinkel, 
(1930). 

4 Turner, Phys. Rev. 38, 574 (1931). 

5 Van Vleck, Phys. Rev. 40, 544 (1932). 


Zeits. f. Physik. 62, 201 
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The Isotopic Weight of H? 


The mass of neutral H? has been measured 
as 2.01353+0.000064 by comparison of the 
mass of (H'H!H?)+ with Het on fourteen 
spectra obtained with the mass-spectrograph 
recently described.! For the purpose of these 
measurements the mass of helium was taken 
as 4.00216, that of hydrogen 1.00778, and 
for the electron 0.00055 units on the O = 16 
scale.2 The equivalent packing fraction of H? 
is 67.6 parts in 10,000. The energy of binding 
corresponds to 0.00203 mass units or ap- 
proximately 2X10° electron-volts if the 
nucleus is composed of three entities, two 
protons and one electron. It is interesting to 
speculate that if the H? nucleus is composed 
of one proton and one Chadwick neurton’ of 
mass 1.0067, then the binding energy would 
be 9.5 X 10° electron-volts and the H? nucleus 
would be loosely bound compared to He and 
other light nuclei. Professor Urey has pointed 
out to the author that for a stable configura- 
tion the mass of a free neutron must be greater 
than 1.00575. 

The probable error of a single determina- 
tion was +0.00031. The observations were 
weighted from one to three depending on the 
definition of the lines of specific spectra and 
the internal consistency of the mensuration 
of the position of those lines. The sum of the 
weight factors was twenty-three. 

(3H')* and He* provided the dispersion 
measurements for these spectra. The presence 
of (H'H?)* can only introduce a_ negligible 
correction. Earlier unreported work on the 
carbon and oxygen hydrides has shown that 
the mass scale is linear to within one part in 
ten thousand over the first eight centimeters 
of the plate, which includes the region used 
in this investigation. 

Lines on the mass-spectrograph plates cor- 
responding to a mass of 4.0285 were attributed 
to (H'H'H?)* because; (1) no lines of com- 
parable intensity appeared in this position 
when ordinary commercial hydrogen was used 
in the discharge tube; (2) the probability of 
the formation of (2H?)* was much less than 
the probability of (H'H'H?)* under the con- 


ditions existing in the discharge tube, i.e., 
(3H')* was present from 100 to 5 percent 
of (2H')*; (3) the mass is less than the mass 
of (4H')* by an amount outside the limit of 
error. 

Two samples of enriched hydrogen were 
used. Dr. Bleakney gave the writer a sample 
of the hydrogen he had examined, and Dr. 
Brickwedde and Professor Urey® sent a sample 
of hydrogen with a slightly higher concentra- 
tion of H?. 

The high value of the mass of H? makes it 
very difficult at this time to resolve the two 
doublets (3H')*, (H'H?)*, and (2H')*, H?*, 
It may be possible to resolve these doublets 
when more concentrated samples of hydrogen 
are available so that the lines of the respective 
ions can be obtained of comparable intensity. 
The doublet (H'H'H?)*, (2H?)* will be diffi- 
cult to resolve in any case as the respective 
masses differ by only one part in 2000. 

The author is particularly indebted to Dr. 
Brickwedde who made this investigation pos- 
sible by his preparation of the samples of en- 
riched hydrogen at the U. S. Bureau of 
Standards, to Professor Urey and Dr. Murphy 
of Columbia University who placed a tested 
sample at the, disposal of the author and 
to Dr. Bramley and Dr. Bleakney for their 
interest and helpful discussions. 

KENNETH T. BAINBRIDGE 

Bartol Research Foundation of the Frank- 

lin Institute, Swarthmore, Pa., 
June 21, 1932. 


1K. T. Bainbridge, Phys. Rev. 40, 130 
(1932). 

2F. W. Aston, Proc. Roy. Soc. 115, 487 
(1927). 

3 J. Chadwick, Proc. Roy. Soc. 136, 702 
(1932). 


4W. Bleakney, Phys. Rev. 39, 536 (1932); 
Phys. Rev. 41, July 1, 1932. 

5 H.C. Urey, F. G. Brickwedde, and G. M. 
Murphy, Phys. Rev. 40, 1 (1932); 39, 164 
(1932). 
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BOOK REVIEWS 


Photoelectric Phenomena. ARTHUR LLEWELLYN HUGHES AND LEE ALVIN DUBRIDGE. 
Pp. 531. McGraw-Hill Book Company, Inc., New York, 1932. Price $5.00. 


The preparation of an up-to-date text on “Photoelectric Phenomena” constitutes an al- 
most hopeless undertaking. So many and varied are the researches in this field published in 
current technical journals that to digest, organize and include the subject matter in an orderly 
discussion would be next to impossible. The task is doubly difficult due to the fact that so little 
is known concerning the nature of the photoelectric effect and so many theories have been ad- 
vanced to explain the oftimes conflicting results obtained by faulty experimental technique. 

It is for these very same reasons, however, that there has been felt so keenly the need of a 
treatise coordinating the data so far obtained; and it is a source of considerable satisfaction to 
investigators in the field of photoelectricity to discover that “Electrical Phenomena” by Hughes 
and DuBridge so admirably meets these needs. The fact that several rather important papers 
published within the year seem to have been overlooked in its preparation, particularly in the 
field of photoconductivity, and the almost reportorial presentation of the data resorted to by 
the authors with little or no attempt to pass judgment on the relative merits of the research 
are to be expected rather than criticized in view of what has already been said. 

“Photoelectric Phenomena” is a comprehensive well written, carefully edited and at- 
tractively published volume, far eclipsing any other volume on the subject published in the 
English language and perhaps being second to none emanating from the continent. 

Defining a “photoelectric effect as one in which radiation causes a separation of electron 
from atoms so that an electric current may be caused to flow,” the authors set out “to give a 
concise yet comprehensive survey of one of the many fields of physics.” And within the fourteen 
chapters comprising the book they accomplish just that. The text, after a brief introductory 
chapter, consists of four chapters on the emission of electrons from metals, one chapter on the 
theories of photoelectric emission, a chapter each on what they have called “volume” photo- 
electric effects, viz., ionization of gases by light, and photoconductivity, two chapters on the 
photovoltaic effect and photoelectric effect in liquids and dielectrics, a chapter on the photo- 
electric effect of x-rays and y-rays, a chapter each on photoelectric technique and applications 
and ends with a discussion of miscellaneous phenomena. The appendix includes a convenient 
table of physical constants, important photoelectric energy relations, a periodic table of the 
elements and a summary of illuminating engineering units and nomenclature. A very complete 
running bibliography is included in the form of footnotes with a name index and subject index 
at the end of the text. The book is printed in easily readable type and the sketches and diagrams 
are many and clear. 

Unlike previous publications in English, too much valuable space is not devoted to a long 
historical account of the pioneering work or to the multitudinous applications of photocells. 
The text is essentially up-to-date and stresses important problems of the present. The order of 
presentation, chapter by chapter, is: first, ample citation of and reference to early related 
work, secondly, a complete resumé of all researches since that time, and finally a summary and 
discussion which generally includes a concise outline of the problems still unanswered. After 
presenting the researches of others with such reportorial excellence, this discussion of the out- 
standing problems is especially commendable. Incidently, it reflects the thorough understand- 
ing of the field possessed by the authors. 

One of the most praiseworthy features of the presentation is the differentiation bet ween 
“spectral selectivity” and “polarization selectivity.” Yet, after keeping these phenomena sepa- 
rated, as they should be, the authors, without justification perhaps, advance the doubtful sug- 
gestion that the latter is always associated with the former. Unfortunately, at the very begin- 
ning of this chapter a confusing error crept into the definition of polarization selectivity. 

There are undoubtedly those who will criticize the order of presentation in the chapter on 
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the selective photoelectric effect and that on photoconductivity, just as they will the arbitrary 
classification of photoelectric phenomena as “surface photoelectric effects,” “volume photo- 
electric effects” and “photovoltaic effects.” This might follow if the reader is accustomed to the 
style of Gudden and Pohl in their volume bearing the same title as the one under discussion. In 
and of itself, however, the subject matter in “Photoelectric Phenomena,” by Hughes and Du- 
Bridge will be found to be wisely classified, and coherently, competently and clearly presented. 
Even the chapter on “Theories of the Photoelectric Effect,” though brief and undertaken 
apologetically by the authors, is entirely adequate for such a volume. Moreover, the discussion 
on photoelectric technique is pertinent and practicable, and can well be accepted as reflecting 
the experience of careful research workers in the field of photoelectricity. 

A. R. OLPIN 

Bell Telephone Laboratories, Inc., 

New York, N. Y. 


Cartesian Tensors. HAROLD JEFFREYS. Pp. 93, The Macmillan Company, New York 
1931. Price—-$1.75. 


According to the author's preface this little book is intended as a frank attempt to popu- 
larize the tensor methods for systems of rectangular axes, even to the extent of eliminating the 
usual extension development of vector analysis. Even within the realm of classical physics, to 
say nothing of relativity, one can find good reason for sympathizing with this plan. As examples, 
one thinks of the inertial ellipsoid of a rigid body and the polarization ellipsoid of a molecule. 
In these, and in other similar instances, the tensor notation seems much preferable to the 
Gibbsian dyadics. 

The main emphasis in the book is on applications of tensors in classical physics. The fact 
that only twenty short pages are needed for the development of the tensor calculus itself 
should recommend it. The applications are made to dynamics of particles and rigid bodies, 
force systems, electrodynamics (potential theory), elasticity, and hydrodynamics. The discus- 
sions of the individual topics are too condensed to give more than a taste but taken together 
they give the student a valuable summary. The explanations are stated and easily followed. 

E. L. HILt 


University of Minnesota 


Eenige Onderzoekingen over Paramagnetisme. E. C. Wiersma. Pp. 122. Paramagnetische 
Eigenschaften von Salzen. C. J. Gorter. Pp. 112. 


It is perhaps well to call attention in these columns to these two interesting Leiden disser- 
tations devoted to topics in paramagnetism. Both dissertations are primarily the outgrowth of 
original experimental investigations, but at the same time enter into theoretical considerations 
with considerable detail. It is to be particularly commended that attention is focused on the 
newer quantum-mechanical developments rather than on obsolete classical models, as is too 
often done. In fact, over half of Wiersma’s monograph is devoted to theory. The remainder is 
devoted to his experimental methods and results. The latter include the behavior of O. and NO 
at low temperatures, which is of importance in connection with theory for gases, and some 
interesting cryomagnetic anomalies in solid iron chloride, in which the susceptibility depends 
upon field strength in a peculiar way. Gorter’s monograph is particularly valuable for the care 
with which it documents the existing determinations of susceptibilities and magneton numbers 
in the rare earth and iron groups. It is probably the best available compendium of such data. 
Gorter’s own measurements have been indispensable in this field. They have corrected anoma- 
lies reported by other observers for hydrated nickel sulphate. Also his observations on Pr and 
Nd salts at low temperatures have enabled theoretical physists to estimate quantitatively the 
crystalline field. 

J. H. Van VLECK 
University of Wisconsin 
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“Einfiihrung in die Mechanik und Akustik” von R. W. Pont, Zweite Auflage, Pg. 251 
tviil, Figs. 440. Julius Springer, 1931. Price RM 15.80. 

This volume constitutes an important contribution to the presentation and the elucidation 
of the introductory principles of mechanics. It is well organized and the author's style of a 
simple, plain unfolding of the facts of nature as revealed by experiment is a very pleasing one 
for the reader. The illustrations are numerous, 440 altogether, and their quality of reproduction 
is of a very high order. The device of silhouetting figures is frequently used and with good effect, 
and the experiments described exhibit methods of highly refined technique. 

The first eight chapters, consisting of 115 pages, treat the mechanics of particles and of 
rigid bodies up to and including rotational motion and accelerated reference systems. The next 
two chapters of 54 pages deal with the mechanics of fluids, and these complete the mechanics 
division of the book. The division on acoustics consists of 75 pages and is divided into two chap- 
ters, the first of which deals with vibration theory and the second with the theory of wave mo- 
tion and radiation. 

As a groundwork for mechanics force is taken as a fundamental idea, and in pleasing con- 
trast with so many texts which adopt this method the definition of force measurement is given 
at the outset. It is true that the author introduces the term mass simultaneously with that of 
force and also a convention for its measurement. Newton's second law is then an empirical 
relationship between force and mass as defined. 

This mode of development of fundamental ideas is a possible one and is one that is very 
commonly used. It has certain advantages and its desirability is a matter for personal opinion. 
On the other hand the author expresses views to which the reviewer does not subscribe. The 
identity between gravitational and inertial mass is stated to be one of the curiosities of classical 
mechanics, a curiosity which was first made clear by Einstein’s general theory of relativity. 
This view has also been expressed by others. However, a distinction between gravitational and 
inertial mass is not inherent in Newtonian mechanicsalt hough it does seem to be inherent in the 
mode of development which regards force as a fundamental idea. The trucks shown in Fig. 46 
would afford an excellent means of introducing the mass concept without referring to forces an 
without introducing such vague terms as inertia or heaviness to account for it. One of the 
dangers in the use of force as a fundamental idea is that of attributing to it a “causal” signifi- 
cance rather than a “descriptive” one. 


J. W. CAMPBELL 
University of Alberta 








CURRENT LITERATURE OF PHYSICS 


Extracts from the Tables of Contents of Some 
of our Contemporaries 


NATURE 
VoL. 129, No. 3265 May 28, 1932 
Letters to the Editor: 
Ionization at High Gas Pressures. By R. M. Sievert... . 2... 2020 ce 792 
Atomic Weight of Fluorine. By H. S. Patterson and W. Cawood............ i 794 
Vibrational Specific Heat of Carbon Dioxide. By H. O. Kneser. . 5 stanton a Rude ats 797 
Types of Iridescent Clouds. By H. S. Wynne-Edwards and T. S. Dymond. rieaed eoeen’ 798 
Winncel Sais. Thy CG. Tso son isk ces nce canceceeeevenaceevens sce seces 798 
Neutron Hypothesis. By D. Iwanenko........ 2.0... 789 
NATURE 

VoL. 129, No. 3266 June 4, 1932 
Letters to the Editor: 
Nuclear Spin of Phosphorus from the Band Spectrum. By F. A. Jenkins and Muriel 

BE ir ig iG t nha a tS hen aN ON 5K te EME ME NE A Caan 829 
Disintegration of Atomic Nuclei. By H. S. Allen... .. eatersaee , tevecsese 88 
Oscillations of the Methane Molecule. By S. Bhagavantam. NCP ere wieptteoscéooen fe 
Errors in Thermal Measurements. By W. E. Garner.................... - , 832 
Free Ethyl. By T. G. Pearson, P. L. Robinson, and E. M. Stoddart................. 832 
Electrostatic Explanation of the Phenomenon of Flotation. By Ian W. Wark ; .. 833 
Ce OE TOR, Bey BI Oe kono do reiksiersaceaewaewen - oka. bthaaut a0 
Isotopic Constitution of Lead. By H. Schiiler and E. Gwynne Jones bu Ae ea stvuvecce ‘ae 


DIE NATURWISSENSCHAFTEN 
VoL. 20, Nos. 22/24 May 27, 1932 
Proceedings and Annual Reports of the Kaiser Wilhelm-Gesellschaft zur F érderung der 
Wissenschaften 
New Electrochemical Processes of Radiochemistry. By Otto Erbacher................ 390 
On the Theory of Atomic Disintegration by Resonance. By H. Kallmann.... piaeem ware 


ANNALEN DER PHYSIK 
VoL. 13, No. 4 May 4, 1932 
Faraday Effect in Fluid Mixtures. By K. Scharf...................... seem aan 
Measurements between Two Perfectly Conducting Planes in the Radiation Field of a 
Vertical Antenna Excited to its Fundamental Vibration. By L. Bergmann and W. 


aD eae a ai te ha Sala elas a arein dichila sh ikea an atinkat es edatide Ske abate eRe RR 409 
Theoretical Considerations about Cosmic Radiation. By W. Heisenberg... . ven a 
Electron Diffraction and Molecular Structure. Il. By R. Wierl....... 453 


On Changes in the Striking Place of Long Electric Sparks without Change of Spark Gap. 


Fe ais ak Rho eee AE TAR ROE ORD RHEE RE REE SeNe Re paworeans wee 
COMPTES RENDUS 

VoL. 194, No. 20 May 17, 1932 

Effect of Hysteresis in Heating by an Oscillating Magnetic Field. By André Blondel..... 1700 

Diameter of CO Molecule. By E. Mathias, W. J. Bijleveld and Ph. P. Grigg......... 1708 


Variation of Drag at Low Velocities under the Influence of the Compressibility. By 
NE IN oie iss a sie gee Weed w MaKe wig e eRe araus aCe A Cerne md ee ae 1720 








On the Role Played by the Nature of the Electrodes in the Conductivity of Liquid Semi- 
conductors. By J. Sambussy 

On the Variable Paramagnetism of Crystallized Iron Perchloride and the Constant 
Paramagnetism of the Fe? Cl® Molecule in the Gaseous State. Constitution of the 
Molecule. By André Lallemand 

Double Salts, Complex Salts and Circular Dichroism. By J. P. Mathieu 

Stroboscopic Method for the Measurement of Electric Double Refraction. By Réne 
Lucas and Marcel Schwob 

Actinium Uranium Ratio in an Old Uraninite, the cleveite of Aust-Agder (Norway). 
By Ellen Gleditsch and Sverre Klemetsen 

Radioactive Phenomena of the Second Order and of Artificial Origin. By G. Reboul 

Effect of Substiuitions on the Vibration Frequencies of Ethylene Compounds. Procedure 
for the Classification of Radicals. By M. Bourguel 


COMPTES RENDUS 
VoL. 194, No. 21 May 23, 1932 
On the Temperature of a Piezoelectric Crystal as a Function of its State of Oscillation. 

By Armand de Gramont and Daniel Beretzki 
Application of the Method of the Variation of Constants to the Problem of Bodies 

with Variable Masses. By M. Mendes 1794 
On Some Problems of Elastic Vibrations. By V. Smirnoff and S. Soboleff............ 1797 
Discussion of a Nernst Bridge for High Precision. By Jacqueline Hadamard 1799 
Influence of Solvents and of the Temperature on the Optical Activity and Rotary Dis- 

persion of Active Substances. By R. Lucas and D. Biquard 1805 
Role of Multiple Reflection in the Magneto-Optic Kerr Effects of Thin Sheets of Iron. 

By Marcel Cau 1807 
Influence of Diffused Light on Photoelectric Measurements. By Trajan D. Gheorghiu... 1810 
Passage of a-Rays through Air and the Theory of Bethe. By G. Mano 1813 
Transparency of the Pure Atmosphere. By J. Duclaux and M. Hugon.. 1842 
Transparency of Air. By P. Chofardet 








